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a  Avail-Able  information  relating  to  cratering  and 
penetration  in  metallic  targets  has  been  gathered.  Impact 
of  compact  particles ,  microparticles ,  and  rods  at  normal 
and  oblique  incidence  on  quasi- infinite  targets  and  on  thin 
targets,  multiple  spaced  targets,  and  shielded  targets  is 

r\  K  C. 

considered.  Experimental  data  -irs-^ presented  in  tabular  aw.' 
graphical  form  for  ready  references.  Available  theories  and 
semi-empirical  theories,  as  well  as  empirical  correlation 
equations  are  summarised  nn  nnd  ari^rrTnpirrd 

with  each  other  and  with  the  experimental  data.  For  normal 
impact  ch  quasi- infinite  targets  -two-A empirics l  correlation 
express idus  are  deduced  which  are  more  generally  applicable 
than  previous  expressions.  A  qualitative  description  of  the 
cratering  process  is  given,  and  realistic  regions  of  impact 
are 'det'icdd.  Recommendations  for  future  ext:  drimental  and 
theoretical  work  are  made. 
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SECTION  I 


INTRODUCTION 


Cratering  and  penetration  of  targets  by 
projectiles  travelling  at  velocities  at  which  the  pro¬ 
jectiles  suffer  severe  deformation  or  breakup  has  recently 
received  a  great  deal  of  attention,  due  to  the  importance 
of  this  phenomenon  in  the  protection  and  lethality  of 
space  vehicles  and  ballistic  missiles „ 

A  great  deal  of  experimental  data  has  been 
published  by  a  number  of  laboratories,  and  considerable 
theoretical  work  has  appeared  during  the  last  few  years . 
Limited  comparisons  of  data  from  one  laboratory  with  that 
from  another,  and  of  experimental  data  with  theory  have 
shown  some  important  discrepancies .  Considerable  confusion 
exists  because,  except  in  a  few  instances,  the  data  from 
one  laboratory  are  not  directly  comparable  with  those  from 
another,  since  identical  materials,  projectile  shapes  or 
velocity  ranges  were  not  used.  As  a  result,  several  more 
or  less  contradictory  empirical  expressions  have  been  pro¬ 
duced  at  different  laboratories  to  fit.  their  own  data.. 

Only  the  lower  end  of  the  velocity  range  of  interest  in 
space  applications  has  been  explored  experimentally  in  any 
detail.  Extrapolations  of  the  various  empirical  data  fits 
to  higher  velocities  leads  to  large  discrepancies. 

Experimental-theoretical  comparisons  have  been 
hampered  by  the  fact  that  on  the  one  hand,  target  and  pro¬ 
jectile  material  properties  appearing  in  the  theories  have 
not  been  measured  or  reported  by  many  laboratories  On  the 
other  hand  most,  of  the  theories  are  simplified  to  the  point 
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where  they  may  be  applicable  to  only  very  limited  ranges  of 
velocity  and  material  properties,  and  the  range  of  applica¬ 
bility  of  these  theories  is  not  clear „ 

In  an  attempt  to  assess  the  present  status  of 
experimental  and  theoretical  knowledge  of  impact  phenomena, 
experimental  data  published  prior  to  31  March  1961  was 
gathered  and  checked  as  far  as  possible  for  consistency, 
and  compared  where  possible  with  the  available  theories. 

An  attempt  was  also  made  to  correlate  the  data  on  the  basis 
of  simple  empirical  expressions. 

This  report  is  limited  to  consideration  of  impact 
into  metallic  targets  under  conditions  such  that  the  projec¬ 
tile  suffers  severe  deformation  or  breakup  in  the  cratering 
process.  Data  on  shaped  charge  jet  penetration  has  been 
excluded  since  this  has  been  covered  adequately  elsewhere.. 
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SECTION  II 


EXPERIMENTAL  DATA 


During  the  data  acquisition  phase,  thirty-nine 
groups  active  in  high  velocity  impact  or  in  developing  nro- 
jection  techniques  were  either  visited  or  contacted  by 
telephone.  Not  every  group  had  performed  cratering  experi¬ 
ments  o  Table  I  lists  the  laboratories  in  possession  of 
experimental  facilities,  together  with  an  indication  of  the 
type  of  work  reported.  The  bibliography  lists  most  of  the 
important  pxiblications  containing  information  on  cratering 
generated  at  these  laboratories. 

Several  different  target  configurations  have  been 
used.  Most  of  the  data  related  to  configurations  shown  in 
Fig  1.1  They  may  be  classified  as: 

(a)  Normal  impact  on  quasi- infinite  targets 

(b)  Oblique  impact  on  quasi-infinite  targets 

(c)  Rod  impact  on  quasi-infinite  targets 

(d)  Impact  on  single  thin  targets 

(e)  Impact  on  quasi-inf inite  targets  protected  by 
a  thin  shield 

(f)  Impact  on  multiple  spaced  thin  targets 

An  attempt  has  been  made  to  extract  the  original 
measured  quantities  for  each  experimental  firing.,  In  some 
cases,  published  reports  contained  tabulations  of  measured 
quantities..  In  many  cases,  the  data  was  plotted  in  non-di¬ 
mensional  form,  and  the  extraction  of  the  measured  quant it ie 


required  considerable  care.  In  many  cases.,  additional  infor¬ 
mation  was  requested,  and  received  from  the  laboratories 
concerned.  In  a  few  cases,  some  difficulty  was  experienced  in 
extracting  measured  quantities  due  to  ambiguities,  or  exces¬ 
sively  small  scale  data  plots,  and  as  a  consequence  some  un¬ 
certainty  is  connected  with  the  corresponding  data. 

Material  properties  for  both  projectile  and  target 
materials  were  quoted  in  the  publications  in  only  a  few  in¬ 
stances.  Some  of  the  data  on  material  properties  which  were 
quoted,  were  handbook  values.  In  only  a  few  instances  were 
material  properties  measured  on  the  actual  test  specimens. 
Where  adequate  data  were  missing,  additional  information  was 
requested  from  the  laboratories  concerned.  In  some  instances 
where  no  information  on  material  properties  was  supplied, 
average  handbook  values  corresponding  to  the  material  speci¬ 
fication  were  used.  In  a  few  instances  no  detailed  material 
specifications  were  supplied,  or  the  material  specifications 
were  insufficient  to  allow  reasonably  reliable  handbook  values 
to  be  chosen. 

The  data  are  tabulated  in  the  Addendum. 

By  far  the  largest  proportion  of  firings  have  been 
made  into  targets  of  quasi-infinite  thickness  at  normal  inci¬ 
dence,  over  1700  data  points  having  been  reported  for  this 
conf iguraion.  The  target  material-projectile  material  com¬ 
binations  which  have  been  used  for  this  configuration  are 
shown  in  matrix  form  in  Table  II.  Inserted  are  the  velocity 
ranges  in  feet  per  second  used  at  the  various  laboratories. 
Data  for  those  material  combinations  for  which  sufficient 
data  points  exist  are  plotted  in  Figs. 1. 2  to  i. 2b.  Cratering 
in  quasi-infinite  targets  by  compact  projectiles  at  normal 
incidence  is  discussed  in  Section  V 
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Only  relatively  few  firings  have  been  conducted 
using  other  projectile- target  configurations.  Presentation 
of  data  and  discussion  relating  to  these  cases  is  deferred 
to  later  sections. 
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SECTION  III 


EXISTING  CRATERING  AND  PENETRATION  THEORIES 


Several  theoretical  analyses  of  high  velocity 
cratering  and  penetration  have  been  made.  These  may  be  con¬ 
veniently  divided  into  four  groups  in  accordance  with  the 
simplifying  assumptions  involved: 

(a)  Rigid  Projectile,  i.e.  the  projectile  is  undeformed 
during  penetration. 

(b)  Hydrodynamic .  i.e.  the  projectile  has  zero  shear 
strength.  Strength  of  the  target  may  be  neglected, 
or  accounted  for  by  a  correction  factor. 

(c)  Thermal  Penetration,  i.  e.  the  material  is  removed 
from  the  target  by  melting  or  vaporization. 

(d)  Explosive  Analogy,  i.e.  the  crater  is  considered 
to  be  identical  to  that  formed  by  an  equivalent 
mass  of  high  explosive  detonated  in  contact  with 
the  surface  of  the  target. 

3.1  Rigid  Projectile 

Several  theories  in  this  category  refer  to  penetration 
of  thin  targets  by  armor-piercing  projectiles,  and  as  such  are 
strictly  not  applicable  here.  They  are,  however,  included  for 
completeness.  These  include  Beth^'s  theory  of  penetration  of 
relatively  thick  target  plates^^,  later  extended  to  apply  to 


The  numbers  appearing  in  superscript  pertain  to  the  references 
at  the  end  of  the  report. 
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relatively  thin  targets  where  the  target  thickens  appreciably 
near  the  hole,  „  Thomson 11  s  theory  of  penetration  of 

very  thin  targets  which  bend  out  of  their  plane  during  pcne- 
tr&tion,  and  Zaid  and  Paul  s  theory  for 

very  thin  plates  which  petal.. 

Three  other  theories  in  this  category  were  proposed 

for  the  penetration  of  space  vehicles  by  meteoroids,  each 
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assuming  that  the  projectile  remains  rigid..  Griiraninger 

assumed  that  the  projectile  was  slowed  by  drag  during  a  ‘'fluid*1 

phase  of  the  penetration,  and  estimated  the  remaining  pene- 

tration  at  low  velocity  with  the  aid  of  an  empirical  armor- 

penetration  law.  Bohn  and  Fuchs  assumed  a  similar  "fluid” 

drag  force,  with  an  additional  static  force  proportional  to 

1 19  1 20 

the  Brinell  Hardness  of  the  target..  Zaid  *  considered 
a  rigid  projectile  which  is  separated  from  the  rigid  target 
by  a  zone  of  incompressible  fluid  composed  of  both  projectile 
and  target  material,  which  flows  radially  outwards,  eroding  the 
target  to  form  a  crater. 

3.1.1  Thick  Plate  Targets,  Bethe"’ 

Bethe  approached  the  problem  of  target  perforation 
by  finding  the  work  required  to  expand  a  hole  to  the  radius 
of  the  projectile  in  an  infinite  plate  of  thickness  T  ..  For 
the  thick  plate,  plane  strain  is  assumed  A  first  approxi¬ 
mation  to  the  dynamic  problem  is  obtained  by  using  the  static 
solution,  making  a  subsequent  correction  for*  dynamic  behavior.. 
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The  two-dimensional  static  stress  distribution  in 
an  infinite  plate  is 
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r  *r, 


is  the  elastic-plastic  boundary,  for  hole  radius  b  ,  which 
satisfies  the  Tresca  yield  criterion 

trs  -  <V  -  >'  (3-2) 

°r  3  are  t^e  an^  circumferential  stress 

respectively  at  radius  r  ,  V  the  Poisson's  ratio  for  the 
target  material  and  Y  the  yield  strength  nf  the  target 
material.  In  the  elastic  zone,  r  ?  r-j  ,  the  radial  dis¬ 
placement  vi  is  given  by 
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The  dynamic  equation  In  the  elastic  region 

2  ,  c rr  -  o~e 
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(3.4) 


is  now  solved  using  the  static  solution  as  a  first  approxima¬ 
tion.,  the  radius  of  the  hole  during  perforation  being  de¬ 
fined  by  the  ogival  shape  of  the  projectile. 
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where  R  is  the  projectile  radius,  L  the  length  of  ogive, 
and  V  the  projectile  velocity  assumed  constant,  during  per¬ 
foration,  A  second  approximation  to  the  solution  of  this 
equation  is  now  made  using  the  dynamic  boundary  conditions, 
in  place  of  the  static  boundary  conditions  initially  con¬ 


sidered  This  effect  will  be  of  the  relative  order 
where 
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In  order  to  complete  the  solution  of  the  problem,  it  is 
necessary  to  solve  the  plastic  equation, 
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where  the  yield  stress  Y  is  considered  constant.,  in  which 
the  earlier  static  and  dynamic  elastic  solutions  are  drawn 
upon  to  specify  the  boundary  conditions .  In  this  manner, 
neglecting  second  order  terms,  the  i-adial  stress  at  the 
hole  during  the  ogival  penetration  is  obtained .  The  total 
resistance  to  penetration  calculated  fox-  this  solution  is 

rR 

W  =  -  T  J  zir  b.o-r(b)  .  Jb 
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Bethe  shows  that  for  Orowan*3  dynamic  yield  condition. 


Yh  *  4  (  l  -  P  H  t) 

where  A  is  the  rate  of  straining,  the  total  resistance 
should  be  calculated  with  a  yield  corresponding  to  a  rate 
of  deformation 

7T  r  ~  2  W  Yo  V 
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(3.8) 
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3,1.2  Thin  Place  Targets,  Befche 
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During  the  perforation  of  a  thin  plate,  the  target 
is  found  to  thicken  appreciably  near  the  holeo  However, 
since  the  plate  is  thin,  plane  stress  is  still  assumed.  Be¬ 
cause  of  this  assumption  two  yield  criteria  are  necessary  to 
define  the  plastic  zone. 


oe  -  or  ”  Y  if  rr0  is  tensile. 


-  ar  “  Y  if  o0  is  compressive 


a  always  being  compressive.  The  analysis  of  the  resulting 
equations  is  far  more  complicated  than  for  the  thick-plate 
target,  Bethe  assumed  that  the  dynamic  correction  to  the 
initial  static  solution  would  be  the  same  as  for  a  thick 
plate  target. 
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G.  I.  Taylor^*"^  has  pointed  out  that  the  original 
thin-plate  solution  formulated  by  Bethe  did  not  use  the 
correct  stress-strain  relations  In  the  plastic  region,,  He 
resolved  the  equations  numerically  using  the  corrected 
relations  and  found  that  the  work  necessary  to  expand  the 
hole  statically  in  the  incompressible  target  material 
was 


W  =  13  3  ir  ft*  T  y 

The  total  resistance  to  penetration  is  the  sum  of  this 
and  Bethe' s  dynamic  correction  for  an  ogival  projectile, 
that  is 


W  =  IT  R'T  {  I.  33  y  t  0-9  72  8  ft 
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Here  the  yield  stress  is  taken  to  be  approximately  the  same 
as  that  for  a  thick  plate,  Le.  corresponding  to  a  rate  of 
deformation 
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3  o  1 . 3 


Very  Thin  Targets, 


Thomson 


107 


Statically  it  has  been  found  that  a  hole  can  be 
enlarged  symmetrically  without  bending  the  plate  until  its 
diameter  is  between  7  and  10  times  the  thickness  of  the 
plate.  At  about  this  stage  the  plate  always  bends  out  of 
its  plane,  and  the  material  of  the  tube  shaped  crater  will 
then  expand  under  the  influence  of  a  tangetial  tension  equal 
to  the  yield  stress  Y  ,  the  stress  perpendicular  to  the 
plane  of  the  plate  being  zero,  Thomson  and  others  consider 
this  mode  of  armor  penetration,  assuming  that  the  displace¬ 
ment  configuration  is  the  same  as  that  under  static  condi¬ 
tions,  but  taking  into  account  the  heating  effect,.  However, 
Thomson  concluded  that  the  energy  dissipation  by  heating  is 
very  small  compared  with  the  energy  dissipated  by  plastic 


deformation.  Thus  has  been  substantiated  by  some  experimental 

49 

work  carried  out  by  Krafft  at  the  Naval  Research  Laboratories. 

For  these  conditions  the  work  done  by  an  ogival 
projectile  in  perforation  is  found  to  be 
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3.1.4  '  Very  Thin  Targets,  Zaid  and  Paul 
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Zaid  and  Faul ! s  theory  considers  the  target  to 
be  of  incompressible  material  of  zero  strength.  In  the 
theory  it  is  assumed  that  the  non- intersected  part  of 
the  target  remains  essentially  undeformed,  as  it  is  con¬ 
sidered  that  the  plastic  waves  will  not  have  had  sufficient 
time  to  radiate  any  appreciable  distance..  The  intersected 
portion  of  the  target  is  assumed  to  take  on  the  shape  of  the 
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projectile's  surface „  Furthermore,  since  there  can  be  no 
tensile  force  between  the  projectile  and  the  petals  formed, 
contact  is  maintained  only  as  long  as  the  petal  momentum  is 
increasing  monotonically,  that  is,  the  petals  separate  from 
the  projectile,  when 


(petal  momentum)  =  0 
(petal  momentum)  0 


In  normal  impact  the  motion  is  axially  symmetric 
and  the  velocity  of  a  particle  In  the  plaLe  will  be 
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where  z.  is  the  coordinate  of  the  plate  out  of  its  plane, 
and  Z  the  distance  travelled  by  the  projectile  at.  constant 
velocity  from  its  initial  position  of  perforation.  The  axial 
momentum  transferred  to  the  plate  is  then 
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where  b  is  the  radius  of  the  deformed  portion  of  the  target. 
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For  a  truncated  projectile  a  further  term  KmqV  is  added 
to  the  above  expression  to  allow  for  the  detached  mass  of  the 


plate  which  I*.S.3  til  <3  Sutu8  radius  r  as  the  projectile 

m 

raeplat  (flat  nose  of  the  truncated  projectile) »  K  is 
an  empirical  constant  with  values  between  1  and  2..  The 
velocity  lost  by  the  projectile  will  then  be  given  by 


b 

j  2  ir  ft  T  V  jj  r  dr  (3„13) 


for  AV  <  <  V 

The  oblique  impact  of  a  truncated  conical  pro- 
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jectile  was  considered  in  a  similar  manner  Perforation 

was  considered  in  four  stages: 

(i)  start  of  the  slug  formation,  petal  element 
begins  to  form, 

(ii)  slug  has  been  completely  formed;  petals 
conform  to  the  conical  profile 

(iii)  petal  formation  continues,  petal  separation 
and 

(iv)  petal  formation  has  ended,  petal  separation 
well  underway. 

The  accompanying  diagram  illustrates  those  phases 
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The  decrease  in  projectile  velocity  may  be 
computed  from  the  relation 


Momentum  transferred  to  target 


m 


where  A  V  << 


H7 

3.1o5  Gritnminger 1  g  Theory 

Grimminger  considered  that  penetration  occurred 
in  two  phases..  During  the  first,  the  rigid  projectile 
moves  through  the  target  material,  which  behaves  like  a 
compressible  perfect  fluid,  exerting  a  drag  force  on  the 
proiectile..  For  very  high  initial  projectile  Mach  numbers 
(ratio  of  velocity  to  the  velocity  of  plastic  deformation 
waves)  the  drag  coefficient  was  assumed  to  be  unity.  For 
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a  spherical  projectile  the  resultant  equation  of  motion  is 
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Assuming  that  this  process  continues  until  the 
projectile  is  slowed  to  a  Mach  number  of  5,  the  penetration 
is 
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where  c  is  the  velocity  of  plastic  deformation  waves. 

Subsequent  penetration  below  a  Mach  number  of  5 
was  assumed  to  be  given  by  the  empirical  armor  penetration 
formula 
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where  K  is  empirically  defined  as 


(3.16) 


Ultimate  tensile  strength  of  Target  Material 

K  „  5  o03  - - - — - - - - - - - - — - 

Ultimate  tensile  strength  of  Copper 

with  units  ft.  lb ./ft.  . 

The  total  penetration  therefore  is 
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where  the  speed  c  of  plastic  waves  was  taken  as  lUOO 
Grimminger  considered  that  the  assumption  tha 


(3.17) 
ft  ./sec.. 
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target  melting  could  be  neglected  was  reasonable,,  since  he 
considered  the  transfer  of  momentum  to  be  a  more  rapid  pro¬ 
cess  than  the  transfer  of  heat.,  and,  therefore.,  the  penetra¬ 
tion  xi/ould  be  completed  before  the  full  effects  of  heating 
would  be  felto 

3,1,6  Bohn  and  Fuchs 8  Theory^ 

Bohn  and  Fuchs  assumed  the  resistance  to  projectile 
penetration  was  made  up  of  two  parts,  a  static  part 

TT  R'H 

and  a  dynamic  part 

jtJjl  r'v" 

z 

where  H  is  the  Brinell  Hardness  number  for  the  target,  R 
the  radius  of  the  projectile  and  f  a  shape  factor  asso¬ 
ciated  with  the  projectile.  This  shape  factor  was  taken  to 
be  a  function  of  the  angle  at  which  a  section  of  the  surface 
of  the  projectile  is  inclined  towards  the  direction  of  flight, 
defined  by 

j  Sm  Cp  f) 

f  =  - (3a8) 

where  dA  is  the  element  of  the  projectile  area  perpendicular 
to  the  direction  of  flight,  and  0  the  inclination  of  that 
element  to  that  direction.  The  following  values  are  given, 

for  micrometeorites  £  =  1 

spheres  f  '=  2/3 

pointed  projectiles  f  «=•  1/3 


The  total  resistance  encountered  by  the  projectile 
was  then  assumed  to  be 


(3.19) 


This  may  be  integrated  to  find  the  expression 
for  the  penetration 
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For  a  spherical  projectile 
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3..  1,7  Zard  s  Theory 


The  most  recent,  theory  to  estimate  penetration 
at  very  high  velocity  has  been  postulated  by  Zaid.  He 
assumed  that  the  shock  waves  in  the  projectile  and  target 
could  not  move  at  high  enough  velocity  to  move  away  from 
the  projectile- target  interface,  and  therefore  the  pro¬ 
jectile  and  target  could  be  regarded  as  rigid,  being 


separated  by  a  tain  layer  of  material  encompassed  by  the  shock 

zone  whi ch 
was  assumed  to 
be  fluid o  The 
fluid  is  assumed 
to  move  radially 
as  the  projectile 
advances  eroding 
the  target  to  form 
the  crater,  Iner¬ 
tial  expansion  of 
the  crater  is  neg¬ 
lected  .. 

The  further 
assumptions  neces¬ 
sary  to  obtain  an 
equation  of  motion  were  (a)  the  thickness  of  the  fluid  zone 
is  a  constant  independent,  of  radial  position,  (b)  at  radial 
distances  greater  than  the  projectile  radius,  the  dispelled 
particles  do  not  interact  with  the  material  in  the  fluid 
zone,  so  that  the  fluid  pressure  at,  the  outer  radius  of  the 
cylindrical  fluid  zone  is  essentially  zero,  (c)  the  fluid-zone 
thickness  is  small  compared  with  the  penetrated  distance,  and 
to  the  radius  of  the  projectile,  which  is  assumed  to  be  cylin¬ 
drical. 

From  these  assumptions,  the  following  equation 
of  motion  is  obtained 


R  R 


0  r 
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where  p  is  the  penetration  of  the  projectile,  (*)  repre¬ 
sents  differentiation  with  respect  to  time,  h  is  the  fluid 
zona  thickness,  tu  the  shear  at  the  proj ectile- fluid  zone 
interface,  and  K  an  averaging  factor  which  is  the  product 
of  two  averaging  factors,  i»e»  the  averaging  factor  for  the 
pressure  gradient  djp/ur  over  the  height  h  of  the  fluid 
zone,  and  the  ratio  of  pressure  in  t.he  fluid  to  that  at  the 
proj  ectile- fluid  zone  interface..  Solution  of  the  equation 
depends  on  assuming  relations  for  the  interface  shear  and  the 
zone  growth.  The  penetration  depends  strongly  on  the  assump¬ 
tions  made,  and  satisfactory  solutions  have  so  far  not  been 
given o  Further  details  are  given  in  Appendix  A. 


3 . 2  Hydrodynamic 
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Birkhoff  et  al  “  suggested  that  since  pressures 
generated  in  high-velocity  impact  are  so  much  greater  than 
the  shear  strengths  of  the  projectile  and  target,  the  shear 
strengths  might  be  neglected  find  the  projectile  and  target 
treated  as  inviscid  fluids.  Under  this  assumption,  the 
resistance  to  penetration  arises  entirely  from  the  inertial 
forces  required  to  accelerate  the  target  material.  Theories 
based  on  this  assumption,  therefore,  have  been  termed  hydro- 
dynamic  theories,. 


The  most  detailed  solution  is  that  of  Bjork^^ 
who  integrated  the  two-dimensional  equations  of  motion  of 
a  compressible  inviscid  fluid,  together  with  an  entropic 
equation  of  state  of  the  solid  material,  using  numerical 
methods..  Analytical  solutions  of  t:h.e  nonlinear  equations 
could  not  be  obtained 


v  H 


Solutions  may  be  obtained  more  simply  if  the 

additional  assumption  is  made  that  the  material  is  incom- 
121 

pressible.  Opik  used  a  simplified  theory  of  this  type, 
including  a  correction  term  for  target  strength.. 

Further  simplifications  arise  if  only  the  longi¬ 
tudinal  portion  or  the  motion  is  considered,  and  the  pro¬ 
jectile  is  regarded  as  a  segment  of  a  fluid  jet  impinging 
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on  a  fluid  target,  e.g.  Birkhoff  et  al  and  Pack  et  al  '  . 

For  these  theories,  motion  continues  only  as  long  as  the 
jet:  continues  to  impinge  on  the  target,  and  inertial  expan¬ 
sion  of  the  target  after  dissipation  of  the  jet  is  neglected. 
One  would  expect,  therefore,  that  these  theories  would  give 
only  approximate  agreement  for  very  high  velocity  projectiles 
of  relatively  great  length  compared  with  diameter,  e.g. 
shaped  charge  jets  or  slender  rods.  Theories  of  this  type 
have  been  used  nevertheless  to  estimate  penetration  by 
short  projectiles. 


3.2.1  Biork's  Theory128 


In  this  approach  to  cratering  of  a  target  by  a 
high-velocity  projectile,  the  behavior  of  the  material  is 
described  by  an  entropic  equation  of  state  that  relates  the 
internal  energy,  pressure,  and  specific  volume,  buL  neglects 
the  shear  strength  of  the  material.  The  problem  solved  by 
Bjork  was  for  a  cylindrical  projectile  impacting  normally 
on  a  semi-infinite  target.  The  equations  of  motion  are 
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P  tT"  +  P  V  8  +  P  cl  tv  V  -  Q 

J  J  J  (3.22) 

P'  Hf.e) 

where  v  is  the  fluid  velocity  vector.,  P  the  pressure, 
e  the  internal  energy,  p  the  density  of  both  the  target 
and  projectile,  and  t  the  time  measured  from  the  time  of 
initial  contact o  The  equations  specify  an  initial  value 
problem  for  which  the  initial  conditions  are  t  »  0,  P  »  0, 
e  *  0  everywhere,  the  initial  velocity  in  the  target  every¬ 
where  zero,  and  the  initial  velocity  V  of  the  projectile 
normal  to  the  target  surface.. 

Equations  of  state  of  solids  under  hydrodynamic 
conditions  are  available  from  work  at  Los  Alamos  and  else¬ 
where'^”1'^  for  pressures  up  to  several  megabars,  and 
from  the  Ferrai-Thomas-Dirac  theory  for  much  higher  pres¬ 
sures  .  A  numerical  solution  was  obtained  by  means  of  a 
high  speed  digital  computer.  The  machine  representation 
is  fairly  rough,  since  a  total  space  mesh  of  25  by  25  cells 
was  used.  However,  the  main  features  of  penetration  are 
displayed  quite  clearly,, 

Bjork's  theoretical  results  are  fitted  approxi¬ 
mately  by 

-f  =  k  V  /JJ  (3  .23) 

where  &  =  (  ^  R  L )  (3,24) 

is  considered  as  the  representative  dimension  of  the  projectile, 
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3.2.2  Opik* s  Theory121, 124 

The  first  theory  put  forward  on  cratering  by 
meteoroids  was  given  by  Opik  in  1936..  In  this  theory.,  an 
allowance  is  made  for  the  yield  strength  of  the  target 
material. 


A  cylindrical  projectile  (L/d  =  1)  of  radius  R  ,  density  Pp  , 
mass  mp  ,  and  velocity  V  impacts  a  semi-infinite  target  of 
density  Pt  »  The  projectile  is  considered  to  be  incompres¬ 
sible,  and  its  shape  to  be  defined  after  time  t  by  p  ,  the 
penetration  below  the  original  surface,  and  r  the  radius  of 
the  crater  occupied  by  the  projectile,  where 


1.  rr  R  3  fp  =  2  ir  r  x  (  f>  -  f  )  f  h  (3  .  25  ) 

Here  p  is  the  depth  of  the  center  of  gravity  of  the  deformed 
project4  In  vc1.''W  the  original  surface.  The  momentum  equation 
then  gives 

yv\  p  -  ■■■  ft  (  I 
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(3  .,26) 


where  P  is  the  pressure  at  the  base  of  the  crater  which 
resists  penetration*  Applying  Bernoulli9s  thereto 

P  =■  i  t-  k  (3-27) 

where  k  is  the  strength  factor  of  the  target,  i*e*  it  is 
the  minimum  pressure  for  penetration  to  occur*  Similarly 
for  the  radial  expansion, 

P'  =  £  rX  +  K  (3*28) 

Opik  relates  P  and  P*  by  applying  Bernoulli's 
equation  in  the  frame  of  reference  moving  with  velocity  p  , 

P  -  P'“  +  r  /V  r  X  <3  29) 

Since  the  flow  in  this  frame  of  reference  is  not  steady, 
this  procedure  can  at  best  be  only  a  crude  approximation.. 

The  differential  equation  resulting  from  the 
elimination  of  the  pressure  terms  from  equations  3.25 
to  3*29  is  solved  subject  to  the  initial  conditions, 
t  =  0,  r  *=>  R,  p  »  0,  p  =  v  Details  of  the  integration 
are  given  in  Appendix  B« 

3.2.3 _ Shaped  Charge  Jet  Penetration 

Several  investigators  (Birkhoff  et;  alx^  , 
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Cook  ,  Pack  et  al  "  )  have  proposed  similar  incompres 
sible  theories  to  explain  the  penetration  produced  by  a 
shaped-charge  jet,.  However,  cylindrical  or  shpericai  pro- 
"  -  f  ;1  ■"  '■mo  'inmetimf'n  ilnn  considered  as  short  iets 

In  this  theory  both  the  jet  and  target  material  are  regarded 


2  7 


as  incompressible  fluids.  The  flow,  when  viewed  from  a 
frame  of  reference  moving  with  the  velocity  of  the  jet-target 
interface  v  ,  is  considered  to  be  steady..  Applying  Bernoulli's 


equation  along  the  axis  of  symmetry,  the  stagnation  pressure 
is  given  by 

P  =  i  ft  -  *  fp  (y~v)  (3.30) 

Originally,  the  strength  of  the  material  was 
neglected  in  comparison  with  the  high  pressures  produced 

IQ/ 

by  the  impacting  jet.  Rostoker  included  the  target 
strength  in  this  equation  by  inciuding  a  factor  k  , 

£  ft  k  =  i  fP  W  ~V)X  (3.31) 

1 2  ^ 

where  k  is  the  target's  dynamic  strength.  Ei  chp'lbercor 
writes  this  in  the  foim 

k  =  if,  Kl  (3.32) 
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where  v  is  the  minimum  velocity  of  the  jet  that  will 
0 

just  cause  a  crater  to  form  in  the  target.. 

For  a  jet  moving  at  constant  velocity  V  and 
of  density  P  ,  and  finite  length  L  ,  most  of  the  pene¬ 
tration  is  considered  to  take  place  during  the  steady  part 
of  the  motion,  so  that  penetration  p  is  approximately 

r  .  /  V.At  (3-33) 


where  is  the  total  time  elapsed  from  the  time  that 

the  jet  first  struck  the  target. 


Hence 


z 

L 

Foe  a  jet  of  high  velocity,  such  that  V  > >  vq 


p  —  fjS  (3,35) 

L  ~~  /  Pt 


which  would  have  been  obtained  directly  from 
which  disregards  the  strength  or  the  tax. get. 


eqnn  t  i  on  3 , 30 , 


Pack  et  al  made  a  corroctf  on  for  the  effort 
of  strength  by  considering  the  penetration  as  a  series  of 
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powers  of  the  nondimenaional  parameter  Y/p^V” 
approximation  being 

jl  _  r~ft  (  I  -  — ) 

L  -/  ft  l  '  r"'"} 
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(3=36) 


where  Y  is  the  dynamic  yield  stress  of  the  target ,  and  ct^ 
is  an  empirical  function  of  the  densities  of  the  jet  and 
of  the  tergeto  For 'armor,  the  correction  term  tt-jY/p  V4 


was 


found  to  be  as  high 

As  has  be 
not  consider  the  in 
jet  ceases  to  imp  in] 


as  0„3„ 


=n  pointed  out  earlier,  this  theory  does 
srtial  expansion  of  the  crater  after  the 
ige  on  the  targe t„  In  their  experiments 
with  shaped  charge  jets.  Pack  et  al  found  that  the  shape 
at  the  bottom  of  the  crater  was  roughly  hemispherical  with 
radius  equal  to  that  of  the  hole  made  by  the  jet  They 
reasoned  that  if  the  previous  pressure  which  suddenly  ceased 
to  apply  after  the  jet  had  been  used  up,  was  uniformly  dis¬ 
tributed  over  the  bottom  of  the  hole,  the  latter  expanding 
freely  adopted  a  hemispherical  shape .  pack  et  al  therefore 
modified  the  expression  for  the  penetration  to 
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where  r  is  the  radius  of  the  hole  made  by  the  jet 
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Cook  has  extended  the  theory  to  estimate  the 
radius  of  the  hole  made  by  the  jet..  He  assumed  that  the 
impulse  imparted  to  the  Largcc  is  us>ed  up  in  cue  iaeerai 
expansion  of  the  hole,  i  e., 

*  ft,  (V-  ■-/  R'  at  =  (r  ft  /-)  vrr'At 
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(3. .38) 


where  R  is  the  radius  o£  the  jet,  and  r  that  of  the 
crater.  The  final  cross-sectional  area  will  bo  given  when 
r  ■=  0  ,  i.e, 

i/— 

/ 

Eliminating  v  by  means  of  Eq,  (3,30),  the  radius  of  the 
hole  is 


(3,39) 


3,3  Thermal  Penetration 

In  the  thermal  penetration  theories  that  have 
been  formulated,  penetration  is  assumed  to  take  place  by 
melting  or  vaporization  of  the  target  by  the  projectile 
as  its  energy  is  dissipated, 

3,3.1  Whipple^  Theory'*' ^ ^ 

Whipple  assumes  that  the  crater  formed  by  melting 
is  a  right- circular  cone  for  which  the  total  apex  angle  is 
60° „  Assuming  that  all  of  the  initial  kinetic  energy  is 
converted  into  heat,  the  depth  of  penetration  is  simply 
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where  Qt  is  the  heat  required  to  remove  a  unit  mass  of 
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the  target  material,  here  taken  to  be  the  latent  heat  to 
fusion-  For  an  equivalent  spherical  projectile  m^  *  ~ 
so  that 


JZ  (  ^ 

*  (  'f-  U  Q-f  '  (3  „  40a) 

3-3-2  Langton8s  Theory^^ 

Langton's  theory  is  similar  to  that  given  by  Whipple 
since  in  this  case  also,  the  energy  of  the  projectile  is 
assumed  to  be  converted  into  heat-  Four  possible  cases  are 
considered, 

(a)  the  projectile  velocity  is  such  that  neither  the 
projectile  nor  target  melt 

(b)  The  projectile  is  melted  but  the  target  is  not 

(c)  The  target  is  melted  but  the  projectile  is  not: 

(d)  Both  projectile  and  target  melt 

Penetration  would  only  occur  for  cases  (c)  and  (d).. 

A  hemispherical  crater  is  assumed  to  be  formed  Case  (d) 
results  in  the  equation 

_  i 
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w!mrr  O  0  'X-.  w  t  h>,  uwut  uei.  uiui.  uiUj;j  u  1  l.ai.  aud 

'C  p  K 

projectile  material  respectively,  necessary  to  raise  the 
temperature  oi  the  material  to  the  melting  point,  and  to 
supply  the  latent  heat  of  fusion  Thus,  the  penetration  is 


given  by 
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For  an  equivalent  spherical  projectile  this  becomes 
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3 . 3 ..  3  Grow  8  s  Theory ' ' 

A  similar  theory  has  been  developed  by  Grow  for 
the  energy  lost  by  a  projectile  In  perforating  a  thin 
target.  The  projectile  was  assumed  not  to  melt  or  deform, 
but  a  hole  equal  in  diameter  to  the  projectile  was  assumed 
to  be  formed  in  the  target  by  melting..  Denoting  the  exit 
velocity  by  V  ,  the  resulting  expression  Is 
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where  T  is  the  target  thickness  and  R  the  projectile 

radius.  Here  Q  is  taken  to  be  the  energy  per  unit  mass 

required  to  raise  the  temperature  of  the  target  to  the 

melting  point  and  supply  the  latent  heat  of  fusion,  but. 

Qp  is  taken  tn  bo  the  energy  per  unit  mass  necessary  to 

raise  the  project  i  n±  uiuL.ei.iul  cu  Lui  gef.  tempera trrM , 

projectile  not  being  assumed  to  melt.  Grow  further  assumed 

that  the  last  term  could  be  approximate  by  kV  ,  varying 

from  zero  for  the  velocity  at  which  negligible  deformation 

r>/»nnvr«»H .  to  n  maximum  0  for  the  velocity  at  which  the 
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projectile  shatters., 


The  theory  may  be  readily  extended  to  cratering 
in  thick  targets  where  no  perforation  occurs  and  the  exit 
velocity  is  zero,,  The  penetration  will  then  equal  T  ,  i.e, 
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For  an  equivalent  spherical  projectile.,  this  becomes 
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3  <.3, 4  Lavrent 8 yev*s  Theory 

Lavreut ' yev  considers  two  cases  of  thermal 
penetration,  a  rod  target  and  a  hemispherical  plate  target, 
both  being  of  incompressible  material ,  The  analysis  is 
idealized  by  considering  the  medii  to  be  made  up  of  a  set 
of  infinitesimally- thin  layers  which  impact  inelastically, 
and  for  which  the  energy  dissipation  distribution  is  evaluated 
If  the  density  of  this  distribution  is  greater  than  the  amount 
required  to  vaporize  the  material,  a  crater  will  be  formed. 


The  first  case  considers  a  plate,  depth  L  ,  of 
initial  velocity  V  ,  impacting  upon  a  rod  of  the  same 
material  and  cross-section  If  a  total  length  z  of  material 
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is  set  in  motion  at  time  t  ,  the  conservation  of 
momentum  gives 

e  <?l i/  +-  v  d  z  ==•  0  (3o45) 

which  on  integration  for  the  initial  condition  z  =  L  , 
v  **  V,  yields 
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The  kinetic  energy  is  given  by 


or  differentiating 

ci  £  *  r  +  z  \S  dif 

which  by  Eq.  (3  <,45)  becomes 


Assuming  that  all  the  kinetic  energy  is  converted  into 
heat,  its  distribution  along  the  rod  will  be 
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If  the  heat  density  necessary  to  vapourise  the  material 
■f  s  0  ,  then  a  length 


(3.46) 
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will  be  vapourised,  this  length  including  the  length 
of  the  projectile.  A  similar  analysis  was  given  assuming 
spherical  symmetry  for  a  hemispherical  plate  projectile = 


3.4  Explosive  Analogy 
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Several  investigators,  Baldwin,  Rinehart  and 
Stanyukovitcii'"'’-'  have  proposed  that  crater  formation  due 
to  high  velocity  impact  is  similar  to  crater  formation  by 
high  explosives.  Allison.^^  has  given  a  review  of  this 
subject,  illustrating  the  similarity  of  craters  formed  in 
lead  targets  by  a  projectile  with  a  kinetic  energy  of  93  x  10' 
Joules,  and  by  a  tetryl  charge  with  an  explosive  energy  of 
63  x  10^  Joules. 


3.4.1  Stanyukovltch* s  Theory^^-* 

Stanyukovltch  considers  the  phenomenon  of 
cratering  to  be  similar  for  both  the  high  velocity  impact 
and  an  explosion.  The  pressure  on  the  front  of  the  shock 
wave  caused  bv 'either  the  explosion  or  impact,  drops  sharply 
with  increasing  distance  from  the  point  of  initiation,  and 
the  mechanism  of  cratering  changes  from  vaporization,  to 
melting  and  pulverization  of  the  material  This  ceases 
when  the  energy  density  at  the  shock  front  becomes  less 
than  a  certain  value.  However,  the  evaporated  and  finely 
fragmented  material  expands,  and  has  an  effect  similar  to 
an  explosion, 

The  equivalent  mass  of  explosive  is  determined 
by  the  relation 
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(3,47) 
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where  mc  is  Line  mass  of  the  equivalent  explosive  material, 
i]  the  explosive  efficiency,  and  Q  the  energy/unit  mass 
of  explosive. 

The  approach  used  by  Stanyukovitch  for  estimating 
the  radius  of  the  crater  is  empirical .  See  Appendix  C 
for  a  more  detailed  account  of  his  analysis. 
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SECTION  IV 


EMPIRICAL  RELATIONS 


Due  to  the  lack  of  adequate  theories ,  many 
attempts  have  been  made  to  correlate  experimental  cratering 
data  with  the  aid  of  empirical  expressions .  Many  of  these 
correlations  have  been  based  on  simple  functions  of  inde¬ 
pendent  quantities,  such  as 

•  i  \y  J 

=  X  V 

where  ,  and  V  are  the  independent  variables  projec¬ 
tile  mass  and  velocity,  respectively. 

Other  correlations  have  been  attempted  using 
noirdimensional  parameters,  so  that  the  constants  to  be 
fitted  empirically  are  non-dimensional 

More  refined  correlation  schemes  have  been 
evolved  by  assuming  some  type  of  resistance  law  based  on 
the  physical  phenomena  considered  to  be  dominant.  The 
resultant  equations  contain  non-dimensional  constants  which 
are  then  fitted  to  the  data. 

In  order  to  gain  perspective,  the  low  velocity 
armor  penetration  work  is  included  first  in  the  following 
discussion.  After  presenting  the  high  velocity  correlation 
equations,  a  general  discu;  '  'n  of  dimensional  analysis  of 
the  cratering  problem  is  given,  and  the  important  parameters 
based  on  elementary  static  material  properties  are  discussed, 


4,1  Early  Armor  Penetration  Work 


The  study  of  ballistic  penetration  is  generally 

considered  to  have  begun  with  the  work  of  the  English 

ballistician  Robins 1742)  whose  ideas  were  elaborated 
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a  few  years  later  by  Euler ~ "  ,  both  in  theory  and  experi¬ 
ment. 

In  armor  penetration,  the  object  has  been  to 
estimate  the  limiting  velocity,  i.e.  the  minimum  velocity  of 
a  projectile  necessary  to  perforate  the  target  completely. 
This  was  done  by  assuming  a  law  of  penetration  based  on  the 
physical  phenomena  assumed  to  predominate,  solving  the 
ensuing  equation,  and  fitting  the  resulting  expression  to 
experimental  data  with  suitable  constants.  Most  of  the 
laws  postulated  are  variations  of  the  form 

V  r.  k  T'1'  Rc  (4,1) 

l-  f 


where 

r 


and 


iv 

k. 


Is  the  limiting  perforation  velocity, 
the  thickne  ss  of  the  target 
the  weight  of  the  projectile 
the  radius  of  the  projectile 
a  constant. 
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In  the  theory  postulated  by  Robins'  the  reals 
tance  to  penetr~at.ion  was  taken  to  be  constant,  hence 


leading  to 


2  hv\ 
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(4.2) 
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Euler  evaluated  k  from  experiments  on  elm  wood 

and  earth. 

Workers  at  the  Naval  Research  Laboratory'1'0"1 
modified  this  equation  in  order  to  allow  for  free  surface 
effects.  The  idea  involved  is  that  the  material  near  the 
free  surface  is  more  easily  pushed  aside  than  material  near 
the  plate  center  line.  This  is  taken  into  account  by  sub¬ 
tracting  from  the  plate  thickness  an  amount  proportional 
to  R  ,  Thus, 

AT.*  K  (4 „ 3 ) 

Poncelot  (1829) '10^J  published  a  theory  which 
combined  the  Euler-Robins  type  resisting  force  with  a  force 
proportional  to  the  square  of  the  velocity.  This  was  pro¬ 
posed  for  application  to  materials  weak  enough  so  that  their 
inertial  resistance  to  motion  is  significant.  The  constant 
force  k  is  interpreted  as  the  pressure  necessary  to  break 
the  cohesion  of  the  target.  Thus,  the  resistance  is  ex¬ 
pressed  by 

W  =  7 r  R  ("  k  +-  T  j°  f  V  j 


which  leads  to  the  differential  equation 


v  ,j,= 


ctv 

r- 


=  -  i r  (  k  + 
\ 


where  p  is  the  distance  penetrated  into  the  target,  The 
integration  of  this  equation  yields 
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a  is  empirically  found  to  be  not  far  from  1/2 . 
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In  1945  Roberson  and  Irwin  worked  out  a  similar 
expression  in  which  account  was  taken  of  the  target  cohesive 
resistance.,  the  inertial  resistance,  and  the  free-surfaee 
effects o  They  presented  an  equation  in  the  form 
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(4,5) 


They  also  took  into  account  the  increase  of  projectile  cross- 
sectional  area  due  to  projectile  deformation  in  the  high-velo¬ 
city  range.  The  value  of  the  area  was  obtained  by  measuring 
the  diameter  of  the  projectile  after  deformation,,  They 
found  that  the  radius  of  the  projectile  could  be  expressed 
empirically  as 


R  -  R( 
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for  impact  velocities  greater  than  2,500  ft/sec. 

Early  experimental  work^0Zf  showed  that  the  force 
required  to  perforate  a  plate  varied  approximately  as  the 
plate  thickness  if  the  diameter  is  constant,  and  as  the 

A  -i  -I  ■T  t*Vrf  1*'  *  f^1r 

*  ^  “  1  ‘  "  "  '*  '  KS 

these  assumptions,  the  force  is  proportional  to  TR,  con¬ 
sequently,  the  work  done  to  perforate  the  plate  is  proportional 
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to  T  R  .  Hence 


(4.7) 


Even  over  a  restricted  velocity  range  it  was  found  that 
there  is  a  variation  in  the  proportionality  constant „ 

For  cast-iron  shot  on  wrought  iron  plate,  Fairbairn^^lSbl) 
derived 


)/  -  16*0 


,  T  R 


(4.8) 


where  T  and  R  are  in  inches,  mp  in  pounds  and  in 

ft/sec „ 

On  the  assumption  that  in  high  velocity  penetra- 

1 

tion  the  force  required  is  propostional  to  TR/V,  Tresidder 
obtained  for  the  perforation  of  wrought  iron 


V/.  «  k 


r  *> 
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where  log  k  ®  3 ,,0473,  for  T  and  R  inches,  in 

pounds,  and  V.  in  ft/sec.  The  Krupps10^  expression  for 

JL» 

wrought  iron  is 


^  -k  TV>  *** 

VL  *  _  .A 


(4.10) 


where  log  k  ■»  3.1397..,  while  De  Marre'*"w  '  postulated 


a-cr-  o  ?/.- 

V  -  k  71 _ 1- 

L  ^  or 


(411) 
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where  log  k  *  3*1874,  for  the  same  material* 

A  later  expression  adopted  by  Krupps^0^  for  hard 
faced  plates  is 


T  R 


m. 


i/x 


(4.12) 


where  log  k  =  3*3271*  And  for  ordinary  steel.  De  Marre^0^ 
modified  his  wrought  iron  expression  to 


Vi  =  k  — 


•  0‘-7  0-7 s~ 
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(4*13) 


where  log  k  <=  3*23522*  In  order  to  correlate  the  resistance 
of  various  targets  of  the  same  thickness,  the  "De  Marre 
coefficient"  was  adopted*  This  is  the  ratio  between  the  vel¬ 
ocity  calculated  from  the  limiting  velocity  expression  for 
ordinary  steel. 
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About  1930,  Robert  Kent  '  put  into  use  a  slightly 
modified  version  of  the  De  Marre  expression*  This  relation 
is 


1/, 


k 


(4*14) 


and  has  been  used  extensively  for  armor  perforation  and  pane- 
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tration  work  At  about  the  same  time,  L..T.E  Thompson"  at 
Dahlgren  introduced  the  expression 
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where  ,0  is  the  angle  of  obliquity  of  the  projectile's 
path  to  the  plate  normal.  This  expression  has  been  used 
extensively  by  the  Naval  authorities  for  the  design  of 
ship  armor  plating. 

The  above  expressions  are  all  applicable  to 
plate  penetration  where  the  projectile  is  essentially  intact 
after  penetration,  but  are  of  no  use  at  higher  velocities 
where  the  projectile  disintegrates. 


4.2  High  Velocity  Penetration  Work 


In  order  to  correlate  the  experimental  data  on 
penetration  depths  and  volumes  of  craters  formed  in  high- 
velocity  impact,  numerous  empirical  expressions  have  been 
proposed.  Some  of  the  most  widely  used  expressions  are 
listed  in  Table  III,  together  with  conditions  for  which 
the  constants  were  evaluated. 


One  of  the  earliest  and  most  frequently  used 
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assumptions  was  put  forward  by  Rinehart  and  Pearson 
They  assumed  that  the  crater  volume  was  proportional  to 
the  projectile  kinetic  energy.  If,  in  addition,  it  is 
assumed  that  the  crater  shape  does  not  change  with  velocity, 
then  taking  the  crater  volume  proportional  to  the  projectile 
kinetic  energy  leads  to  an  expression  of  the  form 


d 


(4.16) 


where  p  is  the  depth  of  penetration,  d  the  diameter  of 
the  projectile,  the  projectile  mass,  and  V  the  pro¬ 
jectile  velocity. 
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14 

On  the  other  hand  Pugh  and  Eicheiberger  found 
that  the  crater  depth  was  proportional  to  the  projectile 
momentum,  leading  to  an  expression  of  the  form 


t  \ 


(4.17) 


Several  other  equations  have  been  put  forward The 
De  Marre  formula  for  plates  has  been  adapted  for  use  with 
semi- infinite  targets  at  low  velocities  of  impact .  If  the 
plate  thickness  is  regarded  to  be  the  penetration  depth 
corresponding  to  the  limiting  velocity,  the  follovzing  expres¬ 
sion  results 
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Partridge  et  al  proposed  that  the  crater  volume  for  lead 
to  load  impact  could  be  represented  by  the  equation 


Vr  '■  V,9  I  I  ~  exp  (~  k  £ 


(4  ,.19) 


where  V  is  the  crater  volume,  V  the  volume  of  the  pro- 
c  p 

jeetile,  and  E  the  projectile  kinetic  energy. 

All  the  equations  as  they  have  been  written  above 
require  a  different  empirical  constant  for  different  projec¬ 
tile-target  material  combinations.  The  next  step,  therefore, 

was  to  relate  the  constants  to  material  properties  in  order 
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to  obtain  a  more  general  correlation.,  Hutn  in  examining 
impact  data  for  projectiles  and  targets  of  similar  material 
proposed  a  possible  dimensionless  relationship  of  the  form 


*  (k/ar 
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where  was  the  sonic  velocity  in  the  target,  Qt 

the  energy  necessary  to  melt  or  vaporizse  a  unit  mass  of 
the  target  material,  and  j=c  the  thermal  diffusivity,  He 
found  that  his  data  correlated  best  on  the  basis  of  V/ct 
where  cfc  is  the  longitudinal  bar  velocity  of  sound. 

The  velocity  range  in  which  the  tests  were  carried  out 
was  0,1  <  (V/Cj.)  <  1,0  »  In  this  range  the  terms  V^/Qt 
and  Vd/K  were  found  to  be  insignificant.  He  also  con¬ 


cluded  that  the  crater  volume  was  not  precisely  propor¬ 
tional  to  either  the  energy  or  the  momentum  of  the  pro¬ 
jectile. 


On  the  other  hand  van  Valkenburg  et  al  '  found 
that  the  volume- energy  relation  could  be  correlated  for 
nonalloy  metals  using  p  and  c,.  ,  that  is 
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The  equation  reduces,  on  assuming 


I 


where  D  is  the  diameter  ot  the  mouth  of  crater,  which 
in  this  case  is  taken  to  be  a  section  of  a  sphere,  to 
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Van 

Valkenburg  et  al 

measured  L,j: 

crater  dimension*: 

to  the  top  of  the  projectile  material  remaining  in  the  crater. 
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The  parameter  (V/e)  has  subsequently  been  used 
quite  extensively,  various  experimenters  finding  different 
exponents  to  fit  their  own  data.  Charters  et  al  added 
a  term  consisting  of  the  ratio  of  projectile  to  target 
density  to  account  for  a  variety  of  projectile  materials 
finding  that  data  in  lead  and  copper  targets  were  best 
fitted  by 


Several  other  groups  used  this  expression  or  expressions 
based  on  these  parameters,  e.g.  at  Pou'iter  Laboratories 
and  GARDE. 
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Partridge  et  al  ~  and  Collins  et  al'  on  the 

other  hand  noted  that  their  data  seemed  to  be  best  fitted 

by  the  first  power  of  velocity.  The  Langley  investigators 
theorized  that  penetration  might  be  proportional  to  the 
momentum  per  unit  presented  area  of  the  projectile, 

r  =  k,  (  /„  V  i.  -  k,)  (4  23} 

where  kz  is  the  momentum  per  unit  area  below  which  no  per¬ 
manent  deformation  of  the  target  was  observed..  They  further 
found  that  the  proportionality  constant  was  a  function  of  the 
projectile  density  and  modulus  of  elasticity  of  the  target, 
finding  that  their  data  could  be  fitted  by  the  expression 

^  „  _ 30  9  6  ( fr  ^  L  *  ) _  (4  ,24) 
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Apart  from  Collins  and  Kinard ’ s  expression, 
none  of  those  mentioned  so  far  have  included  a.  term  which 
allows  for  either  the  target  or  projectile  strength.. 
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Pugh  and  Eichelberger A  first  introduced  the 
effect  of  target  strength  through  the  use  of  the  Brinell 
Hardness  number «  They  found  that  the  onstant  in  the  ex¬ 
pression 

.  % '  ,/  % 
p  =  l<  Wp  V 


could  be  reduced  for  various  target  materials  by  the 
Brinell  Hardness  number  and  the  density,,  and  found  that 
experimental  data  for  several  materials  could  be  corre¬ 
lated  by 
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Experimenters  at.  the  Ballistic  Research  Labor- 
atoryb  subsequently  noted  that  the  energy  per  unit  volume 
tended  to  become  constant  as  the  velocity  increased,  and 
found  that  this  constant  could  be  correlated  roughly  with 
Brinell  Hardness  number..  Assuming  a  constant  crater  shape, 
this  leads  to  an  expression  of  the  form 

V3  Vj 

p  =  k  (426) 

J  Nt 

A  similar  expression  is  also  in  use  at  the  Naval  Research 
Laboratory . 
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Engel  introduced  the  yield  strength  as  a 
correction  term  to  the  equation  which  assumes  penetration 


proportional  to  the  projectile  velocity.  The  equation  was 
of  the  form 


F  * 


k,  V  -  ^ 


(4.27) 


the  yield  being  included  in  the  last  term.  was  developed 

by  considering  the  movement  of  a  cylindrical  core  of  material 
through  the  target  plate  under  the  area  of  contact  involved 
during  the  collision.  This  cylinder  is  restrained  laterally, 
but  is  assumed  to  be  free  to  move  in  the  direction  of  the 
impact  for  the  time  taken  by  a  stress  wave  to  make  a  round  trip 
through  the  projectile,  i„e„  2d,/Cp  .  For  an  elastic  colli¬ 
sion,  this  cylindrical  core  within  the  target  will  traverse 
a  depth 


fir 


1/ 


c*  A  r  cr  fir 


Cppp  and  ctPt  being  the  acoustic  impedances  of  the  pro¬ 
jectile  and  target  materials,  respectively.  Experimental 
data  for  a  number  of  target  materials  impacted  by  mercury 
and  water  drops  indicated  the  need  for  a  correction  factor 
of  3.6.  Thus,  the  equation  is 
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The  intercept  of  this  equation  with  the  velocity 
axis  represents  the  impact  speed  below  which  no  indentation 
would  occur..  This  speed  appeared  to  be  a  function  of  the 
target  strength.  Engel  found  that  it  could  be  accounted  for 
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by  the  dimensionless  product 
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yielding  the  final  equation 


ii  »  2*.  Cp  -fe. _  (--)  -  !3i'Q  -2±2/j.  .1  (4o29) 

^  s  f**-crfp  '  crl  r2  ft  ^ 

where  *  136 8  and  is  the  dynamic  yield  strength 

of  the  target o  Note  that  the  targets  used  for  the  evaluation 
of  ki  and  k„  had  thicknesses  of  between  1,5  times  and  5 
times  the  diameter  of  the  impinging  projectile,  the  back  of 
the  target  being  unsupported „ 


4.3  Dimensional  Analysis 

In  carrying  out  a  dimensional  analysis  of  the 
high-velocity  cratering  problem,  the  main  difficulty  lies  in 
distinguishing  the  main  features  of  the  phenomena,  and  defining 
appropriate  physical  quantities  describing  the  behavior  of  the 
materials  involved „ 

The  geometric  quantities  of  importance  are: 

d  relevant  dimension  of  the  projectile 

p  relevant  dimension  of  the  crater 

V  projecLile  velocity  relative  to  the  target 

$  angle  of  incidence  of  projectile  on  the 
target  surface 

T  target  thickness,  or  other  relevant  target 
dimension 


51 


Where  projectiles  of  different  shapes  are  involved, 

the  relevant  projectile  dimension  is  not  obvious..  Several 

investigators  have  used  different  systems „  Shaped  charge  jet 

penetration  has  been  correlated  on  the  basis  of  the  length  of 

27 

the  jeto  Summers  and  Charters  similarly  used  the  length 
of  the  projectile  in  the  direction  of  initial  projectile 
motion v  Others  have  in  effect  used  the  cube  root  of  the 
mass,  to  define  an  equivalent  diameter 

-  -tt; 

Jfi 

Where  spheres  and  cylinders  only  are  compared, 
the  diameter  of  an  equivalent  sphere  of  mass  equal  to  that 
of  the  cylinder  has  sometimes  been  used  i„e. 

*te  -  (i  (4 -«) 

Where  d  is  the  cylinder  diameter  and  L  its 
128 

length,,  Bjork  used  arbitrarily 

Je  =  (JXL)‘'3  (4-32) 

to  compare  his  theory  for  cylindrical  projectiles  with  experi¬ 
ments  conducted  with  spherical  projectiles,, 

As  long  as  the  crater  shape  remains  identical  as 
other  parameters  are  varied.,  the  crater  size  may  be  repre¬ 
sented  by  a  single  parameter  Throughout  the  range  explored 
experimentally,  however,  the  crater  shape  varies  from  a  cylin¬ 
drical  cavity  of  the  same  diameter  as  the  projectile  at  very 
low  velocities,  to  an  almost  hemispherical  crater  at  the 


highest  velocities  examined.  Of  prime  importance  in 
practice  is  the  penetration,  or  depth  of  crater,  and 
this  parameter  has  been  used  most  in  empirical  correlations 
of  cratering  in  thick  targets  Other  parameters  of  interest 
are  the  crater  volume,  and  diameter  of  the  crater  in  the 
plane  of  the  original  surface. 

Quantities  describing  the  relevant  material 

properties  are  not  so  straightforward „  Cratering  involves 

shock  compression  of  the  materials  to  very  high  pressures, 

possibly  accompanied  by  melting  or  vaporization.  Subsequent 

motion  involves  both  extremely  high  stresses  and  strain 

rates.  The  final  phase  of  the  motion  has  been  shown  to 

7 

involve  elastic  springback  of  the  crater. 

Parameters  describing  the  behavior  of  solids  under 
very  high  stresses  and  strain  rates  have  not  been  adequately 
formulated.  Equations  of  state  of  solids  under  extremely 

high  pressures  have  received  considerable  attention  recently. 
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From  the  work  oi  the  group  at  Los  Alamos  and  else¬ 

where,  it  appears  that  the  shock  Hugoniot  may  be  described 
adequately  by  only  two  parameters  in  addition  to  the  density 
at  zero  pressure-  For  a  number  of  metals,  one  parameter  co¬ 
incides  with  the  adiabatic  sound  speed.  The  physical  inter¬ 
pretation  of  the  other  parameter  is  not  clear,  but  it  appears 
t:o  be  related  to  the  change  of  compressibility  with  pressure 

The  parameters  necessary  to  define  the  strength  of 
the  material  are  also  not  clear..  The  conventional  engineer¬ 
ing  stress-strain  properties  are  not  related  in  a  simple  way 
to  the  stress-strain  properties  of  the  material  under  complex 
combined  stresses  at  high  strain  rates.  Nevertheless,  draw¬ 
ing  on  properties  which  are  defined  under  static  conditions, 
one  has  the  relevant  elastic  modulus  G  and  yield  stress  Y 
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In  addition  there  would  be  quantities  describing  the  strain 
hardening  and  strain  rate  behavior * 

One  quantity  which  has  been  widely  measured,,  but 
which  is  difficult  to  relate  to  other  physical  properties, 
is  the  indentation  hardness.  Its  ease  of  measurement  to¬ 
gether  with  the  fact  that  the  process  of  pushing  a  ball  or 
pyramid  into  a  solid  surface  is  suggestive  of  the  cratering 
process  has  nevertheless  led  to  its  use  in  empirical  expres¬ 
sions.  The  Brinell  Hardness  number  H  is  strictly  defined 
as  the  load  used  divided  by  the  surface  area  of  the  indenta¬ 
tion  produced  by  a  spherical  stylus,  and  thus  has  dimensions 
of  a  stress,  with  units  kilograms  per  square  millimeter. 

Not  all  of  these  elementary  material  properties 
are  independent.  For  example,  the  sound  speed  is  related 
to  the  density  and  relevant  elastic  modulus” 

2  JL 
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There  appears  t.o  be  a  relation  between  Brinell 
Hardness  and  yield  stress.  In  Fig  4.1  Brinell  Hardness  has 
been  plotted  versus  yield  stress  for  a  number  of  representa¬ 
tive  target  materials,  and  it  is  seen  that  Brinell  Hardness 
is  roughly  proportional  to  yield  stress,  at  least,  within 
the  scatter  usually  associated  with  empirical  correlations. 


Using  the  independent  elementary  material  properties 
P  ,  c  ,  and  H  together  with  the  geometric  quantities  allows 
the  following  dimensionless  parameters  to  be  formed. 
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Penetration  ratio 


P 


Density  ratio 
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Mach  number 
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Best  number 


The  last  parameter  appears  frequently  in  one 
or  another  of  its  equivalent  forms  in  both  the  empirical  and 
theoretical  cratering  equations,  and  it  appears  appropriate 
to  give  it  a  name.  Following  a  suggestion  of  Dr.  John  S. 
Rinehart,  it  is  proposed  here  to  call  this  the  Best  number, 
after  the  French  ballistician  Best.  Between  1835  and  1845 
the  French  carried  out  a  series  of  experiments  to  measure 
the  velocity  of  fragments  by  putting  hollow  cannon  balls 
loaded  with  powder  into  a  well  and  detonating  them.  At 
the  bottom  of  the  well  was  a  1  1/2  foot  layer  of  damp  clay. 
They  measured  the  volume  of  the  holes  made  in  the  clay  by 
the  fragments.  In  order  to  calibrate  the  clay,  they  fired 
a  pistol  into  the  clay  alongside  the  fragment  hole,.  Knowing 
the  velocity  of  the  pistol  ball,  and  its  mass,  and  the  volume 
of  the  hole  produced  together  with  the  fragment  mass  and  its 
hole  volume,  the  fragment  velocity  could  be  calculated  by 
assuming  that  the  cratering  efficiency  remained  constant,, 

The  cratering  efficiency  is  defined  as  the  crater 
volume  divided  by  the  kinetic  energy  of  the  projectile,.  It 
has  been  noted  in  some  experiments  that  the  cratering  effi¬ 
ciency  is  proportional  to  target  hardness..  Under  this 
assumption 

I  =  k  H 
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or  rearranging,,  by  noting  that  the  crater  volume  is  proportional 
to  the  cube  of  a  linear  dimension  of  the  crater  if  the  crater 
shape  remains  constant. 
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Thus,  it  appears  particularly  apt  to  apply  the 

name  Best  number  to  the  parameter  /pv£\.. 

V  H  / 

If  melting  or  vaporization  are  considered  important, 
then  the  relevant  quantity  wi 11  be  the  energy  to  melt  or 
vaporize  unit  mass  of  material..  This  will  be  related  to 
specific  and  latent  heat..  Although  melting  or  vaporization 
are  a  result  of  shock  heating,  and  may  be  considered  to 
occur  at  some  time  during  the  high-pressure  phase  of  the 
motion,  the  material  is  finally  ejected  at  ambient  (zero) 
pressure,  and  thus  the  energy  carried  away  by  melted  or 
vapoi'ized  material  may  be  described  by  zero  pressure  data 
on  specific  and  latent  heats. 

If  the  energy  necessary  to  sublimate  or  melt  unit 
mass  of  material  is  denoted  by  Q  ,  an  additional  non-dimen¬ 
sional  parameter  may  be  formed,  i.. e 


Thermal  parameter 


It  is  not  possible  at  this  time  to  form  other 
dimensionless  parameters  based  on  compressibility,  strain 


i 

1 


hardening,  or  strain  rate  effects,  since  suitable  material 
properties  have  not  been  adequately  defined. 

Examining  the  various  empirical  and  theoretical 
equations,  it  may  be  seen  that  the  equations  which  are 
dimensionally  correct  can  be  reduced  to  relations  between 
the  dimensionless  parameters  presented  above. 

It  is  evident  that  the  density  ratio  (p  /p.)  and 

pc 

Mach  number  (V/c^)  have  been  used  in  numerous  empirical 
expressions . 

By  noting  that 

>>  -  k  /V 

the  expressions  used  at  BRL  and  NRL  Eq»  (4,26)  may  be 
reduced  to 


(4,.  33) 


thus  involving  the  Best  number,  while  Pugh  and  Ei ohelberger ' s ^ 
expression,  Eq.  (4.25),  is  very  nearly  equivalent  to 


il  _ 

d  _ 


(4.34) 


except  for  the  exponent  in  Brinell  Hardness.. 

The  expression  used  by  Engel^0.,  Eq„  (4.29),  in¬ 
volves  a  variation  of  the  Best  number .  It  may  be  written 


k, 
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r>>  fp 


(4  35) 
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In  Collins  and  Kinard's"^  expression,  Eq,  (4„24) 


»  =  k, 

r  „  ,  .  \  6-7Q 


(4o36) 


the  constants  are  evidently  not  dimensionless.  Terms  of  the 
form 
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are  involved, 
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By  rewriting  Bohn  and  Fuchs*  theoretical  expression, 
Eq,  (3,21),  slightly,  it  may  be  seen  that  they  also  use  the  Best 
number. 


(4,37) 

It  will  be  remembered  that  shaped  charge  jet  pene¬ 


tration  results,  Eq,  (3,36)  were  extensively  correlated  with 
the  use  of  the  Best  number,  i,e. 


JL 

l 


(4  38) 
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While  no  empirical  fits  have  been  previously 
attempted  using  thermal  quantities,  the  thermal  penetration 
theories  (Section.  3  =  3)  all  involve  only  the  parameters 


fp 

ft 
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4 « 4  Comparison  of  Theories  and  Empirical  Expressions  on 
the  Basis  of  Resistance  to  Penetration 


Some  of  the  penetration  theories  and  empirical 
expressions  were  devised  by  either  deducing  or  assuming  the 
dependence  of  the  force  resisting  penetration  on  the  physical 
parameter s,  and  solving  the  resultant  equation  of  motion 


X  e 

6  Jr 


l-  '  c\y  =  -  F"  .  A  r. 


(4.38) 


where  F  is  the  resistance  force. 

It  is  very  iixteresting  to  compare  these  theories 
and  empirical  expressions  on  the  basis  of  the  resistance  to 
penetration  which  they  implicitly  assume.  For  those  empiri¬ 
cal  expi'essions  which  were  not  derived  by  assuming  the  form 
of  the  expression  for  the  resistance.,  F  may  be  obtained 
simply  by  differentiating  the  expression  for  penetration 


4.4.1  Resistance  Independent  of  Velocity 

t  r  >  i  ,  •  p.  .  •».,«** 

jl.  4»  cut:  -.'JtioLuutfj  auti  Uiiifccci  uj  wu  eilt  O  L 

velocity,  and  simply  dependent  on  the  presented  area  of  the 
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projectile. 


F  =  k  cl'1' 


(4,39) 


which  is  equivalent  to  the  Euler-Robins  penetration  law, 

Eq,  (4,2),  Adding  a  constant  to  allow  for  free  surface 
effects  yields  the  NRL  expression,  Eq,  (4,3)  „ 

f  =  K  fr  V%  +  (4,40) 

If  free  surface  effects  are  taken  into  account 
by  assuming  that  the  resistance  depends  on  the  depth  of 
penetration, 

F-  kp.cl  Z  ~  K  fr  V'  (4,41) 


which  is  eqxiivalent  to  the  expression  introduced  by  Fairbairn, 
Eq..  (4,8),  Assuming  a  stronger  dependence  of  resistance  on 
depth  of  penetration 


F  -  k.Px 


Vi 


(4,42) 


which  is  the  expression  proposed  by  Rinehart,  Eq,  (4.16), 
and  which  is  in  use  at  the  Naval  Research  Laboratory  and  the 
Ballistics  Research  Laboratory,  Eq,  (4,26),  A  slight  varia¬ 
tion  yields 


t  - 


(4,43) 


which  may  be  recognized  as  the  expression  introduced  by 
Charters,  Eq.  (.4 .22) 
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4„4.2  Resistance  Dependent  on  Velocity 


If  the  resistance  is  assumed  to  be  directly 
proportional  to  velocity., 

F  -  k  dxv  f  =  K  pr  V  (4.44) 


An  expression  of  this  type.,  with  the  addition  of  a  constant 
to  allow  for  a  minimum  velocity  to  produce  finite  penetration 
has  been  used  by  Kinard  and  Collins.,  Eq.  (4.24),  and  by 
Engel,  Eq°  (4«28) .  By  allowing  for  an  increase  of  resistance 
with  penetration  depth 


r  --  k  p'v 


(4. .45) 


which  is  the  expression  proposed  by  Pugh  and  Eichelberger, 
Eq  o  (4.17). 

Expressions  in  effect  assuming  an  inverse  depen¬ 
dence  on  velocity  are  the  Tresidder  equation,  Eq  (4.9) 


P.d 

F  -  k  ™ 


J  *  Kfi> 


(4.46) 


and  the  DeMarre  expression,  Eq..  (4.13) 


7  *  K  Pr  Y  <*-«) 


F  =  K  „ 
Jr 


It  may  be  noted  that  all  the  expressions  in  this 
subsection  may  also  be  derived  by  assuming  resistance  inde¬ 
pendent  of  velocity,  but  dependent  on  some  power  of  the 
Denetr.af -Ion  Thus  Kinard  and  Cwlliuo’  amu  Engel's 
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expressions  also  imply 
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Fugh  and  Eichelberger 1 s  expression  implies 
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Tresidder's  equation  yields 


r  -  /<  (/>, V)6 
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and  the  DeMarre  expression  yields 
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(4.43) 


(4..  49) 


(4.50) 


(4,51) 


4,4,3  Resistance  Dependent  on  Inertia 

If  the  resistance  is  assumed  to  be  proportional 
to  the  inertia  forces  exerted  by  the  target ,  the  resistance 
force  may  be  written 

F  =  k  <)'//  (4  52) 


In  this  case  integration  of  the  equation  of  motion 
cl  J  /  cli/  ~  -  k  cl  p,  V  d/>  (4,53) 
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leads  to  the  expression 


J  -  K  U%  V  +  ^  (4o54) 

v/hich  is  the  form  of  expression  used  by  Grimminger,  Eq.  (3»17)„ 
If  the  resistance  is  assumed  to  be  proportional  to  an  inertia 
term,  and  a  term  independent,  of  velocity  to  allow  for  material 
strength,  the  resistance  becomes 

F  =  dX  (  ^  *  fi  yX)  (4.55) 

This  is  exactly  the  type  of  resistance  term  found 
ixi  the  armor  penetration  theories  of  Bethe,  and  Thomson, 
(Section  3.1),  who  derived  their  results  by  solving  the 
dynamic  elastic-plastic  equations  of  motion  of  the  target,, 

Their  resistance  expressions  are  of  the  form 

F  *  !<,  dx  (KY  *  vl)  (4 .56) 


See  Eq„  (3«7),  Eq„  (3. .9)  and  Eq.  (3„ 11)«  Substituting  the 
more  convenient  Brineli  hardness  for  the  yield  strength, 
since  these  have  been  shown  to  be  roughly  proportional,  the 
equation  of  motion  becomes 


fp  -  -  k.d'(K  H  *  ft*')  dr 


When  integrated  this  yields  the  equation 


p 

rl 


/ 

4  k, 


(4.57) 


(4  58) 
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Bohn  and  Fuchs  (Section  3.1.6)  assumed  a  resistance 
force  of  the  form 

F  =  k,  (  Jk  til  +  Jfr  v  )  (4.59) 


Integration  of  the  equation  of  motion  then  led  to  Eq.  (3.20) 


4„4o4  Discussion 


It  may  be  seen  from  the  above,  that  the  penetration 
laws  in  use  may  be  divided  into  two  general  categories.  In 
the  first,  the  assumption  is  made  that  the  resistance  varies 
with  some  power  of  the  distance  from  the  surface.  This 
assumption  then  leads  to  a  penetration  law  of  the  general 
form 


£ 

A 


(4.61) 


If  the  resistance  is  assumed  to  increase  more  rapidly  with 
increasing  distance  from  the  surface  (higher  power  of  p  ) , 
the  resultant  velocity  exponent  n  decreases.. 

In  the  second  category,  the  resistance  is  assumed 

to  be  due  to  the  inertia  of  the  target,  and  therefore  prop-or- 

2  * 
tloxial  to  v  -  The  resultant  dependence  of  penetration  on 
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velocity  is  logarithmic .  When  a  resistance  term 
independent  of  velocity  is  included,  the  resultant 
equation  has  some  very  interesting  properties,  see 
Eq,  (4,58) .  At  low  velocities,  the  contribution  of  the 
strength  term  is  dominant,  i,e,  in  the  neighborhood 


However,  as  the  velocity  increases,  the  strength  term 
becomes  increasingly  less  important,  so  that  at 


Eq,  (4„63)  becomes  very  nearly 


% 


(4,62) 


It  is  clear  that  the  strength  term  now  enters 
as  a  constant  correction,  much  in  the  way  that  it  was 
included  in  Grimminger ' s  theory.  Not  until 


will  the  effect  of  strength  become  negligible. 


It  may  also  be  noted  that  the  logarithmic  function 
exhibits  a  decreasing  slope  with  increasing  velocity.  In 
order  to  i^nrnpf’r113  t-ho  lopari  f  hmic  form  with  the  previous 
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simple  power  law,  Eq.  (4,58)  has  been  plotted  on  log- log 
paper  in  Fig. 4 ,2  denoting  ./Pt ■=  a  .  In  this  way 
the  slope  of  the  curve  will  correspond  to  the  exponent  n 
in  the  expression  p  ■=  kVn  . 

It  is  seen  that  the  slope  decreases  from  about  2 
at  otV  «=  0,3  to  1/6  at  aV  *»  100  „  The  slope  is  approxi¬ 
mately  2/3  in  the  range  3  <  aV  <  6  .  Eq.  (4.54)  ana  Eq.  (4.60) 
are  also  plotted  in  Fig. 4.2,  the  curves  having  been  adjusted 
to  coincide  at  aV  =■  4  ,  i.e.  in  the  velocity  range  where  the 
curves  are  approximated  by  p  «  kV  '  .  Grimminger ' s  form 

Eq.  (4.54)  agrees  well  with  the  previous  curve  at  higher 
velocities,  but  departs  drastically  at  lower  velocities.  How¬ 
ever  Bohn  and  Fuchs'  expression  Eq.  (4.60)  agrees  very  well 
at  lower  velocities,  but  departs  at  higher  velocities,  the 
slope  decreasing  more  slowly,  and  approximating  1/3  in  the 
range  20  -c  aV  <  200  „ 

It  is  possible  to  deduce  approximate  velocity 
ranges  over  which  the  curves  of  Fig..  4.2  may  be  approxi¬ 
mated  by  p  =  kVn  for  different  values  of  n  ,  for  a 
typical  target  material.  For  aluminum  targets  for  example, 
it  is  observed  that  penetration  follows  the  law  p  =  kV^^ 
above  about  10,000  ft/sec.  The  curves  of  Fig.  4?.2  first 
show  agreement  with  this  law  at  about  aV  «  3  .  Thus  a 
may  be  approximately  evaluated,  and  corresponding  velocity 
ranges  found. 

These  are  shown  In  the  table- 
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Approximate  Velocity  Ranges  in  ft/sece 
in  which  Penetration  May  be  Approximate  by  p  <=  kV 


Velocity 
Exponent  n 

Eq»  (4»54) 

Eq.  (4,58) 

Eq«  (4,60) 

v2 

4,000  -  6,000 

1200 

1200 

VV3 

6,000  -  10,000 

2,000  -  6,000 

2,000  -  6,000 

V 

7,500  -  12,000 

5,000  -  10,000 

5,000  -  10,000 

V2/3 

10,000  -  20,000 

10,000  -  20,000 

10,000  -  35,000 

vl/3 

30,000  -  200,000 

30,000  -  200,000 

60,000  -  300,000 

vV6 

600,000 

600,000 
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SECTION  V 

NORMAL  IMPACT  ON  QUASI -INFINITE  TARGETS 


This  section  will  be  devoted  to  a  discussion  of 
the  applicability  of  the  theories  presented  in  Section  III 
and  the  empirical  relationships  presented  in  Section  IV  to 
cratering  in  quasi-infinite  targets  by  projectiles  striking 
the  target  at  normal  incidence.  The  cratering  process  is  an 
extremely  complex  one,  involving  many  different  types  of 
material  fcehavior.  It  is  therefore  necessary  to  first  give 
a  description  of  the  cratering  process.  The  applicable 
theories  are  then  compared  to  the  appropriate  experimental 
data  and  the  fits  of  empirical  expressions  discussed  in 
Section  IV  to  the  data  are  then  examined. 

5  - 1  Description  of  Phenomena 

Several  excellent  studies  of  cratering  by  high 
velocity  projectiles  have  recently  been  carried  out  in 
which  high  speed  photography,  X-radiography  and  other 
techniques  were  used  to  make  active  observations  of  pheno¬ 
mena  accompanying  the  penetration  process.  When  taken 
together  with  post  mortem  observations  of  the  crater  and 
adjacent  material,  and  interpreted  in  the  light  of  present 
knowledge  of  the  behaviour  of  solids  under  high  pressures, 
a  fairly  detailed  qualitative  understanding  of  the  physical 
phenomena  involved  in  cratering  can  be  gained.  Several 
descriptions  of  these  phenomena  have  recently  been  given.  3 
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At  the  instant  of  first  contact  of  the  projectile 
with  the  target,,  local  conditions  at  the  projectile- target 
interface  may  be  estimated  by  one- dimensional  shock  theory > 
Such  theories  both  for  very  high  pressures*  and  for  low 
stresses,  have  recently  received  a  great  deal  of  attention 
The  Hugoniot  relations  for  the  shock  proceeding  into  the 
target  and  the  shock  receeding  back  into  the  projectile 
determine  the  initial  shock  velocities  and  interface  velo¬ 
city  and  pressure..  (See  Appendix  D»)  Interface  pressures 
and  velocities  calculated  on  this  basis.,  and  for  the  acoustic 
approximation,,  i,e«  neglecting  the  change  of  compressibility 
with  pressure  of  the  materials.,  are  plotted  as  a  function  of 
projectile  velocity  in  Fig,  1>  ,  la  and  b  for  several  common 
materials.  It.  may  be.  noted  that  the  acoustic  approximation 
leads  to  very  large  errors  at  even  modest  velocities. 
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The  shock  receeding  into  the  oncoming  projectile 
will  actually  be  carried  below  the  original  target  surface 
when  the  initial  projectile  velocity  exceeds  the  velocity 
of  the  shock  wave  generated  in  the  projectile.  The  criti¬ 
cal  velocity.,  above  which  the  shock  in  the  projectile  is 
carried  below  the  original  target  surface,.  Is  a  function 
of  the  densities  and  compressibilities  of  the  projectile 
and  target  materials.,  see  Appendix  D,  Critical  velocities 
for  several  common  material  combinations  are  shown  below 


Minimum  Projectile  Velocities  for  the  Shock  to 
Move  Below  the  Original  Target  Surface 


'"--..Target 

Project! le  \ 

A1 

Cu 

Pb 

A1 

43 

73 

79 

Cu 

22 

37 

38 

Pb 

12 

20 

21 

velocities  in  thousands  or  tect  per  second 
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The  initially  one “dimensional  motion  is  of  course 
localised  around  the  nose  of  the  projectile,  and  is  almost 
immediately  altered  by  the  geometry  of  the  subsequent  motion. 
Nevertheless  the  above  considerations  are  useful  for  esti¬ 
mating  maximum  pressures,  temperatures  and  shock  velocities 
resulting  from  the  impact. 

At  low  impact  velocities  of  a  few  thousand 
feet  per  second,  initial  interface  pressures  are  of  the 
order  of  the  dynamic  yield  strengths  of  typical  materials, 
as  may  be  seen  from  Fig.  5.1a.  It  is  therefore  clear  that 
the  motion  will  be  strongly  influenced  by  material  strength. 

In  cases  where  the  projectile  strength  is  much 
greater  than  the  target  strength,  as  in  the  impact  of  a 
hard  copper  projectile  on  a  soft  lead  target,  the  pro¬ 
jectile  may  suffer  only  minor  deformation  at  low  impact 
velocities,  and  an  almost  cylindrical  crater  ol  diameter 
equal  to  that  of  the  projectile  is  formed  by  a  process  of 
dynamic  plastic  deformation  of  the  target  material.  As 
the  initial  projectile  velocity  increases,  the  projectile 
suffers  increasing  plastic  deformation.  The  inertia  effects 
also  lead  to  increased  radial  motion  of  the  target  material, 
so  that  the  resultant  crater  is  of  Increasingly  greater 
diameter  than  the  original  projectile  diameter,  until  a 
velocity  is  reached  where  the  inertia  effects  piredominate, 
and  an  almost  hemispherical  crater  is  formed. 

When  the  proiecfile  is  b-riftlR,  a  iienarturo  from 
the  above  behaviour  is  observed.  The  shock  wave  receeding 
into  the  projectile  on  reflection  at  the  free  surfaces  of 
the  projectile  may  lead  to  stress  concentrations  which  are 
sufficient  to  cause  brittle  fracture  at  velocities  where 
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very  little  plastic  deformation  of  the  projectile  occurs.. 
Tungsten  carbide  and  hardened  steel  projectiles  exhibit 
this  type  of  behaviour,,  Just  above  the  velocity  necessary 
to  cause  projectile  fracture,  the  projectile  is  observed 
to  break  info  several  large  fragments.  The  presented  area, 
of  the  fragments  is  considerably  greater  than  that  of  the 
unbroken  projectile,  and  a  reduction  of  penetration  may 
result,,  although  the  crater  volume  continues  to  increase. 

At  increased  the  projectile  suffers  increasing 

fragmentation,,  and  the  penetration  again  increases.  Finally., 
when  inertia  effects  predominate.,  a  nearly  hemispherical 
crater  again  results., 

These  effects  are  illustrated  in  Fig,,  5.2  which 
shows  craters  formed  by  lead,  copper  and  hardened  steel 
spheres  in  soft  lead  targets  at  several  identical  initial 
projectile  velocities.  Figure  .5,3  shows  the  penetration  as 
a  function  of  velocity  for  typical  ductile  projectile  be¬ 
haviour,,  while  Fig,  5,4  shows  the  penetration  as  a  function 
of  velocity  tor?  typical  brittle  projectile  behaviour., 

Three  regions  of  impact  are  commonly  defined; - 

a)  Low1  velocity  region  in  which  the  projectile  suffers 
only  minor  deformation,  and  which  has  also  been  termed  the 
armor  penetration  region  or  undeformed  projectile  region, 

b)  transition  region,  c.)  high  velocity  region  where  nearly 
hemispherical  craters  are  formed,  and  which  has  also  been 
termed  the  hypervelocity  region,  or  the  '’fluid  impact" 
region,,  from  the  fact  that  inertia  effects  predominate. 

It  is  interesting  to  note  that,  the  penetration 
appears  to  be  nearly  proportional  to  the  4/3  power  of 
velocity  in  the  low  velocity  region,,  while  the  penetration 


is  approximately  proportional  to  the  2/3  power  of  velocity 
in  the  high  velocity  region.  In  the  case  of  ductile  pro¬ 
jectiles,  the  high  and  low  velocity  regions  are  connected 
by  a  smooth  curve  in  the  transition  region..  Brittle  pro¬ 
jectiles  show  an  extended  lev/  velocity  region  with  the  drop 
in  penetration.,  discussed  previously,  in  the  transition 
region . 


At  high  impact  velocities,  the  initial  pressures 
are  very  much  greater  than  the  dynamic  yield  strengths  of 
the  materials,  as  may  be  seen  from  Fig.  5.1a,  and  the 
strength  of  the  materials  may  well  be  neglected  at  early 
stages  of  the  cratering  process.  The  projectile  and  target 
materials  are  in  fact  observed  to  flow  very  much  like  fluids 
in  this  region.  For  example,  hemispherical  craters  are 
often  observed  to  be  plated  with  a  thin  layer  of  projectile 
material.  Fig.  5.5.  At  still  higher  velocities  all  the  pro¬ 
jectile  material  apparently  flows  out  of  the  crater. 


As  the  projectile  approaches  the  target,  the 
projectile  and  target  surfaces  will  generally  meet  at 
a  small  angle.  High  velocity  collision  oi  surfaces  at  small 
angles  leads  to  the  phenomenon  of  jetting^^,  material  being 
expelled  along  the  bisector  of  the  angle  between  the  surfaces 
at  considerably  higher  velocities  than  the  closing  velocity 
of  the  surfaces.  The  expelled  particles  may  move  with 
sufficient  velocity  so  that  ablation  occurs  and  the  ablated 
material  burns  if  air  is  present.  The  resultant  luniino- 
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has  been  studied  in  some  detail  by  Grow  et 
7  7  H 1 

Cook'  ’  has  photographically  observed  areas  interpreted 
to  be  clouds  of  metal  vapour  emanating  .from  the  impact  area, 
which  very  likely  are  due  to  the  same  mechanism. 
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As  the  projectile  eontirmes  to  move  into  the 
target,,  the  shock  wave  in  the  target  precedes  the  projectile- 
target  interface.  Rarefactions  from  the  free  .surface  of  the 
target  and  projectile  modify  the  shock  system,  and  the  shock 
rapidly  becomes  nearly  hemispherical »  High  speed  framing 
camera  sequences  of  penetration  in  transparent  plastics 
(Fig.  5.6}  and  flash  X- radiographs  sequences  of  penetration 
in  metallic  targets  (Fig.  5.7)  clearly  show  this  behaviour. 

At  the  target  surface  adjacent  to  the  projectile,  the 
high  pressure  target  material  is  unsupported  and  is  forced 
out  of  the  crater  to  form  a,  spray.  Several  fine  X- radio¬ 
graphs  and  high  speed  framing  camera  sequences  of  the  spray 
have  been  made,,  (Fig.  5.8).  It  is  clear  that  a  large  pro¬ 
portion  of  the  material  originally  in  the  crater  is  expelled 
in  the  spray.,  which  moves  at  relatively  low  velocities. 

Numerical  solution  of  the  differential  equations 
of  motion  of  an  inviscld  fluid,,  with  appropriate  equation 
of  state,,  show  shock  and  rarefaction  systems  closely  re¬ 
sembling  those,  observed  experimentally,,  and  considerably 

clarify  the  understanding  of  the  cratering  process.  Results 
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obtained  by  Rjork  in  the  form  of  pressure  contours 

and  velocity  vectors  for  the  impact  of  a  cylindrical  pro¬ 
jectile  on  a  semi- infinite  target  are  shown  in  Fig.  5.9 
and  5,10.  (It  should  be  noted  that  the  numerical  method 
uses  an  artificial  viscosity  term  to  render  the  solution 
continuous.,  and  shocks  are  therefore  smeared  over  a  finite 
distance  instead  of  appearing  as  discontinuities..  The 
shocks  are  smeared  over  a  fairly  large  distance  in  Bjork's 
solution  since  a  rather  coarse  finite- difference  space 
mesh  was  used  tor  reasons  of,  computer  time  economy.) 
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Figure  5 <>9  shows  results  for  the  impact  of  an  iron  projectile 
on  an  iron  target  at  18 <,000  ft/sec «  The  shock  in  the  projec¬ 
tile  is  seen  to  rise  above  the  original  surface  of  the  target. 
Figure  5.10  shows  results  for  the  impact  of  an  iron  projectile 
on  rock  (tuff)  at  100,000  ft/sec.  The  shock  in  the  projectile 
is  carried  below  the  original  surface  of  the  target*  and  the 
resultant  differences  in  the  flow  pattern  can  be  seen  when 
compared  to  Fig.  5.9. 

If  the  shock  in  the  projectile  is  carried  below 
the  original  surface  of  the  target*  the  projectile  material 
which  has  not  yet  been  reached  by  the  shock  is  carried  well 
below  the  target  surface.,  and  one  may  expect  much  greater 
penetration.  On  the  other  hand*  if  the  shock  in  the  pro¬ 
jectile  receeds  above  the  target  surface*  the  sides  of  the 
shock  zone  in  the  projectile  are  unsupported,,  and  increased 
lateral  flow  may  be  expected,,  with  consequent  shallower 
craters.  In  an  extreme  case*  such  as  the  impact  of  nylon  on 
aluminum.,  the  stresses  in  the  unsupported  shocked  zone  in 
the  projectile  may  be  sufficiently  high  to  cause  the  pro¬ 
jectile  to  flow  out  over  the  surface  of  the  target*  and 
craters  considerably  shallower  than  hemispheres  have  been 
observed  under  such  conditions. 

While  the  initial  phase  of  the  motion  resulting 
from  high  velocity  impact  may  be  adequately  described  by 
hydrodynamic  principles*  the  stresses  rapidly  decay  due  to 
geometrical  divergence  and  dissipation  to  the  point  where 
material  strength  becomes  important.  The  material  being 
ejected  from  the  crater  will  now  remain  attached  to  form  a 
raised  lip.  If  the  target  material  is  brittle, portions  of 
this  lip  may  fracture  and  be  removed.,  The  final  phase  of  the 
motion  will  involve  a  certain  amount  of  elastic  springback. 
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Gehring^ ,  using  flash  X-radiography  has  estimated  that 
the  elastic  recovery  in  crater  volume  may  be  as  high  as  3Q7o 
in  an  aluminum  target  struck  by  a  steel  projectile  at  lb, 400 
ft/ sec  „ 

Material  which  has  been  subjected  to  very  intense 
shock  compression,  followed  by  an  adiabatic  expansion  is 
returned  to  zero  pressure  at  higher  temperature  than  ambient, 
due  to  the  entropy  increase  in  the  shock  process.,  Under 
conditions  of  high  velocity  impact ,  this  temperature  may 
be  above  the  melting  point,.  For  impacts  below  20,000  ft./sec, 
theoretical  estimates  indicate  that  a.  small  zone  of  target 
material  near  the  initial  point  of  impact  may  be  melted  in 
normal  metallic  materials,.  In  low  melting  point  materials 
such  as  lead,  melting  may  be  much  more  extensive.  Since  the 
melted  material  has  no  strength  on  return  to  zero  pressure, 
it  will  be  ejected  from  the  forming  crater  somewhat  more  readily 
than  it  it  had  not  been  melted.,  This  mechanism  may  there¬ 
fore  contribute  to  the  cratering  process  at  high  impact 
velocities..  The  question  whether  the  material  melts  during 
the  high  pressure  phase  of  the  motion  is  academic,,  since  the 
mbin.  effect  wi  11  b®  an  absorpfc ist*.  oi  latent  heat  >.  "which. . 
is  negligible*  Theoretical  estimates  for  minimum  shock 
pressures  necessary  to  cause  subsequent  melting  on  return 
to  zero  pressure  for  typical  materials  are  7.4  ki. lobar s  for 
lead,  223  ki lobar s  for  tin,  325  kilobars  for  cadmium,  and 
1.4  megabars  for  copper. 

The  material  immediately  adjacent  to  the  crater 
sometimes  shows  signs  of  recrystallisation  at  the  highest 
experimental  impact  velocities,  confirming  that  local 
melting  has  occurred.  More  remote  material  usually  shows 


evidences  of  severe  plastic  deformation  in  severely  deformed 
crystal  structure  or  shock  twinning,  and  in  increased  indenta¬ 
tion  hardness,  these  effects  being  reduced  at  greater  dis¬ 
tances  from  the  crater <>  Examples  are  shown  in  Figo  5 oil. 


5.2  Comparison  of  the  Data  with  Theory 

The  penetration  data  for  quasi- infinite  targets 
for  which  sufficient  data  points  exist  have  been  presented 
in  Fig.  1.2  through  1.25.  (See  Section  II.)  The  plots 
represent  data  for  impact  of  a  variety  of  projectile 
materials  on  aluminum,  steel,  copper  and  lead  targets. 

The  penetration  has  been  normalised  with  respect  to  pro¬ 
jectile  diameter  in  these  plots.  For  cylindrical  projec¬ 
tiles,  the  diameter  of  an  equivalent  sphere  of  equal  mass 
was  used,  i.e.  dg  «=  (3/2  .  Drawn  on  these 

figures  are  curves  representing  theoretical  predictions 
of  penetration  according  to  the  applicable  theories  pre¬ 
sented  in  Section  III.  Although  some  of  these  theories 
were  intended  to  apply  to  a  much  higher  velocity  range,  they 
have  nevertheless  been  computed  in  the  experimental  velocity 
range,  and  included  in  the  plots.  Each  of  the  theories  will 
be  briefly  discussed. 

5.2.1  Rigid  Projectile 


The  applicable  theories  in  this  class  are 
Griuuninger ' s  theory  (Section  3.1.5)  and  Bohn  and  Fuchs* 
theory  (Section  3.1.6).  Zaid ‘ s  theory  (Section  3.1.7) 
was  not  included  since  reliable  solutions  could  not  be 
obtained. 


Grimminger 1 s  theory  contains  two  somewhat  arbitrary 
constants;  the  transition  velocity  at  which  the  penetration 
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Piechanism  changes  .trow  a  fluid  drag  mechanism  to  one  of 
armor  penetration,,  and  the  constant  in  Che  armor  penetration 
relatione  For  the  present  comparison,  the  values  tor  these 
two  constants  given  by  Grimminger  were  taken.  The  transi¬ 
tion  velocity  was  taken  as  live  times  the  plastic  wave  velo¬ 
city,  the  latter  being  given  as  1G00  ft/sec  for  all  mate¬ 
rials,,  while  the  armor  penetration  constant  is  defined  in 
Section  3,1.5,  Using  these  values,,  it  may  be  seen  that 
Grimminger * s  theory  generally  underestimates  the  penetra¬ 
tion  below  5000  ft/sec,  and  overestimates  the  penetration 
above  this  velocity.  The  discrepancy  is  exceptionally 
large  in  the  case  of  lead  targets,  while  agreement  is  best 
for  copper  targets, 

A  change  in  the  constants  would  have  the  effect  of 
altering  the  level  of  the  portion  of  the  curve  above  the 
transition  velocity.  The  most,  obvious  improvement  would  be 
to  make  the  transition  velocity  a  multiple  of  the  elastic 
wave  velocity.  Taken  relative  to  the  results  in  copper., 
this  would  lower  the  predicted  penetration  in  lead  consid¬ 
erably  but  slightly  raise  the  predicted  penetration  in 
aluminum  and  steel,,  thus  representing  an  improvement  in 
agreement,  in  the  case  of  lead,  but  a  deterioration  in  the 
case  oi:  aluminum  and  steel, 

Bohn  and  Fuchs*  theory  contains  a  shape  factor, 
which  wan  assumed  to  lie  between  1/5  and  1,  For  the  value 
of  2/3  given  tor  spherical  projectiles..  It  may  be  seen  that 
the  theoretical  predictions  are  in  rather  surprisingly  good 
agreement  with  the  data  at  low  velocities.  It  may  be  seen 
from  Fig.  1.2  to  125  that  tor  most  materials  the  agreement 
is  very  good  tor  low  velocities,  the  theory  overestimating 
prui'frnt  iox).  at  higher  velocities.  For  tungsten  carbide 


the  theory  underestimates 


projectiles  at  low  velocities, 
penetration.  After  projectile  break-up  however,  the 
theory  again  overestimates  penetration., 

An  increase  in  shape  factor  .t  ,  would  reduce  the 
predicted  penetration.  At  higher  velocities  the  penetration 
is  roughly  proportion  to  1/f  ,  thus  a  value  of  f  of  about 
4/5  instead  of  2/3  would  lead  to  considerably  improved  agree¬ 
ment  over  the  whole  velocity  range,  except  tor  the  tungsten 
carbide  projectiles. 

5.2.2  Hydrodynamic  Theories 

Bjork* s  theory  (Section  3.2.1),  Oplk's  theory 
(Section  3.2.2)  and  shaped  charge  jet  theory  (Section  3.2.3) 
are  considered. 

Numerical  solution  of  the  equations  of  fluid 
dynamics  have  been  given  by  Bjork  for  only  two  metallic 
target  materials  at  the  time  of  this  study,  i.e.  aluminum 
impacting  aluminum,  and  iron  impacting  iron.  The  theory 
contains  no  adjustable  constants,  all  material  properties 
apnearing  in  the  theory  having  been  evaluated  in  indepen¬ 
dent  measurements.  Due  to  the  nature  of  the  assumption 
that  material  strength  is  negligible,  the  theory  is  only 
expected  to  be  applicable  at  very  high  velocities. 

The  theoretical  curve  given  by  Bjork  for  aluminum 
impacting  aluminum  has  been  drawn  on  Fig.  1.2  through  1.4 
for  2024-T3,  24ST  and  HOOF  targets  respectively  and  for 
iron  impacting  iron  on  Fig.  1.12  through  1,14  representing 
1020,  1030,  and  4140  targets  respectively.  Bjork  arbitrarily 

defined  the  equivalent  diameter  of  his  cylindrical  projectile 

2  1/^ 

as  d  =>  (d  Ll  ^  .  Definition,  of  the  equivalent  diameter 


on  the  basis  of  equivalent  mass  Eq „  (5„2)  leads  to  a 
reduction  in  predicted  penetration  of  16  per  cent 5  and 
leads  to  the  lower  curves  shown  in  Fig.  1„2  through  1.4 
and  1.12  through  1 „ 14 „ 

I  may  be  seen  that  the  highest  velocity  experi¬ 
mental.  points  for  the  HOOF  aluminum  and  1030  steel  targets 
arc?,  clustered  on  Bjork**  original  curve,  while  experimental 
points  for  the  other  alloys  lie  well  below  the  theoretical 
curve,,  However  the  experimental  points  show  a  steeper 
gradient  than  slope  predicted  by  the  theory,  so 

that  the  experimental  points  disagree  with  theoretical  pre¬ 
dictions  at  lower  velocities „  No  conclusions  regarding  pos¬ 
sible  agreement  at  higher  velocities  may  be  drawn  on  the 
basis  of  the  present  evidence.. 


Opik  s  theory  (Section  3 >2. >2)  contains  an 

2 

adjustable  constant  K.  =  1/2  vq  where  vq  represents 
the  velocity  below  which  no  penetration  will  occur.  In 
calculating  the  theoretical  curves,  vQ  was  found  by  extra¬ 
polating  the  data  to  zero  penetration,.  It  is  clear  from  the 
curves,  (Fig,,  1.2,  1.1,4,  1.15,  1.18,  1.22,  1.24)  that  the 
agreement  between  Opik ' s  theory  and  experiment  is  rather  poor; 
in  particular,  Opik ’ s  theory  predicts  that  the  rate  of  in¬ 
crease  of  penetration  with  increasing  velocity  becomes  very 
low  at  velocities  above  10,000  ft/sec. ,  A  change  in  the 
value  of  v  alters  the  predicted  penetration  but  does  not 
appreciably  alter  the  slope  of  the  curve  at  higher  velocities, 


O  tr  nyniT 


seen  from  Fig  1.18  where  tlicujTctital  curves  for 


v  =  200  ft/scc  and  300  ft/sec.  are  shown,, 


First  order  shaped  charge  jet  theory  predicts  that 
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penetration  will  be  independent  of  velocity  (Section  3.2,3) 


p  _  /  ' 

l  '  /  ;V 

which  may  be  expected  to  apply  only  at  very  high  velocities. 
Including  a  correction  for  target  strength  reduces  the  pre¬ 
dicted  penetration  at  low  velocities.,  the  theoretical  curve 
becoming  asymptotic  to  the  maximum  value  predicted  by  first 
order  theory.  Horizontal  lines  representing  first  order 
theory  have  been  drawn  on  each  of  Fig.  1.2  through  1.25. 

It  may  be  seen  that  there  is  a  consistent  trend  for  the 
penetration  to  be  considerably  greater  than  the  predicted 
value  at  higher  velocities.  Inclusion  of  the  correction 
for  target  strength  will  lead  to  lower  predicted  penetra¬ 
tions.  No  attempt  has  been  made  to  calculate  these  curves. 

The  lack  of  agreement  is  not  surprising,  since  the 
major  contribution  to  the  penetration  at  higher  velocities 
may  be  expected  to  be  due  to  the  inertial  expansion  of  the 
crater.  The  correction  for  inertial  expansion  was  not 
included  since  its  value  depends  strongly  on  the  velocity 
vq  for  which  no  penetration  occurs.,  and  which  must  be 
found  empirically.  This  could  not  be  done  accurately  with 
the  present  data. 

5 . 2  ■■  3  Thermal  Penetration 

The  applicable  thermal  penetration  theories 
are  those  of  Whipple,,  (Section  3.3.1)  Langton  (Section  3.3.2) 
and  Grow  (Section  3.3.3).  verves  representing  the  three 
theories  are  drawn  on  Fig.  1.2  through  1.25  for  comparison 
with  the  data. 
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Whipple's  and  Lang  ton "  a  theories  differ  mainly  in 
the  proportionality  constant  between  the  penetration  and  the 
2/3  power  of  the  velocity,  due  to  the  different  shape  assumed 
for  the  crater..  Langton  in  addition  includes  a  term  to  account 
for  the  melting  of  the  projectile,,  which  becomes  relatively 
less  important  as  the  velocity  increases,  Whipple's  theory 
overestimates  penetration  in  almost  ail  cases  except  for 
penetration  of  tungsten  carbide  projectiles  in  the  low  velo¬ 
city  region.  On  the  other  hand,,  Langton1 s  theory  consider¬ 
ably  underestimates  the  penetration  in  all  cases.  As  may  be 
expected  from  the  dependence  of  penetration  on  the  first  power 
in  velocity,  Grow's  theory  grossly  overestimates  penetration 
at  all  but  the  lower  velocities  in  all  cases  for  which  it  was 
computed , 

None  of  the  thermal  penetration  equations  take 
into  account  the  strength  of  the  target,  since  they  were 
intended  to  apply  only  to  extremely  high  impact  velocities,, 
and  they  can  therefore  not  be  made  to  fit  data  in  targets  of 
the  same  material  butt  different  hardness  within  the  experi¬ 
mental  velocity  range. 


Empirical  Power  Law 


In  this  section  it  is  assumed  that  a  penetration 
law  of  the  form 

p  k  v" 

is  valid.  Sufficient  experimental  data  exists  so  that:  the 
dependence  of  penetration  on  each  experimental  parameter  may 
te  investigated  separately,  and  a  non-dimensional  penetration 
law  is  devised  Crater  diameters  and  volumes  are  then  consi¬ 
dered 


5.3.1  Consistency  of  the  Data 


Several  projectile  material- target  material 
combinations  have  been  used  at  different  laboratories  in 
the  same  velocity  range*  thus  allowing  a  direct  comparison 
of  the  data  from  one  laboratory  with  that  from  another. 

Except  for  some  early  work.,  the  experimental  data  agree 
within  the  experimental  scatter*  as  may  be  seen  in  Fig.  1.2 
through  1.25 o  Apparent  disagreements*  as  in  the  case  of 
aluminum  projectiles  impacting  aluminum  targets  (Fig.  1.2* 
1.3,  1.4)  may  be  traced  to  differences  in  strength  properties 
of  the  targets  used  at  the  various  laboratories.  In  subse¬ 
quent  statistical  analyses*  all  of  the  data  was  weighted 
equally. 


5.3.2  Dependence  on  Projectile  Dimension 

For  large  number  of  projectile  material- target 
material  combinations*  firings  have  been  made  with  a  variety 
of  projectile  sizes*  and  with  both  shperlcal  and  cylindrical 
projectiles.  In  plotting  Fig.  1.2  through  1,25*  the  penetra¬ 
tion  has  been  normalised  to  the  projectile  diameter.  In  the 
case  of  cylindrical  projectiles*  the  diameter  of  an  equiva¬ 
lent  sphere  of  mass  equal  to  that  of  the  cylinder  was  used* 
i  o  e . 


However*  it  should  be  pointed  out  that  nearly  all  of  the 
projectiles  considered  had  length  to  diameter  ratios  (I/d) 
near  unity*  and  the  equivalent  sphere  diameter  did  not 
differ  greatly  from  the  cylinder  diameter.  Almost  equally 
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good  correlation  could  be  obtained  by  using  the  cylinder 
diameter  or  length  directly,,  For  L/d  far  from  unity.,  this 
scaling  is  not  likely  to  hold.  (See  Section  6,3)  Insuffi¬ 
cient  data  exists  at  present  to  de'fine  the  limits  of  L/d 
over  which  the  correlation  on  the  basis  of  equivalent  sphere 
diameter  can  be  expected  to  hold. 

Except  for  a  consistent  tendency  in  the  data 
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reported  by  Kinard  and  Collins  for  the  penetration  to  be 
slightly  greater  for  1/2  inch  diameter  projectiles  than  0..22 
inch  diameter  projectiles  no  sise  effect  could  be  detected, 
Nearly  all  of  the  data  was  in  the  range  of  projectile  diam¬ 
eter  of  1/16  to  1/2  inch., 

5o3o.3  Dependence  on  Projectile  Velocity 

If  it  is  assumed  that  the  penetration  depends  on 
the  velocity  through  a  relation  of  the  form  p  =  kVn  this 
type  of  dependence  can  be  most  conveniently  investigated  by 
means  of  log- log  plots  of  penetration  versus  velocity  ..  Typi¬ 
cal  curves  of  this  type  for  ductile  projectiles  impacting 
aluminum,  copper  and  lead  targets  are  shown  in  Fi.g„  5, .12, 

3.13  and  5..  14  and  for  brittle  projectiles  in  Fig..  3, ,13.,  5,16 
axrd  5. ,17.  Although  there  is  considerable  scatter,,  the  curves 
all  show  a  tendency  to  vary  from  a  velocity  exponent  of  4/3 
at  low  velocities  to  about  2/3  at  high  velocities,  (except 
for  lead.  Fig-  3. 14  which  does  not  extend  to  low  enough  velo¬ 
cities  to  show  a  4/3  power  dependence)  For  the  high  strength 
projectiles,,  the  4/3  slope  is  maintained  to  velocities  at  which 
the  proj ectile  shatters,  this  effect  being  most  pronounced  in 
the  softest  material.  Lead,. 
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5.3.4  Dependence;  on  Projectile  Properties 


Relatively  little  information  is  available  on  the 
effect  of  projectile  strength  on  cratering.  Maiden  et  al^" 
shot  series  of  steel  projectiles  which  had  received  differ¬ 
ent  heat  treatments  into  1030  steel  targets  in  a  velocity 
range  from  4,000  to  13,000  ft/sec.  Results  (Fig.  5.18) 
indicate  that  no  significant  differences  in  penetration 
are  observed  for  projectile  Brinell  Hardness  Numbers  of 
210,  290  and  580.  It  is  clear  that  these  data  lie  in  the 

high  velocity  region,  since  they  conform  to  the  2/3  velo~ 
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city  exponent.  On  the  other  hand  Abbot  found  signifi¬ 
cant  differences  in  penetration  in  steel  targets  by  two 
types  of  steel  projectiles,  having  Brinell  Hardness  Numbers 
of  140  and  705  respectively,  up  to  7,000  ft/sec,  i.e. 
well  into  the  high  velocity  range,  although  the  differ¬ 
ences  were  small  above  5,000  ft/sec,  (Fig,  5.19.) 

Q  £ 

Grow  et  al  shot  two  types  of  steel  projectiles  into 
lead  targets,  (Fig  5.20.)  up  to  a  velocity  of  about  5,000 
ft/sec.  Significant  differences  in  penetration  were  ob¬ 
served  up  to  a  velocity  of  about  4,000  ft/sec,  which  from 
Fig,  5.17  is  seen  to  correspond  to  the  upper  limit  of  the 
transition  region, 

It  may  be  surmised  that  the  projectile  strength 
does  not  affect  penetration  above  the  transition  region, 
particularly  for  ductile  projectiles  It  Is  clear  that 
the  same  will  be  true  of  the  undeformed  projectile  region 
tor  high  strength  projectiles,  where  no  deformation  of  the 
projectile  occurs. 

If  this  is  true,  then  the  only  factor  to  account 
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for  differences  in  penetration  in  a  given  target  material 
by  different  projectile  materials  is  the  projectile  density 
By  analysing  proportionality  constants  k  in  p  *  kV 
obtained  from  plots  of  the  type  of  Fig.  5.12  through  5.17, 
k  is  found  to  be  very  nearly  proportional  to  the  2/3  power 
or  the  projectile  density  and  independent  of  projectile 
strength  in  me,  high  velocity  region,  for  a  given  target 
material,. 


5.3. 5  Dependence  on  Target  Properties 

While  penetration  may  be  expected  to  depend  on 
target  density,  there  is  strong  evidence  no  suggest  that 
penetration  also  depends  on  target  strength  over  the  entire 
vs rue xiv  range  covered  by  the  present  data.  Summers  has 
reported  firings  of  copper  projectiles  into  two  sets  of 
copper  targets  heat  treated  to  Brinell  Hardness  of  65  and  35 
respectively  (Fig.  5.21)  Significant  differences  in  pene¬ 
tration  were  observed  over  the  whole  velocity  range  which 
was  explored,  extending  from  a  tew  hundred  ft;/ sec.,  to  10,000 
ft/sec  Figure  3  22  shows  penetration  results  for  steel  pro¬ 
jectiles  impacting  four  different,  types  of  steel  target  ranging 
in  Brinell  Hardness  from  ill  to  302..  Differences  in  pene¬ 
tration  arc:  again  evident  up  to  at  least  8,000  ft/sec,  the 
highest  velocity  at  which  a  comparison  can  be  made. 

A  comparison  of  cratering  into  hard  (2024T-3)  alu¬ 
minum  and  soft  (HOOF)  aluminum  targets  by  aluminum  projectiles 
can  also  bo  made  (Fig  5  23)  Penetration  in  the  2024T-3 
aluminum  targets,  having  a  Brinell  Hardness  ot  about  120  is 
observed  to  be  about  60%  of  that  in  the  softer  HOOF  alloy 
up  to  the  highest  experimental  velocity  oi  16,000  ft/sec. 


The  hardness  of  the  HOOF  material  was  unfortunately  not 
measured,  but,  according  to  handbook  values,  lies  in  the 
range  23  to  44. 

It  is  therefore  clear  that  target  strength  will 
have  to  be  taken  into  account  over  the  whole  of  the  experi¬ 
mental  velocity  range.  It  has  already  been  shown  that 
penetration  is  proportional  to  the  2/3  power  of  projectile 
density  and  velocity.  The  following  procedure  was  adopted 
to  verify  this  dependence  on  projectile  density  and  velo¬ 
city,  to  obtain  reliable  values  of  the  proportionality 
constant  for  each  target  material,  and  to  determine  the 
lower  limit  of  velocity  in  which  the  2/3  power  law  is 
valid  for  each  projectile  material- -target  material  com¬ 
bination. 
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2/3 

Large  scaLe  plots  of  ij  -  p/d  (p^V)  '  versus  velocity 
wore  prepared  for  each  target  material.  If  the  penetration 
did  in  tact  follow  the  law 


F 
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the  resultant  plot  should  be  a  horizontal  straight  line 
as  shown  in  sketch  (I)  above-  If  the  penetration  were  pro¬ 
portional  to  some  other  exponent  of  velocity,  the  resultant 
plot  would  be  a  curve  as  shown  in  sketch  (ii)  above-  At 
high  velocities,  the  curve  should  reach  the  2/3  power  law, 
and  thus  become  a  horizontal  straight  line..  The  velocity 
above  which  the  curve,  becomes  horizontal,,  then  corresponds 
to  the  lower  limit  of  the  high  velocity  region..  It  is  to 
be  expected  that  different  projectile  materials  will  produce 
different  curves  at.  low  velocities  as  shown.,  If  on  the 
other  hand  the  penetration  were  proportional  to  some  other 
exponent  in  projectile  density.,  the  resultant  curves  for 
different  projectile  materials  should  tend  to  horizontal 
lines  which  are  separated.,  as  shown  in  sketch  (iii)  above. 
Curves  for  3  types  of  aluminum,  4  grades  of  copper,  4  types 
of  steel,  lead., magnesium,  silver,  zinc,  cadmium  and  tin 
targets  impacted  by  a  variety  of  projectile  materials  were 
prepared  Examples  for  2024T-3  aluminum, copper  of  BilN  b5, 
and  tor  lead  targets  are  shown  in  Fig  3.24,  5.23  and  5,26., 
By  tar  the  greatest  variety  of  projectile  materials  were 
used  with  these  three  target  materials.  It  may  be  seen  that 
there  are  relatively  very  few  experimental  points  in  the 
horizontal  region  of  the  curves,  but  that  the  curves  for  the 
various  projectile  materials  do  converge  on  the  horizontal 
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portion  within  the  experimental  scatter.  Only  very  rough 

values  of  the  limiting  velocity  could  be  obtained.  However, 

the  value  of  p/d  (PpV)*^  =  corresponding  to  the 

horizontal  portions  of  the  curves  could  be  obtained  fairly 

accurately  for  several  of  the  target  materials.  Least 

squares  fits  of  the  values  of  ,  together  with  standard 

mean  deviations,  were  obtained  using  experimental  points  in 

the  horizontal  portions  of  the  curves.  These  are  listed  in 
Table  IV. 

The  effect  of  target  strength  could  be  estimated 
separately,  since  k^  could  be  evaluated  for  several  types 
of  aluminum,  copper,  and  steel  targets,  which  varied  in 
Erinell  Hardness,  but  had  closely  similar  densities  and  sonic 
velocities.  Values  of  k^  found  from  least  squares  fit  of 
the  data  were  plotted  versus  the  Brinell  Hardness  Number  1^ 
on  log- log  paper.  (Fig.  5.2'/)  The  flags  represent  one 
standard  mean  deviation  in  and  the  uncertainty  ranges 

in  Hardness.  Straight  lines  of  slope  -1/3  could  be  fitted 
through  the  points  representing  the  different  types  of  each 
material.  Thus  the  target  strength  can  be  approximately 
accounted  for  by  including  a  term  i.e. 


d 


a*  vv> 

Ht  6 


the  values  of  k9  being  found  from  the  previous  values  of 
k]  by 


The  effect  of  density  may  now  be  investigated 
by  plotting  resultant  values  of  k^  versus  target  density 
on  log-log  paper  (Fig.  5.2b)  The  uncertainty  flags  shown 
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are  a  combination  of  the  mean  deviations  in  k-j  and  the 
uncertainties  in  „  A  straight  line  of  slope  “1/3  has 
been  approximately  fitted  to  the  points,  again  using  a  least 
squares  fit,  leading  to  a  non-dimensional  penetration  law 

f  -  <5-» 


While  most  materials  fell  within  the  limits  given  above, 
penetration  for  some  "anomalous"  material  combinations  may 
depart  from  the  above  expression  by  as  much  as  35  per  cent.. 

It  might  be  noted  that  a  slightly  higher  exponent 
in  pt  might  be  expected  to  fit  slightly  better „  However  it 
was  decided  t.o  retain  the  advantages  of  a  non-dimensional 
fit.  Small  changes  (+  10  "Z  )  in  exponents  of  the  non-dimen  ¬ 
sional  parameters  in  Eq„  (5.1)  did  not  significantly  alter 
the  mean  deviation. 

it  should,  be.  noted  that  the  empirical  expression 
Eq„  (<> I)  is  necessarily  a  crude  approximation,  since  it 
disregards  the  many  complex  phenomena  actually  occurring,. 


S„ 3., 6  Limits  of  Validity  of  the  Empirical  Power  Law 

The  empirical  law  deduced  In  the  previous  section 
is  clearly  limited  to  a  restricted  velocity  range  for  each 
combination  of  projectile  and  target  materials..  Nothing  can 
be  said  at  this  time  about  the  upper  velocity  limit,  on  the 
basis  of  the  experimental  data.  The  penetration  appears  to 
be  proportional  to  the  2/3  power  in  velocity  up  to  the  high 
est  velocity  data  points  included  here,  within  the  experimental 
scatter,.  This  point  will  be  discussed  later,.  However,  the 
lower  limit  of  validity  of  the  empirical  law  may  be  examined 
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in  some  detail* 


2/3 


A  number  of  interesting  trends  appear  from 
Fig.  5.24,  5=25  -and  5*26.,  in  which  the  ratio  p/d  f  (p^V) 
is  plotted  vs.  velocity  for  aluminum,  copper  and  lead  targets. 
Curves  have  been  roughly  faired  through  some  of  the  points 
for  the  sake  of  clarity.  It  is  seen  that  the  penetration  is 
generally  less  for  the  low  density  projectiles,  and  greater 
for  the  high  density  projectiles  at  low  velocities  than  pre¬ 
dicted  by  the  empirical  penetration  law  Eq„  (5.1),  the  differ¬ 
ences  diminishing  as  the  velocity  increases,  until  the  various 
curves  combine.  While  the  high  velocity  data  is  limited  al¬ 
most  exclusively  to  aluminum  projectiles,  the  data  repre¬ 
senting  other  projectiles  do  appear  to  become  asymptotic  to 
the  horizontal  lines  faired  through  the  high  velocity  points 
for  aluminum  projectiles. 


On  the  basis  of  the  qualitative  model  of  cratering 
presented  in  Section  5.1,  it  is  possible  to  suggest  the 
reasons  for  this  behaviour.  Reference  to  Fig.  5.1  shows  that 
the  initial  maximum  pressure  generated  at  the  interface  is 
less  for  the  lower  density  projectiles,  and  greater  for  the 
higher  density  projectiles. 


One  might  therefore  suppose  that  the  strength  of 
the  target  material  may  have  a  relatively  greater  effect  in 
the  case  of  the  lower  density  projectiles,  and  therefore 
restrict  the  flow  and  result  in  smaller  craters  in  a  given 
target  material  than  predicted  on  the  basis  of  the  empirical 
law  Eq.  (5*1).  Clearly  at  some  low  velocity,  the  initial 
interface  pressure  will  be  so  low  that  no  target  deformation 
will  occur.  This  velocity  will  then  be  higher  for  the  lower 
density  projectiles.  The  opposite  will  be  true  of  higher 
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density  projectiles <•  One  might  suppose  that  the  target  defor¬ 
mation  is  related  to  the  quantity  (P  ~  i<,e„  the  initial 

interface  stress  in  excess  of  the.  dynamic  yield  of  the  material o 

At  the  same,  time,  there  will  be  a  similar  effect 
on  the  projectile,.  The  degree  of  deformation  of  the^projee™ 
tile  tnay  also  be  roughly  related  to  the  quantity  (P  -  Y^)  - 
Thus  „  although  the  initial  interface  pressure,  is  low  for 
nylon  projectiles ,  nylon  also  has  very  low  strength,  and  a 
nylon  projectile  may  be  expected  to  suffer  severe  deforma¬ 
tion  with  consequent  increase  in  presented  area  and  decrease 
in  penetration  at  even  moderate  velocities »  On  the  other 
hand  tungsten  and  uranium  projectiles „  while  leading  to  much 
higher  interface  pressures  at  the  same  velocities  have  rather 
high  strength.,  Thus  deformation  of  the  projectile  will  be 
inhibited  and  penetration  enhanced,,  leading  to  "brittle  pro¬ 
jectile"'  type  of  behaviour „  (Section  5° 1) « 

On  the  basis  of  this  qualitative  model,  one  might 

expect  that,  the  empirical  penetration  law,  Eq„  ( 5 -  1 )  .<  will 

be.  valid  when  p”*  Y  s  Y..  „  While  the  limiting  velocities 

can  be  estimated  only  very  approximately  from  curves  such  as 

Fig,,  5„24  through  it  would  appear  that  a  reasonable 

~  '  k 

criterion  is  given  by  P  >  100  Y^  where  Y^_  is  the.  static 
yield  ,  for  material  combinations  not  exhibiting  "brittle 

projectile"  behaviour  Noting  that  the  Brine  11.  Hardness  is 
roughly  equal  to  j .  7-  from  Fig,,  ..  1  v  the  criterion  reduces  to 


For  high  strength  projectiles,  the  velocity  at. 
which  breakup  is  complete  may  be  higher  than  the  velocity 
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indicated  in  Eq„  (5.2).  This  is  particularly  true  in  lead 
and  other  low  strength  targets .  Unfortunately  the  hardness 
of  projectiles  used  in  experiments  has  been  reported  in  only 
very  rare  instances,  so  that  almost  no  information  exists 
to  guide  formulation  of  a  criterion  similar  to  Eq„  (5«2)» 

Only  uue  piece  of  information  is  available,  i.e.  stainless 
steel  projectiles  shot  into  lead  targets  at  Utah  University 
(Fig.  5.20).  In  this  one  case  it  also  appears  satisfactory 
to  use  the  criterion 

P *■  >  32  (5.3) 

While  it  seems  logical  to  extend  Eq.  (5.3)  to  other  materials, 
there  is  presently  no  other  evidence  to  support  such  an  assump¬ 
tion.  If  however,  the  assumption  is  made,  the  criterion  be¬ 
comes 

>  32  and  >30  (5.4) 

where  P  and  H  are  in  identical  units. 

An  attempt  has  been  made  to  apply  empirical  correc¬ 
tions  to  the  low  velocity  data  in  order  to  account  for  the 
strength  effects  by  the  use  of  correction  terms  of  the  type 

Je 

P 

—ac -  ,  where  k  is  a  constant  to  allow  for  strain  rate 

P  -  kH  (. 

effects.  This  is  very  difficult  because  the  target  and  pro¬ 
jectile  strength  effects  apply  simultaneously,  and  generally 
tend  to  act  in  the  same  direction.  Furthermore ,  the  uncer¬ 
tainty  in  the  data,  (generally  +  .05  in  p/d  or  greater)  is 
such  that  at  low  velocities,  and  hence  low  penetrations,  the 
percentage  error  is  large.  It  is  therefore  impossible  to 
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evaluate  the  effects  of  target  and  projectile  strength  separately 
with  the  present,  data* 

5 .>3<.  7  Crater  Diameters 


Diameters  of  craters  in  the  plane  of  the  original 
tetgct  i-izrfate  have  not  been  reported  in  all  cases  c  Figure 
5.2.9  through  5.32  show  typical  trends  in  the  data  for  aluminum, 
copper ,  steel  and  lead  targets .  Shown  are  the  ratios  of 
penetration  depth  to  crater  diameter  as  a  function  of  velocity. 
This  ratio  is  in  effect  a  shape  factor,  which  should  equal  0„ 5 
for  hemispherical  craters,  higher  values  indicating  craters 
deeper  than  h<-  mi  spherical. 

U.  may  be  seen  that  projectiles  such  as  copper  and 
lead  show  a  tendency  for  the  ratio  p/D^  to  increase  with 
velocity,  and  become  asymptotic  to  the  value  0  „  c Reference 
to  Fig,,  5.,  12,  5nl3,  5  <>14  indicates  that,  the  penetration  be¬ 
comes  a  function  of  at  about  the  same  velocity  where 

the  ratio  p/D,  reaches  O.5,  These  therefore,  represent 
examples  of  ductile  projectile  behaviour,  described  by  many 
authors . 

Of  interest  is  the  behaviour  of  aluminum,,  magnesium 

and  nylon  projectiles.  In  each  case,  the  ratios  p/Df  appear 

to  become  asymptotic  to  a  value  of  0.4,  and  possibly  less  in 

the  case  of  nylon,,  in  the.  velocity  region  where  the  penetra- 

2/1 

tion  is  proportional  to  V 

It  may  also  be  seen  that  projectiles  such  as  hardened 
steel  and  tungsten  carbide  show  a  tendency  for  p/D  f.o 
increase  with  velocity  to  a  maximum  much  greater  than  0  ,  p , 
arsd  then  to  decrease  and  in  some  cases  show  a  tendency  to 
became  asymptotic  to  0„c>..  Ihe  velocities  at  the  maxima  of 


the  p/D  curves  may  be  correlated  with  the  velocities 
at  the  maxima  of  the  p/d  curves  Fig.  5.15,  5. 16  and  5.17. 
These  therefore  represent  examples  of  brittle  projectile 
behaviour.  However,  in  several  eases,  the  ratio  p/°c 
has  approached  a  value  considerably  greater  than  0.5  in  the 
velocity  region  corresponding  to  a  penetration  dependence 

rt 

on  V4/:5  „ 

This  type  of  behaviour  is  exactly  what  one  might 
expect  on  the  basis  of  the  qualitative  model  discussed  in  the 
previous  section.  (Sect.  5.3.6)  The  question  as  to  the  vel¬ 
ocities  at  which  the  craters  become  hemispherical,  or  whether 
the  craters  will  become  hemispherical  at  all,  for  the  very 
light  and  very  heavy  projectiles  cannot  be  answered  by  the 
present  data.  However  it  is  clear  that  the  craters  are  not 
necessarily  hemispherical  in  the  high  velocity  region  covered 
by  the  empirical  penetration  law  Eq.  (5.1). 

5.3.8  Crater  Volumes 

Several  experimenters  report  measured  crater  volumes „ 
The  volumes  have  generally  been  measured  by  machining  the 
crater  lip  flush  with  the  original  target  surface  and  filling 
the  crater  with  a  liquid,  wax,  or  powder,  the  volume  of  which 
is  subsequently  measured. 

Of  considerable  interest  is  the  ratio  E/Vc  where 
E  is  the  initial  kinetic  energy  of  the  projectile,  and  V 

c 

the  crater  volume.  The  reciprocal  quantity  V  / E  has  some¬ 
times  been  termed  the  cratering  efficiency,  despite  the  fact 
that  the  quantity  is  dimensional..  Typical  plots  of  E/Vc 
vs.  projectile  velocity  are  presented  in  Fig„  5.33  through 
5° 36"  For  several  material  combinations  a  trend  has  been 
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noted  by  many  authors  for  12/V^  to  become  constant  as  the 
velocity  increases.  Several  empirical  penetration  laws  are 
based  on  this  assumption*  For  many  material  combinations s 
data  exists  only  at  low  velocities s  however  several  material 
combinations  in  Fig„  5»33  through  5.36  do  show  this  trend  at 
higher  velocities?  it  may  be  noted  that  the  E/V _  curves 
become,  horizontal  at  roughly  the  same  velocities  that  the  pene 
tration  becomes  proportional  to  (Fig0  5 „24  through  5,26 


Asymptotes  of  E/V  curves  were  estimated  from 
plots  of  the  type  of  Fig»  5„  33  through  5,36-  There  were 
generally  insufficient  points  in  the  horizontal  portions  of 
the  curves  to  warrant  statistical  fits? 

A  tendency  for  the  asymptotic  values  of  E/V  to 

depend  on  the  density  ratio  p  /pt  was  notec*  by  Summers  J  „ 

Such  a  trend  is  also  obvious  from  Fig,,  5 , 33  through  5„  36., 

Values  of  E/V  were  obtainable  for  several  different  pro- 
C 


jectile  materials  impacting  lead,  copper „  steel  and  aluminum 
targets „  Consequently,  asymptotic  values  of  E/V  were 
plotted  versus  density  ratio  Pf/Pp  00  log  log  paper  (Fig, 
5.37)»  The  points  for  each  target,  material  may  be  approxi¬ 
mately  fitted  by  straight  lines  of  slope  2/3*  indicating  0 
dependence  of  E/V.,  on  (p^/ pp)^^  » 

In  order  to  investigate  the  dependence  on  target 


properties  j,  values  cf 


.£  (fty 
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were  plotted  vs*  Brine  11  Hardness*  (Fie  A  « <»ht 

line  may  be  approximately  fitted  to  the  data  leading  to  a 


96 


final  correlation 


/  3  •  /  Jt 

V 


0 


Vi 


H. 


(5*5) 


where  the  constant  is  dimensionless,  and  independent  of 
projectile  and  target  material  properties . 

It  is  interesting  to  note  that  if  the  craters 
are  assumed  to  be  hemispherical  in  the  high  velocity 
region,  a  law  of  the  form  of  Eq.  (5.5)  may  be  derived 
directly  from  Eq„  (5,1)  by  noting  that 

E  -  r  ¥  J*fp  V' 


K 


A. 


-tir  .  3 
T  P 


resulting  in  an  expression 


t 

14 


(5=6) 


This  differs  slightly  from  Eq,  (5,5)  in  both  the  constant 
and  exponent  in  density  ratio.  However,  in  the  previous 
section  (Sect,  5„3 • 7)  it  was  shown  that  craters  formed  by 
low  density  projectiles  were  shallower  than  hemispheres, 
while  craters  formed  by  high  density  projectiles  were  deeper 
than  hemispheres.  Thus  Eq„  (5,6)  will  predict  smaller  cra¬ 
ters  than  actually  formed  for  aluminum  projectiles,  and  will 
overestimate  craters  actually  formed  by  tungsten  and  tung¬ 
sten  carbide  projectiles.  The  change  in  exponent  of  the 
density  ratio  just  compensates  for  this  effect,  as  may  be 
seen  from  Fig,  5,37.. 
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5a 3. 9  Comparison  with  Previous  Empirical  Expressions 

In  this  section  it  is  intended  to  compare  the  em¬ 
pirical  expression  Eq.  (5»1)  with  empirical  expressions 
obtained  previously  by  other  authors  and  listed  in  Table  III 
Many  of  these  expressions  were  devised  on  the.  basis  of  por¬ 
tions  of.  the  experimental  data  used  to  derive  Eq ,  (3,1),  and 
should  therefore  be  comparable  over  the  limited  ranges  in 
experimental  parameters  considered  in  their  derivation.. 

Several  empirical  expressions  (Table  III)  are  based 
on  a  dependence  of  penetration  on  velocity  to  a  power  greater 
than  2/3.  It  is  clear  from  previous  discussion  that,  these 
expressions  were  derived  from  data  in  a  lower  velocity  range 
than  that  covered  by  Eq.  (5«1)0  In  Section  5.3.6  the  compli¬ 
cations  arising  from  projectile  and  target,  material  strength 
have  been  discussed.  in  general ,  both  very  low  and  high 
density  projectile  data  have  been  Ignored  in  the  derivation 
of  these  expressions,  thus  allowing  a  reasonable  correlation 
on  the  basis  of  elementary  material  properties. 

A  number  of  empirical  expressions  are  based  on  the 

2g 

form  proposed  by  Charters  and  Locke 


Attempts  were  made  to  reduce  the  data  on  this  basis. 
However  ,  since  the  sonic,  velocity  cf.  does  not  vary  appreciably 
with  hardness  in  a  given  target  material  .,  expressions,  of.  this 
type  are  inadequate  to  correlate  data  from  firings  into  tar 
gets  of  different:  hardness,  but  identical  material..  The 
assumption  made  in  the  derivation  of  Charters  and  Locke's 
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expression;  vis.  that  the  effect  of  hardness  is  negligible 
in  the  high  velocity  range,  is  contradicted  by  the  data  up 

to  the  highest  velocities  considered  here,  and  significantly 

«  2/? 
in  the  region  where  the  penetration  is  proportional  to  V  ' -  . 

Speculation  as  to  the  role  of  strength  at  much  higher  velo¬ 
cities  will  be  deferred  to  a  later  section.  Equation  (5.1) 
fits  the  data  in  half -hard  copper  and  lead  targets  about  as 
well  as  Eq.  (5.7),  used  by  Summers^  to  correlate  penetra¬ 
tion  in  these  materials .  It  is  interesting  to  note  that 
P-  v  Hugoniots  (Appendix  D)  are  almost  identical  for 
copper,  iron  and  lead,  so  that  behaviour  associated  with 
low  and  high  density  projectiles  is  absent  when  these  mater¬ 
ials  only  are  involved.  Furthermore,  it  might  be  pointed 
out  that  data  for  annealed  copper,  and  copper  of  greater 
Brinell  Hardness  cannot  be  correlated  by  Eq„  (5.7),  so  that 
the  choice  of  half-hard  copper  as  the  basis  of  Eq„  (5.7) 
must  be  considered  fortuitous. 


Expressions  derived  on  the  assumption  that  cratering 
efficiency  is  a  constant  proportional  to  target  hardness  and 
that  the  craters  are  hemispheres  in  the.  high  velocity  region, 
are  equivalent  to  the  form  Eq„  (4. 33) 


/t  /  v  H  t  J 


This  expression  differs  from  Eq„  (5„1^  only  in  the  exponent 
in  density  ratio  and  in  the  value  of  the  constant,  Attempts 
to  fit  the  high  velocity  data  used  in  the  derivation  of 
Eq.  (5.I)  to  Eq.  (5,8)  by  the  method  of  least  squares  led 
to  a  considerably  increased  mean  deviation. 

It  is  interesting  here  to  compare  the  expression 
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for  E/V  derived  by  Feldman*' 


i: 

K 


jl-6S  hi. 


with  the  egression  derived  in  Section  5-3o8 


(5. 9) 


(5-5) 


These  also  differ  in  the  value  of  the  constant 
and  exponent  in  density  ratio „  Attempts  to  fit  the  data  used 
in  the  derivation  of  Eq0  (5.5)  by  Eq.  (5..  9)  again  led  to  a 
considerably  increased  mean  deviation*.  Equation  (5-8)  was 
derived  from  Eq„  (5,9)  by  assuming  that  the  crater  is  hemi¬ 
spherical  <>  This  assumption  has  been  shown  to  be,  unjustifi¬ 
able  in  the  velocity  range  covered  by  the  experimental  data 
in  Section  5„3<,7  and  5,3„8«  Thus  it  is  not.  surprising  that 
Eq„  (5-8)  and  (5-1.)  do  not  correspond- 

Tjork  has  recently  suggested  that  the  penetration 
in  a  given  target  material  by  projectiles  of  identical  mass 
and  velocity,  but  different  materials,  is  proportional  to  the 
initial  interface  velocity  induced  by  the  impact-  For  pro¬ 
jectiles  of  identical  mass  and  velocity  in  3  given  target, 

1/‘J 

Eq„  (5,.  I.)  reduces  to  p  =  kp  ~  „  In  Fig,,  5-39  values  of 

p  1  j  > 

the  interface  velocity  V  have  been  plotted  against  p  - 

for  several  common  target  materials  at  a  projectile  velocity 

of  5  km/ sec.  The  interface  velocity  is  very  nearly  pro- 

portions  I  to  p  Since  the  V  vs.,  V  curves.  Fig,,  5,1, 

are  very  nearly  straight  lines,  a  similar  degree  of  scatter 

may  be  expected  at  other  velocities.  (See  Appendix  I>)  Thus, 

Bjork  a  suggestion  is  essentially  equivalent  to  Eq„  (5-1) 


To  be  pub 11  shed 
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within  the  experimental  accuracy.  Experimental  scatter 

i  /q 

does  not  permit  a  choice  between  V  and  p  ,D  for  cor- 

P 

relation  purposes,  and  the  density  dependence  has  been  re¬ 
tained  for  the  sake  of  simplicity,, 


5.4  Empirical  Logarithmic  Law 


In  this  section,  the  fit  of  the  experimental 
data  to  a  penetration  law  of  the  form 


(5° 10) 


is  examined,,  Based  on  the  discussion  of  Section  4.4,  it 
may  be  seen  that  this  expression  (Eq.  4.58)  and  that  of 
Bohn  and  Fuchs  (Eq.  4.60}  may  be  made  to  almost  coincide 
over  the  experimental  velocity  range  by  suitable  choice  of 
constants „  in  which  case  serious  disagreement  between  the 
two  expressions  would  only  arise  at  very  high  velocities. 

It  Is  therefore  to  be  expected  that  it  will  not  be  possible 
to  choose  between  these  two  expressions  on  the  basis  of  the 
present  experimental  data  alone.  Equation  (5. 10)  was  there¬ 
fore  chosen  as  the  basis  of  correlation  for  its  relative  sim 
plicity. 


b. 4.1  Data  Fits 

Individual  least  squares  fits  were  obtained  to 
the  experimental  data  for  each  projectile-target  material 
combination,  fitting  kT  and  in  Eq„  (15.10)  by  using 

a  standard  iteration  method  on  a  high  speed  digital  com¬ 
puter.  Material  combinations  exhibiting  brittle  projec¬ 
tile  behaviour  were  excluded,  since  in  these  cases  the  pene¬ 
tration  is  not  a  monotonically  increasing  function  of 
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velocity,  and  therefore  could  not  be  fitted  by  an  expression 
of  the  type  Eq,  (5»1Q),,  Results  together  with  standard  mean 
deviations  are  presented  in  Table  V. 

It  may  be  seen  that  the  individual  fits  are  very 
good,  mean  deviations  in  many  cases  being  no  larger  than  the 
scatter  in  the  data„  However  the  constants  and 

vary  for  the  different  projectile  target  material  combina¬ 
tions  o  Attempts  were  made  to  obtain  correlations  of  k^ 
and  k.£  with  elementary  material  properties <>  No  correlating 
could  be  found  with  H  or  c£  .  and  the  best  correlations 
obtained  are  shown  in  Fig„  5o^0  and  5«JH  respectively  on  the 
basis  of  the  density  ratio ,  leading  to  an  approximate  expres¬ 
sion  of  the  form 


£ 

A 


fa  *  i) 


(5.11) 


It  is  evident  that  the  correlation  is  rather  poor„ 
This  is  hardly  surprising „  The  effects  of  projectile  and 
target  strength  in  the  transition  region,  discussed  in 
Sections  5°  3-6  and  5., 3 7,  strong3.y  influence  k,  and  k0  , 
since  the.  fits  were,  made  using  data  from  the  entire  experi¬ 
mental  velocity  range.  It  could  therefore  hardly  be  expected 
that  a  simple  correlation  with  powers  of  elementary  material 
properties  would  be  obtained »  Before  embarking  upon  an  attempt 
to  correlate  and  with  more  complex  functions  of 

material  properties,  it  would  be  desirable  to  predict  reason 
abLe  forms  of  such  functions  on  theoretical  grounds,, 

When  material  combinations  exhibiting  extreme 
target  and  projectile  strength  effects,  such  as  magnesium  and 
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aluminum  with  copper,  steal  and  lead,  are  omitted,  a 
somewhat  improved  correlation  for  k-j  is  obtained ; 


However,  the  correlation  for  k9  is  not  improved. 

S.4.2  Theoretical  Considerations 

Iti  deriving  Eq.  (5. 10),  Section  4.4. 3,  the 

9 

instantaneous  interface  pressure  was  taken  to  be  1/2  pfcV  , 
assuming  the  target  to  be  an  incompressible  fluid,  and  the 
projeetlle  to  be  rigid  and  incompressible.  It  might  be 
expected  that  a  better  approximation  might  be  obtained  by 
assuming  that  the  interface  pressure  P  is  related  to  the 
interface  velocity  v  by  the  shock  relation  (Appendix  D) 

P  “  ft  /  (  v)  (5,13) 

This  assumes  that  the  conditions  between  the  shock  and 
interface  are  uniform,  which  would  appear  to  be  a  reasonable 
approximation,  except  in  the  later  stages  of  the  motion,  from 
pressure  contours  obtained  by  Bjork,,  Fig.  5.9  and  5. 10. 

Thus,  assuming  the  retardation  force  on  the  pro¬ 
jectile  to  be  of  the  form 

F  -  k{  (p  +  +  cl  (5.  is-) 

integration  of  the  equation  of  motion 

-  F  c\t  *  1/  c\v  (5° 15) 
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between  the  limits  x  =  0  ,  P  =  P  ,  v  =  V'  ;  s  *  p  „  P  a  0 
v  =  0  leads  to 


m 

(k,St  Vft  J 


P* 

I  f  _C_ 
K  H+ 


-  R 


(5° 16) 


where  R  depends  on  V  *"  but  may  be  shown  to  become  a 
relatively  small  constant  at  high  velocities » 


It  is  interesting  to  compare  Eq„  (5„16)  with 
Eq0  (5.11).  Using  the  approximate  expression  for  interface 
velocity,  obtained  in  Appendix  D 


1/ 


(bo  17) 


the  interface  pressure  may  be  expressed  as 


(5.18) 


When  V  »  ,  the  first,  term  of  Eq..  (5,2.8)  becomes  small 

compared  with  the  second,  leading  to  an  approximation  for  the 
interface  pressure 


P 


(5. 19) 


Values  of  are  all  fairly  close  to  1,5  for  the 

materials  considered  here.  Thus  the  logarithmic  term  in 
Eq,  (5„11)  is  a  representation  of  that  in  Eq0  (.5,16)  at  ex¬ 
tremely  high  velocities.  At  velocities  comparable  to  the 
speed  of  sound  c  ,  Eq,  (5. 19)  is  a  very  poor  approximation, 
in  fact,  the  first  term  in  Eq,  (5, 18)  becomes  dominant  for  V  <• 
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Thus  some  of  the  large  scatter  in  k.j  may  be  ascribed 
to  the  incorrect  use  of  Eq.  (5. 19)  at  low  velocities. 

In  the  derivation,  an  implicit  assumption  in 
Eq.  (5.17)  is  that  the  projectile  is  decelerated  as  a 
rigid  body.  Thus,  although  the  interface  pressure  is 
found  correctly  by  taking  the  compressibility  of  the 
materials  into  account  ,  the  transit  times  of  shock  and 
rarefaction  waves  across  the  projectile  are  neglected. 

This  is  probably  not  a  bad  assumption,  since  the  waves 
have  time  to  reverberate  through  the  projectile  many 
times  during  the  flow  period.  In  fact,  detailed  shock 
effects  in  the  projectile  are  also  smeared  in  Bjork*s 
calculations . 

However,  the  assumption  implicit  in  the  above 
derivation,  that  the  projectile  remains  undeformed,  is 
obviously  incorrect  at  high  velocities.  X-radiographs , 
and  Bjork’s  analysis  show  that  the  projectile  flows  so 
that  the  motion  is  more  nearly  spherically  symmetric 
than  linear. 

If  the  projectile  is  assumed  to  expand  radially 
as  an  incompressible  fluid,  the  interface  force  may  be 
represented  by 

F  =  kt  (p  f  kx  Ht)  2  ir  x*  (5.20) 

which,  upon  integration  of  the  equation  of  motion  Eq.  (5.I5) 
leads  to  the  penetration  law 


(.5*  21) 
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where  R  again  becomes  a  small  constant  at  high  velocities „ 

Thus,  it  may  be  expected  that  the  penetration,  law 
will  change  gradually  from  a  form  approximating  Eq0  (5*18) 
at  very  low  velocities  (where  R  may  be  expected  to  vary 
with  V  )  to  a  form  approximating  Eq„  (5„21)  at  very  high 
velocities. 

It  may  be  argued  that  since  high  velocity  effects 
are  of  prime  concern,  Eq,,  (5„21)  should  be  used  to  fit  the 
data  in  the  high  velocity  region  only,  disregarding  the  low 
velocity  data.  However,  present  experimental  data  in  the  high 
velocity  region  covers  at  most  velocities  varying  by  a  factor 
of  two.  This  is  insufficient  variation  to  warrant  a  two  para¬ 
meter  logarithmic  fit,  and  it  has  already  been  shown  that  the 
one  parameter  fit  Eq„  (fhl)  is  adequate,  over  this  velocity 
range 

Thus,  the  individual  fits  to  Eq„  (5.10)  must  be 
regarded  as  convenient  two  parameter  fits  over  the  whole 
experimental  velocity  range,  without  attaching  too  much  phy¬ 
sical  significance  to  the.  form  of  Eq.  (5,.  10)  „  The  fact  that 
there  is  some  correlation  of  k^  and  with  density  ratio 

must  be  ascribed  to  the  fact  that  there  is  also  some  correla¬ 
tion  between  dynamic  compressibilities  and  strength  effects 
with  density  ratio,  The  correlation  is  bound  to  be  rather 
poor,  since  it  may  be  observed  that  some  light  materials  such 
as  aluminum  may  have  high  strength,  while  some  heavy  materials 
such  as  lead  have  very  low  strength.  As  may  be  expected,  alu¬ 
minum  projectiles  striking  lead  targets  therefore,  show  very 
poor  correlation  with  other  material  combinations,  particularly 
in  Fig..  5.40..  Similar  remarks  may  be  made  for  other 

"anomalous"  material  combinations. 
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It  ■would  appear  to  be  impossible  to  find  a 
penetration  law  which  would  adequately  cover  both  a  wide 
range  of  materials  and  velocities  without  taking  the  correct 
geometrical,  compressibility  and  strength  effects  into  account. 
This,  of  course,  requires  a  full  theoretical  treatment  of  the 
problem o  It  would  therefore  be  appropriate  at  this  time  to 
carry  out  further  theoretical  work  before  attempting  further 
empirical  correlations. 


5.4.3  Comparison  with  Empirical  Power  Law 

When  the  material  combinations  exhibiting  extreme 
projectile  and  target  strength  effects  are  omitted,  and  the 
value  of  k-^  from  Eq.  (5.12)  is  used,  Eq.  (5. 11)  becomes 


/ 

4U 


(5°  22) 

O 

It  has  been  shown  that  the  function  y  =  k  log  (1+x  ) 

9  /  ® 

may  be  approximated  by  the  function  y  3  kx'J  over  the  range 
2=7  <  x  <  5  .  (Fig.  4=2)  Thus  Eq.  (5=22)  above  may  be  approxi¬ 
mated  in  the  high  velocity  region  by 


(5=23) 


107 


which  may  be  compared  with  Eq„  (5,1),  Comparisons  for  individual 
material  combinations  over  appropriate  velocity  ranges,  using 
correct  values  of  and  k2  in  Eq„  (5„10)  from  Table.  IV 

are  generally  even  better,  since  both  Eq„  (5„1)  and  (5„10) 
are  good  fits  to  the  high  velocity  data. 

Although  Eq,.  (5, 11)  or  Eq„  (5,22)  are  very  approxi¬ 
mate  due  to  the  poor  correlation  of  and  k2  with  material 

properties,  the  individual  fits  to  Eq,  (5° 10)  with  constants 
from  Table  V  are  very  good,  and  may  therefore  serve  as  very 
useful  empirical  expressions  suitable  for  interpolation  and 
prediction  purposes  over  the  whole  experimental  velocity  range, 

5, 5  Extrapolation  of  Empirical  Penetration  Laws 

Since  there  is  great  interest  in  penetration  effects 
at.  velocities  higher  than  those  Investigated  experimentally, 
it  is  interesting  to  investigate  the  consequences  of  extra¬ 
polating  the  empirical  penetration  laws,  Eq„  (5,,  1)  and  Eq„ 

(5.,  10')  Table  IV,  to  higher  velocities. 

Figure  5,42  shows  experimental  data  and  penetra¬ 
tions  predicted  by  Eq,  (5,1)  and  Eq„  (5,10)  for  aluminum 

projectiles  and  targets.  Original  constants  found  from  the 

2/q 

individual  least  squares  fit  to  p/d-  kV  ,  using  only 
points  in  the  high  velocity  region,  were  used.  For  the  lo¬ 
garithmic  penetration  law,  values  of  constants  k-^  and  k2 
found  by  individual  least  squares  fits,  and  listed  in  Table  V, 
were  nee d  logarithmic  scales  were  used  to  accomodate  a  large 

velocity  variation.  Also  shorn  are  the  three  points  obtained 

128 

theoretically  by  Bjork  It  should  be  noted  that  Bjork ‘s 

definition  of  equivalent  diameter  has  been  retained,.  Defini¬ 
tion  of  equivalent  diameter  on  the  basis  of  equivalent  mass 
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would  reduce  the  penetration  predicted  by  Bjork  by  16  per 
cento 

Firstly,  it  might  be  noted  that  both  penetration 
laws  are  good  fits  to  the  experimental  data  in  a  limited 
velocity  range..  This  velocity  range  will  now  be  defined  as 
the  high  velocity  region. 

Secondly,  it  may  be  seen  that  the  predictions  of 
penetration  by  the  two  penetration  laws  diverge  at  higher 
velocities.  The  logarithmic  law  predicts  lower  penetrations, 
and  over  a  wide  range  of  high  velocities  may  be  approximated 
by  a  law  of  the  type  p/d=  kV^3  o  This  velocity  range  will 
now  be  defined  as  the  hypervelocity  region . 

It  may  be  seen  that  the  logarithmic  law  shows 
surprising  agreement  with  the  theoretical  points,  computed 
by  Bjork,.  From  Fig.  5 .42,  it  is  seen  that  the  curve  fitted 
to  the  soft  HOOF  aluminum  data  points  fits  Bjork "s  points 
almost  exactly,,  In  view  of  the  fact  that  the  logarithmic 
penetration  law  has  limited  physical  significance,  this 
agreement  may  be  fortuitous .  It  may  be  noted  that  the  curve 
fits  Bjork' s  points  better  than  an  expression  of  the  form 
p/d  =  kV1/3  „ 

Furthermore  it  may  be  seen  that  the  logarithmic 
law  for  the  much  harder  2024-T3  aluminum  alloy  predicts 
lower  penetrations  than  Bjork' s  analysis,  the  offset  de¬ 
creasing  with  increasing  velocity.  (This  is  still  true  on 
a  linear  scale).  This  behaviour  is  exactly  what  might  be 
expected  for  a  high  strength  target. 

Figure  5„43  is  similar  to  Fig.  5„42,  but  for  steel 
projectiles  and  targets.  In  this  case,  agreement  is  not 
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quite  as  goocL  The  Tits  to  the  logarithmic  penetration  law 
were  obtained  from  points  limited  to  quite  low  velocities. 

The  steel  targets  had  high  strengths,  thus  all  of  the  data 
used  for  the  fits  were  in  the  transition  region  influenced 
by  projectile  and  target  strength.  The  transition  region 
strength  effects  tend  to  reduce  penetration,  the  effect 
being  greatest  at  low  velocities.  Thus,  the  tendency  will 
be  to  increase  the  initial  slope,  and  increase  the  curvature, 
of  the  empirical  fit.  The  two  effects  partially  offset  one 
another  as  may  be  seen  from  Fig,  5,43,  Agreement  of  the 
logarithmic  law  with  Bjork's  points  is  better  at  5,5  and  72 
km/sec,  than  at  20  km/sec.  Actually,  one  would  expect  that 
the  penetration  should  be  less  than  that  computed  by  Bjork, 
even  at  72  kra/sac,  particularly  for  the  very  hard  30  RC 
steel  alloy.  This  extrapolation  of  the  logarithmic  fits  in 
this  case  cannot  be  considered  to  be  realistic. 

The  same  effects  are  of  course  present  in  the 
logarithmic  fits  for  the  HOOF  and  2024-T3  aluminut  alloys, 
Fig,  5.. 42,  Due  to  the  low  hardness  of  the  HOOF  alloy,  and 
the  wide,  velocity  range  in  the  experimental  data,  relatively 
much  greater  reliance  may  be  placed  on  an  extrapolation  for 
this  material  and  the  agreement  with  Bjork's  calculations  may 
be  considered  to  have  some  significance  The  extrapolation 
for  the  2024-T4  alloy  however  must  be  considered  suspect, 
because,  although  data  to  quite  high  velocities  is  included, 
the  transition  region  also  extends  to  quite  high  velocities, 
and  transition  region  strength  effects  have  fairly  strong 
influence  on  the  fit 
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5.6  Definition  of  Regions  of  Impact 


Based  upon  the  information  in  the  previous  sections , 
it.  is  possible  to  postulate  regions  of  impact  within  which 
certain  types  of  behaviour  may  be  expected.  Since  defini¬ 
tion  of  the  last  two  regions  is  necessarily  based  on  the 
logarithmic  penetration  law  (F,q „  5.10),  which  has  no  firm 
theoretical  basis,  this  section  is  somewhat  hypothetical. 
Nevertheless,  the  behaviour  should  be  qualitatively  as 
described.  Quantitative  definitions  of  the  limits  of  the 
regions  may  be  refined  as  more  data  becomes  available  at 
high  velocities o 

5.6.1  Low  Velocity  Region 

This  region  is  only  important  for  projectiles 
which  have  higher  strength  than  the  target,  and  therefore 
suffer  little  or  no  deformation  at  low  velocity.  The 
region  is  of  no  interest  in  the  present  study,  and  will  not 
be  considered  here.  The  upper  limit  of  the  region  may  be 
expected  to  be  related  to  the  ratio  of  the  maximum  stress 
induced  in  the  projectile  P  ,  and  the  strength  of  the 
projectile,  represented  by  ,  but  it  is  likely  that 

the  projectile  shape  will  be  important  in  producing  stress 
concentrations  leading  to  fracture  of  brittle  projectile 
materials. 

5.6.2  Transition  Region 

Within  the  transition  region,  the  effects  of 
target  and  projectile  St*.  Sugtu  ale  in 

crater  shape  and  penetration  depth.  These  effects  have 
been  discussed  in  Section  5.3-6.  If  attempts  are  made  to 
fit  the  penetration  in  this  region  by  a  power  law. 
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Che  velocity  exponent  varies  rapidly  from  2  at  very  low 
velocities  down  to  2/3  at  the  high  velocity  end  of  the 
region.  Data  for  "ductile”  projectiles,  which  show  a  mono¬ 
tonic  increase  in  penetration  with  velocity,  may  be  fitted 
empirically  by  a  two  parameter  logarithmic  fit  Eq„  (5,10),. 
This  is  not  necessarily  physically  significant . 

The  region  extends  from  the  velocity  at  which  the 
projectile  suffers  serious  deformation  to  a  velocity  such 
that 


or  if  H  >  Ht 


(See  Section  5,3.6)  i.e.  a  velocity  at  which  the  initial 
interface  pressure  is  much  greater  than  the  yield  strength 
of  both  the  projectile  and  target. 

5 „  b o ,3  II Igh  Veloc  ity  Region 

Within  this  region,  strength  effects  are  no  longer 
dominant,  but  are  still  not  negligible. 

Crater  shapes  approach  hemispheres,  but  strength 
effects  are  sufficient  to  cause  craters  formed  by  low  or  high 
density  projectiles,  compared  to  the  target  density,  to  be 
shallower  or  deeper  than  hemispheres  respectively.  Cratering 
efficiency  V  /E  is  approximately  constant  within  this  region. 

Within  this,  region,  the  logarithmic  penetration 
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where  X  =---  p  /pfc  is  the  density  ratio  and 
is  the  Best  Number,  may  be  approximated  by 
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or  Eq„  (5.1).  The  two  equations  may  be  seen  to  be  equivalent 
roughly  between  the  limits 
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from  Fig.  4.2.  Particularly  for  high  strength  projectiles, 
the  limit  Eq.  (5.4)  may  occur  at  a  higher  velocity  than 
Eq.  (5.25)  above  and  Eq.  (5-4)  should  therefore  be  -used  as 
the  lower  velocity  limit. 


5.6.4  Hyperveiocity  Region 


Within  this  region,  the  strength  effects  are 
negligible,  and  Bjork's  theory  may  be  expected  to  be  valid. 
Craters  may  be  expected  to  be  hemispheres.  Within  this 
region,  the  logarithmic  penetration  law  may  be  approximated 
by 
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i.e.  the  penetration  is  approximately  proportional  to  V  ' 
From  Fig.  4.2  the  approximate  limits  of  the  approximation  are 


6  4-  < 


*  Vj  a 
K 


*  3  i  00 


(5-27) 


113 


However,  additional  effects  not  considered  previously 
may  enter  strongly  into  this  region.  In  particular,  target 
melting  may  become  important.  It  is  difficult  to  assess  the 
magnitude  of  this  effect  at  the  present  time,  but  it  may  be 
expected  to  lead  to  an  increase  in  penetration  and  crater 
size.,.  Thus,  the  above  remarks  concerning  the  hypervelocity 
region  must  be  considered  speculative  at  this  time. 

The  Best  Number  may  be  seen  to  have  the  meaning 

of:  the.  ratio  of  a  pressui*e  (p.  V'  )  to  material  strength  (H.)  „ 

**-  2/^ 

It  can  be  seen  that  at  extremely  high  velocities  X  /JptV '^2 
approximates  the  initial  interface  pressure.  P"  for  typical 
values  of  given  in  Table  IV.  (See  Section  5.4,3)  It 

may  be  preferable  to  define  a  non-dimensional  parameter  P  /II 
and  find  least  squares  fits  to  Eq„  (5,,  1.6) 

J  =  >'  1 0  it  (  1  +  ~Hf.  )  (5-28) 

This  would  permit  redefinition  of:  the.  limits  of.  the.  regions 
of  impact  Eq.  (5,25)  and  (5,. 2.6)  in  terms  of  p‘/Ht  ,  which 
will  be  preferable  from  the  point  of  view  of  physical  inter 
pretation,. 

b, 7  Discus s ion 

A  quite  detailed  qualitative  understanding  of 
cratering  has  been  gained  from  experimental  studies  using 
X-ray  and  high  speed  optical  photography  and  other  active 
measuring  techniques,  and  from  detailed  theoretical  analyses 

i  o  ft 

such  as  that  carried  out  by  Bjork  „  When  data  from  a  large 
variety  of  material  combinations  and  over  a  wide  velocity  range 


L 1  ^ 


are  assembled  and  compared,  the  effects  of  the  complicated 
phenomena  accompanying  penetration  on  the  crater  depth  and 
shape  can  be  seen. 

It  is  not  surprising  that  simplified  theories, 
presented  previously,  are  unsuccessful,  since  in  each  case 
they  neglect  many  phenomena  important  Co  the  cratering 
process. 

Even  Bjork's  theory,  which  correctly  takes  geo¬ 
metrical  and  compressibility  effects  into  account,  cannot  be 
considered  to  be  realistic  in  the  high  velocity  region,  ex- 
cept  for  low  strength  alloys,  since  the  shear  strength  of 
the  materials  is  neglected.  The  velocity  region,  in  which 
strength  effects  are  in  fact  negligible,  is  evidently  much 
higher  than  previously  supposed. 

It  is  also  not  surprising  that  so  much  confusion 
and  contradiction  exists  among  previous  empirical  fits 
derived  from  limited  sections  of  the  data,  covering  different 
limited  ranges  in  material  properties  and  velocities,  since 
nearly  all  of  the  experimental  data  lies  in  the  transition 
region  where  strength  effects  lead  to  such  complex  behaviour, 
and  even  the  limited  amount  of  data  in  the  high  velocity 
region  is  affected  by  material  strength  to  some  extent. 

Two  empirical  penetration  laws  have  been  derived 
in  this  study  which  are  more  generally  applicable  than 
previous  empirical  expressions.  Each  is  strictly  limited 
in  the  range  of  parameters  covered,  and  in  physical  signi¬ 
ficance.  This  is  a  property  generally  associated  with  em¬ 
pirical  expressions. 
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The  empirical  power  law 


p  /  \  ,  „ 

(  -  /Con  ±  <o ■  07 j  m  5  (5o28) 

2 

where  A  =»  p  /p^  1?.  the  density  ratio  and  B  =  p£V  /Ht  is 

the  Best  Number ,  is  limited  to  the  high  velocity  region  de¬ 
fined  in  Section  5 , h  „  h  .. 

On  the.  other  hand*  the  empirical  logarithmic  law 

£  =  k,  toU  [  /  +  r,  ]  (5"29) 

with  constants  kj  and  tabulated  in  Table  V  may  be 

expected  to  fit  only  those  material  combinations  for  which, 
and  are.  determined  empirically,  hut  over  the  entire 

experimental  velocity  range,,  The  approximate  equivalence, 
of  Eq„  (5.28)  and  Kq.  (fj. 2.9)  where,  their  ranges  of  validity 
coincide,  has  been  demonstrated  in  Section  5..<4«3  and  5.5., 

Extrapolation  of  the  empirical  penetration  laws 
Eq.  (  5.  2.8)  and  l, 5. .29)  ,  either  to  other  materials  or  to  higher 
velocities,  has  little  theoretical  or  other  justification. 

in  particular  there  is  no  justification  whatever 
for  extrapolating  the  empirical  power  law  Eq„  (3,. 28)  to 
higher  velocities The  argument  has  been  advanced  that  the 
cratering  efficiency  V  /E  will  remain  constant  as  the 
velocity  increases..  This  in  effect  implies  that  the  partition 
of  energy  in.  the  various  modes  such  as  heating,  energy  trans¬ 
mitted  to  other  parts  of:  the  target  by  the  stress  wave  system, 
kinetic  energy  carried  away  by  the  material  thrown  out  of  the 
crater,  and  energy  dissipated  In  plastic  deformation,  remains 
constant.,  Such  an  assumption  cannot  be  substantiated  either 
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from  experimental  evidence,  or  from  theoretical  considerations, 
at  the  present  time. 

It  appears  to  be  fortuitous  that  the  transition 
region  strength  effects  act  in  such  a  way  as  to  partially 
cancel  each  other,  so  that  even  in  the  case  of  high  strength 
steel  targets  the  extrapolated. logarithmic  law  Eq.  (5„29) 
shows  reasonable  agreement  with  the  results  of  Bjork's  cal¬ 
culation,  up  to  72  km/seco ,  Fig*  Bjork's  theory  may 

be  expected  to  provide  an  upper  limit  to  the  penetration, 
since  strength  effects  neglected  in  the  theory,  but  operative 
in  the  high  velocity  region,  may  be  expected  to  reduce  pene¬ 
tration.,  Thus  the  logarithmic  law  Eqc  (5.29)  may  be  expected 
to  overestimate  penetration  at  velocities  below  72  km/sec, 
and  to  underestimate  penetration  above  72  km/ sec,  for  the 
particular  steels , investigated „  No  general  conclusions  about 
extrapolations  for  other  materials  can  be  drawn,  and  in  fact 
some  other  materials  may  well  behave  in  such  a  manner  as  to 
lead  to  very  large  errors  when  Eq„  (5„29)  is  extrapolated  to 
high  velocities „ 

Extrapolations  can  only  be  carried  out  with  some 
confidence,  when  the  penetration  law  is  based  on  sound  theo¬ 
retical  grounds,  and  there  is  some  assurance  that  phenomena 
neglected  in  the  theory  do  not  become  operative  at  higher 
velocities . 

The  comparative  success  of  the  two  penetration 
ietws  £q„  and  ,  based  on  elementary  material 

properties,  must  be  ascribed  to  the  fact  that  there  exists 
a  rough  correlation  between  dynamic  compressibilities  and 
strength  properties,  and  static  properties. 
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A  penetration  law  to  cover  both  a  wide  range  in 
material  properties  and  velocities  necessarily  involves  a 
much  more  complex  functional  dependence  of  penetration  on 
material  properties  chan  Eq„  (5,28.)  and  (5* 29) «  It  appears 
to  be  a  hopeless  task  to  deduce  such  a  functional  dependence 
from  the  present  penetration  data,  both  because  the  present 
data  has  large  uncertainties  associated  with  it,  and  because 
the  data  covers  insufficient  variation  in  material  combina¬ 
tions  and  velocities o  In  particular,  very  few  material  com* 
binations  have  been  investigated  in  the  high  velocity  region., 

It:  seems  appropriate  therefore  that  further  theo- 
ratical  work  be  undertaken  in  order  to  guide  the  formulation 
of  more  realistic  penetration  laws,,  The  understanding  of 
the  penetration  process,,  provided  by  Bjork's  analysis  cannot 
be  underestimated „  It  Is  to  be  strongly  urged  that  theoreti¬ 
cal  work  be  directed  towards  extending  an  analysis  of  the 
type  used  by  Bjork  to  include  shear  strength  of  the  target 
and  projectile,,  Since  melting  may  become  extensive  at  very 
high  velocities,  and  will  have  the  effect  of  reducing  the 
shear  strength,  this  should  also  be  considered,. 

Such  an.  analysis  will  provide  quantitative  infor¬ 
mation  about  the  various  strength  effects  operative  in  both 
the  transition  and  high  velocity  regions,  and  thus  aid 
in  interpretation  of  experimental  results  and  guide  formula¬ 
tion  of  realistic  penetration  laws,.  Furthermore,  such  an 
analysis  will  provide  a  more  rational  basis  for  extrapolation 
to  velocities  higher  than  those,  which  may  be  investigated 
e  xpe  r i mental iy . 


SECTION  VI 


OTHER  PROJECTILE-TARGET  CONFIGURATIONS 


6.1  Microparticle  Impact 


Although  penetration  appears  to  scale  quite  well 
with  projectile  size  over  the  range  of  projectile  sizes 
conveniently  fired  in  guns,  or  by  cavity  charge  techniques, 
there  is  a  question  if  this  scaling  holds  true  for  very 
small  projectiles,  approaching  the  size  of  the  grain  struc¬ 
ture  . 


While  several  laboratories  have  been  actively 
engaged  in  developing  projection  techniques  for  micro¬ 
particles  only  three  groups  have  reported  cratering 
results  by  projectiles  under  200  micron  diameter-  Anderson 
et  at  Stanford  Research  Institute  projected  100  - 

150  micron  stainless  steel  spheres  at  velocities  up  to 
12,500  ft/sec.  against  a  variety  of  targets-  The  particles 
were  individually  projected  by  cylindrical  high  explosive 
charges.  (Fig.  6.1)  Friichtenicht  et  al^  at  Ramo  Wooldridge 
used  an  electrostatic  acceleration  technique  to  project  1  - 

10  micron  carbonyl  iron  particles  at  velocities  up  to  6000 

2  8  10 

ft/sec  against  lead  targets.  (Fig-  6-2)  Gehring  et  al  '  * 
at  Ballistics  Research  Laboratory  used  a  specially  designed 
shaped  charge  with  cast  iron  liner  to  project  a  cloud  of 
particles  of  50  to  200  micron  size,  a  special  shutter  being 
used  to  allow  only  particles  with  a  velocity  in  the  region 
of  about  32,000  ft/sec  to  reach  the  target-  The  size 
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distribution  of  resultant  craters  (about  500  per  target)  was 
then  correlated  with  the  size  distribution  oi:  particles  obtained 
by  crushing  a  similar  cast,  iron  liner  with  pestle  and  mortar - 
Although  several  target  materials  were  used,  only  crater  dimen¬ 
sions  for  lead  and  copper  targets  have  been  reported  at  this 
time.  Grater  diameters  were  measured  -  However,,  based  on 
Gehring's  observation  that  the  craters  were  nearly  hemispher¬ 
ical,  corresponding  penetrations  are  plotted  in  Fig..  6»3„ 

The  microparticle  results  show  a  large  scatter, 
associated  partly  with  the  experimental  difficulties  of 
accurately  measuring  projectile  and  crater  sizes  and  pro¬ 
jectile  velocity.,  However  some  scatter  is  probably  due  to 
the  effect  of  grain  boundaries  and  local  anisotropies  in 
the  target  material.  Gehring^  has  shown  contours  of  craters 
produced  in  coarse  grained  copper  and  copper  single  crystals 
by  200  micron  cast  iron  projectiles  .  The  penetration  appar¬ 
ently  is  affected  to  some  extent  by  grain  orientation,  and 
the  presence  oi  grain  boundaries.,  Anderson*31'1  has  observed 
similar  effects 

On  each  of  Fig.  6.  1  to  6., 3  curves  representing 
28 

Charters  and  Locke * s'  empirical  expression 


have  been  inserted..  This  expression  adequately  fits  penetra¬ 
tion  data  in  Lead  and  copper  targets  tor  large  projectiles  in 
the  velocity  range  ot  interest  here  The  microparticles  results 
lie  consistently  below  the  empirical  curve-  However,  the  dis¬ 
crepancy  is  probably  mainly  due  to  the  fart  that  all.  of  the 
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microparticle  craters  were  measured  to  the  top  of  the 
projectile  material  remaining  in  the  crater  (except  Gehring's 
results,  in  which  crater  diameters  were  measured).  Anderson’s 
results  further  may  be  affected  by  mass  lost  by  the  particles 
during  explosive  acceleration  and  subsequent  ablation. 

One  may  conclude  that  within  the  present  large 
experimental  uncertainty  it  is  impossible  to  distinguish 
if  microparticle  cratering  follows  a  different  penetration 
law. 


6.2  Oblique  Impact 

Cratering  due  to  projectiles  approaching  the 
target  at  oblique  incidence  is  cf  great  practical  interest. 
However,  oblique  impact  has  received  relatively  little 
attention. 

In  the  high  velocity  region,  in  which  nearly  hemi¬ 
spherical  craters  are  formed!  in  normal  impact,  several 
experimenters  have  noted  that  nearly  hemispherical  craters 
are  also  formed  in  oblique  impact  up  to  some  critical  angle. 
(Fig.  6.4)  Beyond  the  critical  angle,  the  crater  becomes 
asymmetrical,  with  greater  depth  and  steeper  sides  towards 
the  direction  from  which  the  projectile  approached.  At 
grazing  incidence,  the  crater  becomes  very  elongated,  and 
it  has  been  noted  that  the  projectile  ricochets,  and  is 
capable  of  producing  further  craters  in  adjacent  targets. 
Spray  patterns.  Fig.  6.5,  and  shock  patterns.  Fig.  6.6, 
also  show  symmetry  below  the  critical  velocity,  asymmetry 
above.  The  critical  velocity  is  clearly  a  function  of 
velocity.,  Fig.  6,7  shows  craters  in  lead  produced  by 
steel  projectiles.  At  7,000  ft/sec,  a  symmetrical  crater 
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is  formed  at  30°  incidence.,  but  not  at  60°  incidence 
An  increase  in  velocity  to  9800  ft/sec  however  results  in 
a  symmetrical  crater  at  60°  incidence.  No  systematic  experi 
ments  nave  been  carried  out  to  determine  the  dependence  of 
the  critical  angle  ot  incidence  on  velocity  tor  various 
materials  * 

30 

Data  on  cratering  has  been  reported  by  Summers 
iox  copper  projectiles  impacting  copper  targets  at  7 .,000  and 
11,000  it/sec..  ,  (Figo  6  ,8),  and  by  Kineke*’  ,  (Fig.  6  9)  ,  tor 
steel  discs  impacting  lead  targets  at  16,400  it/ sec.  Both 
experimenters  noted  that;  the  data  tor  oblique  impact  com¬ 
pared  very  well  with  that  for  normal  impact,  it  penetration 
versus  the  normal  component  of  velocity  is  plotted.  The 
oblique  impact  data  departs  from  the  normal  impact  data  at 
about  the.  same  critical  angle  at  which  the  crater  ceases  to 
be  hemispherical. 

r  m 

Anderson  ,  (Fig..  6.10).,  also  report,  microparticle 
cratering  data  at.  oblique  angles.,  The  scatter  is  such  that 
it  is  not  possible  to  determine  it  the  60°  and  30  incidence 
data  departs  from  the  normal  incidence  data  when  plotted  on 
the  basis  of  normal  component  of  velocity.. 

In  the  high  velocity  region  (Sect..  5-1)  it  has 
already  been  noted  that  crater  volume  is  approximately 
proportional  to  projectile  energy  at.  normal  incidence,.  Thus 
below  the  critical  angle  one  might  expect  that  crater  volume 
would  be  proportional  ro  the  normal  component.  u.i  projectile 
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energy  However  Bryan  '  has  noted  that  the  ratio  of  crater 
voLume  to  projectile  energy  appears  to  be  o.  linear  function 
ot  cos  /’  .,  1  e 

r,v  +  b 

r  • 

c 

•  1  > 
f.e.  d 
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to  angles  beyond  the  critical  angle.  (Fig.  6.11) 
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Bryan  has  also  given  a  simplified  theory  for 
predicting  the  ratio  of  major  to  minor  diameter  of  the  crater 
mouth  for  angles  of  incidence  beyond  the  critical  angles. 

It  is  assumed  that  the  projectile  simply  generates  a  spher¬ 
ical  hole  concentric  with  itself,  which  expands  at  constant 
velocity  .  Choosing  moving  co-ordinates  such  that  the 
projectile  approaches  the  target  vertically  at  velocity 

V  cos  p,  and  the  target  moves  laterally  at  velocity  V  «* 

s 

V  sin  p,  the  crater  will  cease  to  be  spherical  when  V  >  V 

S  X 

The  ratio  of  major  to  minor  axis  is  then  simply 


P 


I/5  *  Vr 

Z  Vr 


Assuming  V  “  k  V  cos  p 

|  /  t  a  r>  4* 

P  *  S'  (  — 
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and  the  critical  angle  at  which  the  crater  begins  to  elongate 
is  (at  p  =  1) 

k  =  tan  (f)Q  (6.3) 


A  refinement  was  introduced  to  allow  for  the 
finite  size  of  the  projectile.,  of  radius  R  s  which  by 
similar  reasoning  results  in  the  equation 
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where  D  is  the  minor  diameter  of  the  final  crater.  The 
cm 

value  of  0  must  be  supplied  by  experiment. 

The  theoretical  curves  for  both  cases  are  plotted 

in  Fig,,  6.12„  Shown  also  are  experimental  points  obtained 
5 

by  Kineke.  Each  point  represents  an  average  of  10  firings 
using  1/2  inch  diameter  by  0.040  inch  thick  steel  discs 
accelerated  by  air  cavity  changes  to  16 „ 400  ft/sec,  the  angle 
of  incidence  being  varied  from  0°  to  69°  „  The  value  of 
the  critical  angle  used  in  the  experiments  was  obtained  from 
Kineke's  experimental  data.  The  projectiles  spalled  during 
acceleration,  and  a  diameter  of  0.3  inches  was  used  in  the 
theory  The  experimental  points  lie  between  the  curves  ne¬ 
glecting  the  projectile  size,,  and  taking  the  projectile  size 
into  account,.  Since  the  projectiles  tumbled  in  flight.,  and 
the  attitude  at  the  time  of  impact  was  not  known,  such  scat¬ 
ter  may  be  expected,  and  the  agreement  is  therefore  quite 
encouraging 

0.3  Rod  Impact 

Relatively  Little  experimental  work  has  been  done 
with  projectiles  ot  length  to  diameter  ratios  much  greater 
chan  unity,,  despite  the  evident  advantage  tor  offensive 
systems.  Figure  6  13  graphically  illustrates  the  greatly 
increased  penetration  per  unit  mass  obtained  by  rods  (moving 
in  the  direction  ot  the  longitudinal  axis)  as  compared  to 


12 '4 


penetration  obtained  by  spheres  and  cylinders  of  unit 
length/diameter  ratio.  Unfortunately  direct  comparison  of 
penetration  in  the  two  cases  is  uncertain  since  different 
aluminum  alloys  were  used  as  targets,  and  rod  velocities 
were  rather  low . 

No  systematic  study  of  penetration  as  a  function 
of  length/diameter  ratio  has  been  carried  out,,  and,  in  fact, 
at  this  time,  results  of  only  three  studies  have  been  pub¬ 
lished.  Summers  and  Niehaus^"  used  rods  of  L/d  between 

6  and  10  at  velocities  up  to  10.000  ft/sec.  against  lead, 

48 

copper  and  steel  targets.  Allen  and  Rogers  used  rods  of 

L/d  between  6  and  12  at  velocities  up  to  10,000  ft/sec 

23 

against  7075-T6  aluminum,  and  Slattery  and  Clay  used 
rods  of  L/d  between  8  and  15  at  a  constant  velocity  of 
8,. 500  ft/sec  against  aluminum.  Results  are  plotted  in 
Fig o  6  o 14  through  6.20. 

The  phenomena  accompanying  rod  impact  differ 
somewhat  from  those  accompanying  impact  of  projectiles  ot 
about  unit  L/d  ratio.  It  may  be  expected  that  the  shock 
propagating  upwards  into  the  projectile  is  rapidly  attenu¬ 
ated  by  rarefactions  from  the  unsupported  sides  o.t  the  rod., 
which  lead  to  rapid  lateral  flow  of  rod  material.  Thus  the 
disturbance  travelling  back  into  the  rod  would  be  expected 
to  move  at  sonic  velocity.  Two  types  of  behaviour  may  then 
be  expected  according  to  whether  the  disturbance  is  carried 
above  the  original  target  surface.,,  i „ e ,  the  impact  velocity 
is  below  sonic  velocity,  or  the  disturbance  is  carried  below 
the  original  target  surface,  i.e.  supersonic  impact  velocity 
In  the  latter  case,  the  flow  after  an  initial  transient 
period  may  be  expected  to  be  almost  pseudo-stationary  until 
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all  of  the  rod  material  has  been  reached  by  the  disturbance., 
There  would  then  be  a  final  transient  period  during  which 
the  projectile  and  target  material  move  under  the  action  of 
inertial  forces,,  resisted  in  the  latter  stages  by  the  strength 
of  the  material  involved., 

The  difference  between  rod  and  compact  projectile 
behaviour  therefore  lies  principally  in  the  existence  of  the 
intermediate  pseudo- steady  flow  period  at  supersonic  velocities., 
which  is  absent  .in  the  case  compact  projectiles „  (See  Section 

5.1) 

23 

Slattery  '  has  noted  different  crater  shapes, 
depending  on  whether  the  impact  velocity  was  above  or  below 
the  sonic,  velocity  in  the  projectile.  Subsonic  velocities 
were  observed  to  lead  to  cone  shaped  craters,  while  super¬ 
sonic  velocities  lead  to  cylindrical  craters.,  Such  behaviour 
may  he  expected..,  since  in  the  subsonic  case,,  the  disturbance 
reaches  the  back  end  of  the  rod  and  is  able  to  decelerate 
the  rod  The  latter  part  of  the  penetration  process  then 
occurs  at  lower  velocity..  Hence  lower  lateral  velocities 
are  induced,  and  less  expansion  of  the  crater  may  be  expected,, 
Figure  6 .21  shows  a  crater  in  which  the  impact  velocity  was 
so  low  that  the  rod  was  decelerated  to  a  velocity  where  material 
strength  halted  the  flow  before  gross  deformation  of  the  back 
end  of  the  rod  occurred 

48 

Allen  and  Rogers  attempted  to  use  shaped  charge 

jet  theory.,  (Section  3.2)  essentially  fitting  their  data  to 

2 

Eq~  (3.34)  by  suitable  choice  of  the  constant  k  «=  1/2  0  V 

_  P  o 

The  semi -empirical  curves  are  inserted  in  rig.  6,14  and  6,15., 

The  value  of  the  constant  k  was  found  to  compare  well  with 
expected  yield  strength  of  the  7075-T6  aluminum  alloy  under 
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dynamic  conditions. 

In  the  case  of  gold  rods.,  the  penetration  was 
found  to  exceed  the  predicted  values  at  the  higher  velocities. 
It  was  found  that  the  assumed  residual  velocity  of  the  projec¬ 
tile.,  (2v  -  V)  where  v  is  the  calculated  velocity  of  the 
target-projectile  interface, was  greater  than  vq  ,  the 
velocity  below  which  no  penetration  occurs  due  to  target 
strength.  Thus  it  was  assumed  that  secondary  penetration 
was  caused  by  the  projectile  material,  of  magnitude  pre¬ 
dicted  by  shaped  charge  jet  theory,  Eq.  (3.34),  with  initial 
projectile  velocity  taken  as  (2v  -  V).  Results  of  this 
calculation  are  inserted  on  Fig.  6.15. 

31 

On  the  other  hand  the  data  of  Summers  and  Niehaus  " 
in  all  cases  lie  well  above  semi-theoretical  curves  derived 
in  the  same  manner,  and  in  fact,  at  the  higher  velocities., 
lie  well  above  the  penetration  predicted  by  simple  jet 
theory  p/L  Jp^  /Pt  which  sets  an  upper  limit  to  Eq.  (3.34). 
See  Fig.  6.17  to  6 . 19 .  Secondary  penetration  which  is  only 
possible  at  these  velocities  in  the  case  of  tungsten  car¬ 
bide  into  copper.  Fig.  6.18,  and  tungsten  carbide  into  lead. 
Fig.  6.19,  is  insufficient  to  account  for  the  large  discre¬ 
pancies,  Thus  it  is  premature  at  this  point  to  place  re¬ 
liance  on  jet  theory.  Figure  6.20  indicates  that  there  is 

little  correlation  on  the  basis  of  density  ratio  at  a  con* 

iB 

stant  velocity  of  8,500  ft/sec  when  the  results  of  Allen 
23 

and  Slattery  are  compared. 

/,u  23 

Both  Allen  and  RogerstTU  and  Slattery  and  Clay^* 
observed  that  the  penetration  of  rods  was  affected  very 
little  by  the  angle  of  obliquity  of  the  target,  provided 
the  rod  remained  unyawed  (i.e.  moved  in  the  direction  of 
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its  axis) .  However  Slattery  and  Clay"  have  shown  that  the 
penetration  is  critically  dependent  on  the  angle  of  yaw. 
Figure  6„22  shows  a  typical  crater  produced  by  a  yawed  rod 
in  aluminum.  The  section  was  taken  in  the  plane  of  the  yawed 
rod..  At  right  angles  of  this  section  the  crater  is  two  or 
three  times  the  rod  diameter  in  width,  and  is  thus  a  narrow 


elongated  slot. 


Slattery  and  Clay  have  suggested  a  simplified 
explanation  for  the  crater  shape 


The  total  length  of  the  rod  does  not  contribute  to  the 
penetration,  The  crater  expands  radially.,  and  thus  is  inclined 
as  shown,,  When  the  rod  has  penetrated  to  a  depth  such  that 
the  rod  intersects  the  crater  llp,»  it  is  assumed  that  primary 
penetration  ceases.  Tbe  remainder  of  the  rod  then  presumably 
continues  to  excavate  material  from  the  intersected  side  of 
the  crater.,  The  length  of  rod  active,  in  primary  penetration 
depends  on  the  radius  of  the  crater.  The  following  assumptions 


Private  Communication 


128 


were  made  to  estimate  this  radius 

a)  the  interface  moves  at  constant  velocity  v  . 

b)  the  crater  on  any  plane  parallel  to  the  surface 
is  cylindrical  of  radius  r  ,  centred  at  the 
point  at  which  the  rod  penetrated  the  plane. 

c)  the  volume  of  the  crater  is  proportional  to  the 
kinetic  energy  of  the  rod. 

For  these  assumptions  the  radius  of  the  crater  r  will 
be  given  by  the  equation 

TTT  K' fi,  (y-*)  *  k4  IT  rx  V  (6.5) 

where  R  is  the  radius  of  the  rod,  and  Kt  is  the 
energy  volume  ratio  E/V  found  for  compact  particle 

v 

impact.  Values  given  by  Feldman  (2.6  Hfc)  were  used, 
but  results  have  also  been  computed  using  K found  in 
Section  5*3.  <3.06  (pt/pp)2/3Ht) 

If  the  interference  occurs  at  a  depth  p  ,  then 
the  portion  of  the  rod  already  in  the  crater  will  continue 
to  penetrate.  The  final  depth  of  the  crater  is  then 

r>f  *  p  (}  +  vh)  <6°6) 

Further  for  a  constant  yaw  angle  0  ,  from  the 

diagram 

<P?  (l-  ~)  *  <--  R  (6.7) 
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which  using  Eq.  (6.5)  and  (606)  becomes 


n 


K  _!!_  /  v  EE_  v-v  ' 

I j)  V-v  l  j2Kt  * 


In  normal  impact 


V” 

]7^7 


(6.8) 


(6.9) 


so  that  Eq.  (6.8)  may  be  evaluated  using  empirical  values 
of  p/L  found  in  normal  impact  experiments. 

The  crater  shows  a  residual  penetration,  Fig.  6.22^ 
assumed  to  be  due  to  the  cross  hatched  part  of  the  rod  in  the 
diagram  above.,  of  length  2R//  „  This  is  assumed  to  lead  to 
an  added  penetration 


2.  R  is 

<P  V  —  v 


(6.,  10) 


so  that  the  total  penetration  is  pfc  =  pf  +  ,  A  compari¬ 

son  of  measured  penetrations  and  penetrations  calculated 
on  the  basis  of  Feldman's  and  that  given  in  Section  5.3 

are  shown  in  the  table  below.  Values  of  p/I>  for  normal 
impact  have  been  taken  from  Ref.  23, 
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Calculated 

Calculated 

Rod 

Yaw 

(Feldman) 

(Sect 

.  5.3) 

He  a 

sured 

Material 

Pt 

pf 

Pt 

Pf 

Pt 

Pf 

0  »  40 

0.469 

0.217 

0.600 

0.277 

0.344 

0 . 125 

6° 

0.312 

0,145 

0.400 

0.185 

0.250 

0.156 

Lead 

o 

o 

I— 1 

0.188 

0.087 

0.240 

0.112 

0.187 

0 . 152 

G 

Oi 

o 

0.573 

0.229 

0.631 

0.252 

0.500 

-- 

6° 

0.287 

0.115 

0.316 

0.127 

0.250 

0.125 

Copper 

10° 

0.173 

0..069 

0.190 

0.076 

0.156 

0.120 

Penetration  Depths  in  Inches 
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6.4  Finite  Thickness  Targets 


Cratering  and  penetration  of  thin  targets  by 
high  velocity  projectiles  has  not  received  as  much  atten¬ 
tion  to  date  as  cratering  in  quasi- infinite  targets.  While 
a  number  of  studies  have  been  reported,  none  of  these  are 
directly  comparable  since  different  experimental  conditions 
were  used*  However,  sufficient  work  has  been  done  to  dis¬ 
tinguish  some  important  phenomena  in  thin  target  cratering. 


u»4o  i 


Relatively  Thick  Targets 


Cratering  in  relatively  thick,  targets  which  are 
not  completely  penetrated  may  give  rise  to  damage  at  the 
reap  surface  of  the  target  plate.  (Fig.  6.23.)  This  damage 
may  range  from  the  appearance  of  a  zone  of  slipped  and 

34 

deformed  grains  near  the  rear  surface  opposite  the  crater 

to  the  appearance  of  internal  cracks  and  bulging  of  the 

rear  surface,  to  complete  detachment  of  a  segment  of  target 

12 

material,  termed  spall  „  The  damage  is  undoubtedly  due 

to  the  reflection  and  interaction  at  the  rear  surface  of 

1.2 

the  wave  propagated  into  the  target.  Maiden  et  al  report 
spall  thickness  as  a  function  target  thickness  at  constant 
projectile  velocity  and  as  a  function  of  projectile  velocity 
at  constant  target  thickness  tor  steel  targets  at  velocities 
up  to  13,000  ft/sec. 

The  presence  of  damage  at  the  rear  surface 
generally  is  accompanied  by  a  greater  crater  depth  than  in 
a  corresponding  quasi. -Infinite  target.  Kinard  et  al'^  have  re¬ 
ported  a  study  of  the  effect  of  target  thickness  on  penetra¬ 
tion  of  steel  and  aluminum  projectiles  into  aluminum  targets 
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at  impact  velocities  between  5,000  and  13,000  ft/sec. 

They  conclude  that  for  the  velocities  used  in  their 
study,  the  penetration  is  affected  by  target  thickness 
provided  that  the  penetration  is  greater  than  20  percent 
of  the  target  thickness .  It  was  found  that  a  projectile 
could  completely  penetrate  a  target  whose  thickness  is 
approximately  one  and  one  half  times  as  great  as  the 
penetration  of  a  similar  projectile  into  a  quasi- infinite 
target*  Some  results  of  this  study  are  plotted  in  Fig*  6*24* 

6*4,2  Thin  Targets 

Thinner  targets  which  are  completely 
penetrated  exhibit  different  behaviour  depending  on  the 
projectile  velocity,  target  thickness,  and  materials 
involved. 

At  relatively  low  velocity,,  penetration  is 
accomplished  by  modes  familiar  in  armor  penetration  work. 

The  resultant  hole  is  nearly  the  same  size  as  the  projectile, 
which  suffers  very  little  deformation  in  the  penetration  pro¬ 
cess,  ^  As  the  impact  velocity  increases,  the  projectile 
suffers  increasingly  severe  damage;  plastic  deformation  and 
''mushrooming"  if  ductile,  or  small  pieces  being  spalled 
from  the  projectile  if  brittle.  The  resultant  hole  in  the 
target  plate  consequently  becomes  appreciably  larger  than 
the  original  diameter  of  the  projectile,  and  may  become 
conical,  increasing  in  diameter  towards  the  rear  surface 
as  the  projectile  continues  to  mushroom,  or  spall  during 
the  penetration  process  .  As  the  impact  velocity  in¬ 
creases,  more  aevere  breakup  or  flow  of  the  projectile  is 
observed  with  consequent  further  enlargement  of  the  target 
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hole,  until  at  relatively  high  velocities,,  both  target 
material  and  projectile  material  are  ejected  from  the  rear 
of  the  target  in  a  spray  of  small  particles  distributed  over 
a  solid  angle  of  90  degrees  or  more,  with  a  considerable 
amount  of  material  alto  being  ejected  from  the  front,  of  the 
target  plate^' 5  ,  The  resultant  hole  in  the  target  plate 
is  typically  several  times  the  diameter  of  the  original  pro¬ 
jectile.  These  effects  are  illustrated  in  Fig.  6.25  and 
6„26  by  frames  from  high  speed  framing  camera  records  of  the 
penetration  process. 

These  phenomena  are  evidently  analogous  to  those 
occuring  in  semi-infinite  targets  and  have  led  several  inves- 
tigato-'s  to  distinguish  similar  regions  of  impact,  namely; 
low  velocity  region  (undetormed  projectile,  armor  penetration); 
transition  region;  and  high  velocity  region,  (hypervelocity 
or  fluid  impact  region) .  Several  dif ferent  types  of  experi¬ 
ments  have  been  carried  out,  as  illustrated  in  Fig.  1.1  d,  e, 
f.  Each  of  these  will  be  briefly  discussed. 

6.4.3  Complete  Penetration  of  a  Single  Thin  Target. 

Several  extensive  studies  have  been  carried  out 
involving  complete  penetration  of  a  single  thin  target  in 
the  armor  penetration  region  where  the  projectile  is  rela¬ 
tively  undetormed.  0t  prime  concern  are  the  kinetic  energy 
lost  by  the  projectile  during  penetration  and  threshold  velo¬ 
cities  corresponding  to  certain  confidence  levels  that  the 
target  will  be  just  penetrated,  (i„e.  probability  of  penetra¬ 
tion  0,  1  or  0.5,  the  velocity  corresponding  to  the  latter 
being  known  as  the  ballistic  limit),,  This  type  of  study  is 
not  considered  here. 
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A  preliminary  study  of  hole  diameters  and  spray 

particle  distributions  for  single  thin  targets  penetrated 

53 

in  the  high  velocity  region  have  been  reported  at  BRL 

at  11,100  and  16,500  ±t/sec,  Results  of  hole  diameters 

as  a  function  of  target  thickness  are  replotted  in  Fig.  6.27. 

A  statistical  study  o.t  hole  diameters  produced 

in  thin  Mylar  and  Testar  foils  by  cast  iron  particles  of 

less  than  100  micron  size  at  about  36,000  it/sec  has  been 

11 

reported  by  Kienards  and  Gehring 

6.4.4  Thin  Shield 

Three  studies  oi  penetration  in  a  semi- infinite 

target  protected  by  a  thin  shield  located  some  distance 

24 

in  trout  o±  the  target,  have  been  reported,  Olshaker 

studied  lead  targets  and  shields  impacted  by  lead  and  steel 

39 

projectiles  at  about;  6,000  tt/sec.  Funkhauser  studied 

aluminum  targets  and  shields  impacted  by  copper  projectiles 

at  velocities  ranging  from  1,000  to  12,500  tt/sec,  Wallace 
22 

et  al  studied  aluminum  targets  with  a  variety  ox  shield 
and  projectile  materials  at  velocities  up  to  about  17,000 
tt/sec.  Unfortunately  the  data  tor  the  latter  study  was 
not  made  available  tor  comparison  with  the  other  two 
studies.  Some  results  ot  the  first  two  studies  are  plotted 
in  Figs.  6.2b,  to  6.30. 

6 .4.5  Multiple  Thin  Targets 

A  number  ot.  studies  of  penetration  of  multiple 
spaced  targets  have  been  reported.  Some  ot  the  early 
studies  have  been  largely  qualitative.  Ny smith  and  Summers’^ 
tired  glass  spheres  at  aluminum  plates,  varying  the  number  ot 
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plates  and  plate  thickness  such  that  the  total  thickness  of 
the  combined  plates  remained  constant „  They  found  the  thresh¬ 
old  velocity  to  just  penetrate  all  of  the  plates  for  each  con- 
figurations  (Fig.  6o31.)e  Halperson  and  Fuller^  found  the 
total  number  of  1/16  inch  thick  aluminum  plates  spaced  4 
inches  apart  penetrated  by  1/4  inch  steel  spheres  as  a  func¬ 
tion  of  velocity,  (Fig„  6.32)..  Projectile  breakup  is  indica¬ 
ted  by  a  reduction  in  the  number  of  plates  penetrated  as  the 
velocity  increases. 

6.4.6  Thin  Target  Theory 

The  numerical  solution  of  Bjork  (Section  3.2.1) 
is  equally  applicable  to  thin  targets  as  to  qua si- infinite 
targets.  All  that  is  required  is  the  addition  of  a  new 
boundary  condition,  specifying  the  rear  surface  of  the  target. 
While  several  such  cases  have  been  computed,,  published  results 
were  not  received  for  inclusion  in  this  study. 

A  simplified  theory  for  the  angle  of  spray  produced 
by  the  impact  of  a  very  high  velocity  projectile  on  a  thin 
shield  has  been  given  by  Lull 

The  projectile  is  assumed  to  punch  out  a  section 
of  target  plate  of  area  equal  to  the  presented  area  of  the 
projectile,.  Shortly  after  impact,  the  projectile  and  punched 
out  section  of  target  plate  are  assumed  to  have  moved  out  of 
contact  with  the  remaining  target,  and  no  further  interaction 
with  the  remaining  target  is  expected,  At  this  point  the  only 
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part  oi;  the  remaining  target  which  has  been  affected 
by  the  impact  is  that  behind  the  shock  vmve  proceeding 
from  the  edges  of  the  hole-  Any  enlargement  of  the  hole, 
spall  or  lip  formation  must  be  a  result  of  energy  stored 
behind  the  shock ,  For  a  thin  plate  this  energy  is  assumed 
to  be  negligible  compared  to  the  energy  contained  in  the 
projectile  and  punched  out  segment  of  target „ 

The  momentum  balance  then  gives 

rvh  K<  =  (mi  K  (6,11) 


The  energy  balance  gives 
i  K  +  A  £ 


(6 


12) 


where  A  E  is  the  increase  in  internal  energy  of  the  pro¬ 
jectile  and  target  segment „  Thus 


a  e  -- 


2 


i’r  ,  rr^ 


(6.13) 


The  projectile  and  target  segment  are  assumed 
to  shatter.,  and  all  of  the  increase  in  internal  energy  is 
assumed  to  go  Into  acceleration  of  the  material  in  a  radial 


direction.  The  projectile  and  shield  material  are  therefore 
as  sowed  to  be  distributed  throughout:  an  expanding  sphere  of 
radius  R  }  the  center  of  gravity  of  which  is  moving  at  velo¬ 
city  V0  .  The  kinetic:  energy  due  to  radial  expansion  of  an 
element  of  the  sphere  is 


where  is  the  radial  velocity  at  the  surface. 

r  A 

A  E  =  i  I-*-  f  f  r  + cl  r 

o 

Noting  that  the  mass  of  the  sphere  is 

4*  O  3 

-  y  Tr  /-\  jj 

AE  =  lo  ("•  +  •*1)  K* 

Thus  from  Eq.  6.13 


Hence 


(6.14) 


(6.15) 


(6. 16) 


The  leading  hemisphere  of  material  is  moving 
at  higher  velocity  than  the  trailing  hemisphere;,  and  con¬ 
tains  most  of  the  kinetic  energy „  It  would  intersect,  a 
circular  area  on  a  second  target  plate  subtending  a  half 
angle  at  the  shield 

re  t  a  n  -  a  ret  5n.  j  —  — ^  (6.17) 


The  trailing  hemisphere  would  affect  a  larger 
area  of  target.,  but  may  be  expected  to  produce  relatively 
less  damage  since  it  contains  less  energy.  The  circular 
area  intersected  by  the  trailing  hemisphere  subtends  a 
half  angle  9-^  at.  the  shield 


W 

9,  -  5  resin 


n= 

2)  re  S  in  j  -- 
\!  " 


(6. io) 


Those  are  inserted  in  Fig,,  6.30  ter  comparison 
with  experimental,  data. 

6.4.7  Discussio n 


None  of  the  st.udi.es  reported  to  date  are  directly 
comparable.  The  darn,  is  far  too  scanty  to  allow  an  attempted 
correlation.  Some  qualitative  remarks  may  however  by  made. 

Projectile  breakup  is  indicated  by  a  drop  in 
total  penetration  in  a  given  projectile'  target  configuration 
as  the  velocity  i.s  increased..  Similarly  a  relatively  sharp 
drop  in  total  penetration  occurs  as  the  plate  thickness  is 
increased  at  constant  velocity,  indicating  projectile  break¬ 
up.  Thus  conditions  roughly  corresponding  ro  projectile 
breakup  may  be  estimated  for  several  of  the  studies  considered. 
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The  velocity  au  which  projectile  breakup  occurs 
appears  to  be  a  function  of  plate  thickness  as  well  as 
projectile  and  target  material,,  thicker  targets  requiring 
a  lower  velocity  for  projectile  breakup. 

There  appears  to  be  some  correlation  between  the 
diameter  of  the  hole  produced  in  a  thin  plate.,  and  the 
cone  angle  over  which  the  plate  and  projectile  material  is 
ejected  from  the  rear  of  the  plate  (Figs,  6,27  and  6,30), 
both  increasing  with  plate  thickness  at  a  given  velocity,. 

This  again  implies  that  the  projectile  experiences  more 
severe  breakup  as  the  plate  thickness  is  increased,  at 
constant  velocity.  Lull’s  theory  shows  the  same  trend. 

As  the  plate  thickness  is  increased  further  (beyond  about 
half  the  length  of  the  projectile  in  Fig,  6,30),  the  spray 
angle  no  longer  increases.,  and  later  decreases,  Undoubtedly 
the  spray  angle  continues  to  decrease  until  the  plate  is  so 
thick  that  it  will  not  longer  be  completely  penetrated.  In 
Lull’s  theory  the  angle  continues  to  increase  with  plate 
thickness.  The  discrepancy  arises  since  his  assumption 
that  the  target  i s  thin  is  violated. 

The  total  depth  of  penetration  in  a  shielded 
target  or  multiple  target  is  also  a  function  of  pellet 
breakup.  As  the  shield  plate  thickness  increases  (Fig  6„28), 
and  the  spray  particles  are  distributed  over  a  wider  area, 
the  penetration  of  the  prime  target  is  reduced.  At  about 
the  shield  thickness  at  which  the  spray  angle  becomes  con¬ 
stant.,  the  penetration  in  the  prime  target  also  becomes 
constant,  The  Coal  penetration  (shield  plus  prime  target) 
therefore  goes  through  a  minimum,  and  increases  with  further 
increase  in  shield  thickness,,  The  optimum  shield  thickness 
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appears  to  be  about  0.4  times  the  projectile  length,  from 

22 

Fig.  6 o 28 o Wallace'  states  that  the  optimum  In  his  study 

also  lies  at.  about:  the  same  shield  thickness  ratio. 

Total  penetration  is  seen  to  Decrease  with  increase 
in  shield  spacing  (Figo  6. 29)  until  at  some  optimum  spacing 
no  further  reduction  in  penetration  is  achieved.  It  appears 
that  beyond  this  spacing  each  particle  in  the  spray  strikes 
a  different  portion  of  the  prime  target.,  while  at  smaller 
s pacings  several  spray  particles  may  strike  the  same  loca- 
tion  on  the  prime  target.,  leading  to  greater  penetration., 

The  fact:  that  Funkhauset  and  Olshaker  find  different  opti¬ 
mum  spacing®  and  minimum  total  penetration  suggests  that 
the  projectiles  used  in  the  two  studies  suffered  different 
degrees  of  breakup  at  the  particular  velocities  used.. 

These  effects  are  illustrated  in  Fig.  6.33  where 
it  may  be  seen  that  an  increase  in  shield  thickness,  or  an 
increase  in  shield  spacing  lead  to  increased  angle  of  spray 
and  decreased  penetration.. 

Wallace  noted  chat  the  spray  angle.,,  calculated 
on  the  basis  of  spray  diameter  at  the  prime  target,  and  the 
shield  spacing.,  appeared  to  decrease  with  increasing  shield 
spacing v  a. If  other  parameters  being  maintained  constant.. 

This  suggests  that  the  spray  particles  actually  originate 
in  an  area  in  the  projectile  above  the  level  of  the  shield, 
an  effect  also  noted  by  Bjork  in  his  machine  calculations.. 

It,  is  clear  that,  much  mo/.e  experimentation  is 
required,,  before  quantitative  trends  may  he  distinguished, 

j.,  _  """  =  ‘  "  :~  — 
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Ideally *  the  particle  size.,  mass  and  velocity  distributions 
in  the  spray  ejected  from  a  thin  target  should  be  measured 
as  a  function  of  target  thickness  and  velocity*  for  given 
projectile  and  target  materials.  The  effect  of  a  given 
spray  on  a  subsequent  target  should  ideally  be  investigated 
separately.  In  this  way  the  possible  parametric  combina¬ 
tions  to  be  investigated  are  minimised*  and  comparison  with 
possible  theoretical  analyses  facilitated. 
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SECTION  VII 


CONCLUSIONS 


Available  information  relating  to  cratering  and 
penetration  by  projectiles  in  the  velocity  range  where  the 
projectile  suffers  major  deformation  or  breakup  has  been 
surveyed.  Much  of  this  information  is  collected  in  this 
report  for  ready  reference. 

While  a  great  deal  of  experimental  data  has  been 
generated,  most  of  this  data  shows  considerable  scatter,  and 
in  most  cases ,  insufficient  information  about  target  material 
properties  has  been  reported. 

Most  of  the  experimental  work  relates  to  normal 
impact  on  quasi- infinite  targets.  Although  a  very  wide 
variety  of  target  and  projectile  materials  have  been  used, 
projectile  velocities  for  most  of  these  have  been  restricted 
to  the  transition  region.  Almost  no  experimental  data  exists 
above  7  km/sec„  Two  empirical  expressions  have  been  deduced 
which  are  more  generally  applicable  than  previous  expressions 
The  power  law  (Sect,  5.3)  is  restricted  to  the  high  velocity 
region,  defined  in  Section  5,6,  while  the  logarithmic  law 
(Sect.  5o4)  fits  the  experimental  data  in  both  the  transi¬ 
tion  and  high  velocity  regions  for  "ductile"  projectile 
behaviour.  Both  expressions  are  good  fits  to  the  data  when 
constants  appropriate  to  the  projectile  and  target  material 
are  used,  and  are  therefore  useful  for  interpolation.  Extra¬ 
polations  to  higher  velocities  may  be  expected  to  lead  to 
large  errors. 

Only  very  rough  correlations  of  the  constants 


in  the  empirical  expressions  with  elementary  material 
properties  can  be  found,  so  that  predictions  of  penetra- 
tion  for  projectile  and  target  materials  other  than  those 
for  which  sufficient  experimental  data  exists,  may  be  expec¬ 
ted  to  be  subject  to  very  large  errors,.  That  some  correla¬ 
tion  exists  at  all  must  be  ascribed  to  the  fact  that  some 
correlation  exists  between  dynamic  compressibility  and 
strength  properties,  and  elementary  static  properties  of 
the  materials  involved,.  Furthermore,  in  the  case  of  the 
logarithmic  law,  it  may  be  expected  that  a  much  more  complex 
functional,  dependence  on  material  properties  will  be  neces¬ 
sary  in  order  to  accomodate  the  many  complex  transition 
region  strength  effects „  It  appears  to  be  a  hopeless  task 
to  deduce  such  a  functional  relationship  from  the  data,  due 
to  the  large  uncertainties  in  the  data,  ana  the  limited  velo¬ 
city  range  of  the  data  for  most  material  combinations „  An 
adequate  functional  relation  must  be  sought  on  theoretical 
grounds „ 

The  available  theories  applicable  to  high  velocity 
cratering  have  been  found  to  be  inadequate  in  the  experimental 
velocity  range „  Most  of  the  theories  were  intended  to  describe 
cratering  at  much  higher  velocities  and  are  therefore  very 
simplified,  generally  neglecting  strength  effects,  or  intro¬ 
ducing  a  simple,  correction  for  strength  effects.  It.  has  been 
found  that  strength  effects  are  important  at.  much  higher  velo¬ 
cities  than  previously  supposed..  In  particular,  strength 
effects  are  important  in  the  high  velocity  region,  which  is 
of  prime  interest  for  missile  and  space  vehicle  applications. 

A  good  qualitative  understanding  of  physical  phenomena  involved 
in  the.  cratering  process  has  been  gained  It  is  appropriate 
therefore  to  attempt  to  develop  a  more  realistic  theory  at 
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this  time.  The  theory  will  necessarily  have  to  include 
realistic  nonlinear  compressibility,  geometric  and  strength 
effects  in  order  to  provide  a  realistic  analog  to  the  cra¬ 
tering  process. 

For  projectile-target  configuration  other  than 
normal  impact  of  compact  particles  on  quasi- infinite  targets, 
relatively  little  experimental  work  has  been  reported.  In 
general  sufficient  experimental  information  exists  to  dis¬ 
tinguish  qualitative  physical  behaviour,  but:  insufficient 
data  exists  to  attempt  empirical  correlations,,  Theoretical 
work  has  generally  been  limited  to  simple  semi -empirical 
treatments  which  have  shown  limited  success  in  some  cases. 

It  must  be  concluded,  therefore,  that  at  present, 
no  reliable  predictions  of  cratering  and  penetration  can  be 
made  for  velocities  materials  and  configurations  other  than 
those  for  which  experimental  data  exists. 

Further  experimental  work  is  required,  particularly 
for  impact  by  projectiles  of  high  or  low  length/diameter 
ratios,  by  projectiles  at  oblique  incidence,  and  for  impact 
on  thin  targets..  Little  is  to  be  gained  from  further  experi¬ 
mental  work  on  quasi- infinite  targets,  except  at  velocities 
over  6  km/sec,  However,  care  must  be  exercised  in  the  de¬ 
sign  of  experimental  programs,,  Since  additional  independent 
variables  have  been  introduced,  it  will  clearly  be  impossible 
to  carry  out  complete  empirical  surveys  in  which  all  of  the 
experimental  parameters  are  varied  over  the  entire  range. 

This  would  obviously  lead  to  an  experimental  program  of 
prohibitive  magnitude. 
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It  seems  imperative  to  carry  out  theoretical  work, 
so  that  adequate  scaling  laws  may  be  deduced.  Once  scaling 
laws  have  been  established,  intelligent  experimental  pro¬ 
grams  may  be  designed  to  supply  the  necessary  constants,  and 
the  experimental  work  may  be  reduced  to  manageable  propor¬ 
tions  . 


146 


APPENDIX  A 


1  r>7 

Zaid8s  Theory'1' 


The  fluid  in  the  zone  between  the  penetrator  and 
target  is  at  time  t  situated  in  a  cylinder  radius  r  , 
and  height  h  .  At  time  (t  4-  dt)  this  will  be  a  cylinder 
of  radius  (r  +  i.dt)  with  thickness  (h  -  p„dt),  where  p 
is  the  velocity  of  advance  of  the  projectile =  For  the 
conservation  of  mass  of  this  incompressible  fluid 

-i r  rx  h  =  ir  (  r  +  ,dt)(  h,  ~  f>.dt)  (A .  l) 

where  r^  is  the  average  fluid  velocity  over  depth  h  • 

This  reduces  to  the  equation 


from  which,  on  dif ferentiating. 


r 


(A  „  2 ) 


(A. 3) 


An  element  defined  by  r  ,  rde  ,  h  in  the  fluid- 
zone  will  lead  to  the  equation  of  motion 

#  0  h  r  .  J  S> .  Jr-t  -j-  r.  dr.  d6> .  dt 

d  r 


jOt>  ft  +■  dtj.dfl  ft  *-<di),(t~  -f-rdt)  ~  ft  ^  dfi(t) 
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urn 

where  7  -r-*1  is  the  average  of  t?  over  the  depth  of  the 
f  or  0  dr 

fluid  zone,  P  is  the  pressure  within  the  fluids  tu  is  the 
shear  at  the  projectile- fluid  interface,  and  dA  an  element 
of  area  in  the  radial  direction.  Using  the  conservation 
of  mass  equation,  the  relation 


h »  kip)  ~  P  dt  +  f  A* 
y  -  a  +  p 


(A. 4) 


where  q  is  the  velocity  of  penetration  into  the  target, 
the  above  expression  reduces  to 


Jr  K 


oi  r 


+  ni 


-  -  ft h  [*  * 


A 


+ 


1  _£xl 

2  h  J 


(A. 5) 


Integrating  within  the  limits  r  and  R  (7  independent 
of  r). 


r  P,  (r) 


(A.  6) 


Defining  the  ratio  of  the  average  pressure  to  that  at  the 
penetrator  fluid  zone  interface  as 


=  o<  Pu. 


(A,  7} 


the  equation  becomes 

-7  -  --  r /^.^r  +  (Ri-r')('p  *  *  ~p' ) 


s  ’  P  ir) 


(A. 8) 
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The  retarding  force  on  the  penetrator  is 


f/p. 


2  ir  I  i  R.(r) 
J  '• 

0 


Pu  (R)  )  r.  dr 


hence.,  the  equation  of  motion  of  the  penetrator  is 


R  R 


%  P 


2  it  f  f  (  ft  R  ( ”  ^  a  Jl  j  l  A 

*  rt  j  J  J  K-^.rJr  +  ~ -  [r  +  2.  h  J  J  s  0 

^  o  r 


(A .  9 ) 


- 


Zaid  assumed  shear  relationships  of  the  form., 

K  tr,  -  k'-jf  (A.  10) 


and  two  different  fluid-zone  growth  relationships , 


ft 

r 

where 


r 


-  a.. 


(A. 11) 


(A. 12) 


are  the  average  values  of  t-,  ,  r,  over  a  plane  surface 

L  a 

area  of  radius  R  , 
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Three  cases  of  fluid -zone  growth  relationship 
have  been  considered,  n  =  1  ,  2  ,  m  =  1 


Jb 


1  - 


9  h 


II 


(A.  13) 


__9t  D  T 

jr  R  p 


in 


These  lead  to  penetration- fluid  thickness  relations: 

’L  k'ftV/x 


-  h  ♦  Ap-  ( 
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i  (k^} 

°%  n 
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b  k'Ri'/i 


-  R 


A.  14) 


a  -  K  + 


R  a. 


1e 


R  a. 


R  a, 
3 


III 


If  the  steady-state  penetration  dominates,  a  good  approxi¬ 
mation  is  p  ^  pQt  o  Using  this  relation  Case  II  has  a  solution 


(A.  15) 
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This  equation  estimates  that  977,  of  the  steady-state 
fluid- zone  is  established  within  the  time  the  penetrator 
has  traversed  one-half  the  depth  of  that  zone.  Thus,  with 
the  exception  of  the  transient  effects,  a  reasonable  solu¬ 
tion  is  obtained  on  assuming  steady-state  conditions.  In 
the  above  equation  this  would  be 


(A.  16) 


Generalizing  this  to 


L 

f  I 


i  k'R)^ 

[~r  ) 


(A. 17) 


and  substituting  in  the  equation  of  motion  of  the  penetrator, 
results. for  the  approximation 


hf-  3  .)* 

f  V  k'R  1 


(A.  18) 


(which  Zaid  considers  to  be  valid  in  the  range  of  interest) , 
in  a  solution 


Assuming  that  the  fluid  zone  was  of  constant 
thickness,  the  motion  of  the  penetrator  would  be  given  by 
the  equations 


± 

A 


exp 


(A. 20) 
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(A. 20) 


(A.  21) 


A. 22) 


In  order  to  evaluate  these  equations  it  is  necessary 
to  enumerate  k  ,  k'  ,  ,  bg  ,  a  ,  and  in  the  case  of  con¬ 

stant  h.  the  value  of  h  itself.  Values  of  these  constants 
have  not  been  given. 

However,  Za:Ld  does  consider  a  hypothetical  case  in 
order  to  estimate  the  magnitude  of  the  penetration.  For  this 
the  flow  is  assumed  to  be  parabolic  with  a  shear  stress  rela¬ 
tion  at  the  projectile-fluid  zone  interface 


where  (x 


6  r„ 

s-  — 


is  the  viscosity  of  the  fluid. 


The  estimated  ranges 
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of  material  and  other  properties  used  were  as  follows: 

(3  x  10“3)(0.06  x  IQ-3)  ?  <  .06  x  lo"3  slug/ft. sec. 

(water  vapor)  (molten  iron) 

Pp  =  Pt  -  15.2  slug/ft.2 
R  =  L  ■»  1/4  inch 
P  -  1 

1.0  £  R/h  >  10‘3 

For  an  impact  velocity  of  10,000  ft/sec,  the  theory 
which  assumes  a  fluid  zone  of  constant  height,  predicts 
a  crater  depth  anywhere  from  3.8  x  lO-^’  inches  to  4.5 
inches,  depending  on  the  viscosity  and  fluid  zone  thick¬ 
ness  assumed . 
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Opik*s  Meteor It ic  Impact  Theory 


In  solving  the  equations  set  out  in  Section  III; 
it  is  easier  to  introduce  the  dimensional  variables  P  P/2R  ; 


m 


*  in  ' 

p/2R  ,  r  **  r/R  and  $  «*  fp T/(pI  +  P  )'  .  The  problem 

in  c  t  p 


is  now  reduced  to  solving  the  equation 


On  ~  4- 


I  s 


\  •  X 

Prv>  ^ r* 


pr)  l 


for  the  intial  condition  t 


//  -  -UltLl1)  *  4-  J  (b. l) 

f  o/  On2  j 


0  ^  r  •»  1  .  p  -  V/2R  ,  and 
m  ' 


^  ip)  = 


<5  Rm  (°) 
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2.  S 
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I  -t-  2.  <f 


(B.2) 


where 


<P  - 


1  k 


ft 


R  ; 


2_k 

R2 


fa.  t  A  J* 


(B.3) 


The  nonlinear  differential  equation  could  not  be  integrated 
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^  PaTent  c° 

"ot  filme9des  ti,t 


directly.  However,,  the  reverse  differential  equation  obtained 
by  the  transformation. 


dt 


I 


M  y  dt  “  3 

is  easier  to  haridle.  The  transformed  equation  will  be 


(B.4) 


drL 


.  -L/_X_]/_jr±_)  / , .  itlh. £1  +  d>^\ 
f  (  '-V/(  /  f  if  /  [  r„‘  y 


(B.5) 


for  the  initial  conditions 


rm  *  I 


This  has  the  solution 


'i'f1) 


R  I  +  1o 

v  5 
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where 

/(r0  = 


r»  fri  J  [  i 

f  (r-J 
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(B  •  7 ) 
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9(r-y  “ 
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I  *  H 
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(B.8) 
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The  maximum  penetration  is  obtained  when  r  «  p  *  0 
The  radius  of  the  meteorite  at  this  time,  R,  is  obtained 
by  placing  rm  *=  0  ,  hence 


\  2.  £  i 

'  v  ) 

l  R  / 

(l  +  lf ) 

1*  1 

^  Vo  J 

(B.9) 


and  the  maximum  penetration,  p  „  ,  is  obtained  from  r 

max  max 


ZR 
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(B. 10) 


157 


This  Document 

Reproduced  From 
Best  Available  Copy 


APPENDIX  C 


St anyukovitch  Equivalent  Explosive  Analogy 

It  is  assumed  that  a  powerful  shock  wave 
emanates  from  the  point  of  impact  of  a  projectile  onto  a 
target  which  is  capable  of  completely  destroying  the  solid 
media  within  a  certain  volume,  and  converting  it  into  a 
strongly  compressed  gas.  The  pressure  at  the  front  of  the 
shock  will  be  given  by  the  relation 


(C.l) 


where  r  is  the  distance  from  the  point  of  impact,  n  the 
polytropic  index  of  the  expanding  gas;  P  the  pressure  on 
the  front  of  the  shock  wave,  and  ■q  a  parameter  depending 
on  r.  such  that  at  n  «  1,  7/5,  3  or  5,  n  is  equal  to  1/3, 
i/5j  1/6  and  1/7  respectively.  For  strong  condensations 
of  a  gas  at  pressures  above  a  megabar  over  a  wide  range  of 
pressure,  it  is  assumed  that  n  =  3, 

The  Rankine-Hugoniot  relations  across  the  shock  are 
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v;here  e  is  the  internal  energy,  the  target  density 

behind  the  shock,  p  the  target,  density  in  front  of  the 
shock,  U  the  velocity  of  the  shock  wave  front,  and  v  the 
velocity  of  the  gas  behind  the  shock  wave  front „  This  mode 
of  cratering  will  continue  until  the  energy  on  the  shock 
wave  front  is  less  than  £  , 

5  «  e  + .  —  --  V2-  (0.3) 


where  £  is  the  internal  energy  required  to  disintegrate 
the  target  material.  For  the  earth's  crust  Stanyukovitch 
estimates  £  as  approximately  3  cal/gm«  From  these  equ¬ 
ations,  the  radius  r  of  the  hemispherical  crater  formed 
during  this  mode  may  be  found. 
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Now  the  mean  density  of  the  energy  throughout  this  crater  is 
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which  for  n  «  3  (q  =  1/6),  gives  <£>  =  £/2  .  This  energy 

density  for  impact  into  the  earth  is  of  the  order  of  the 
initial  density  of  high  explosives,  therefore  the  entire  mass 
of  the  substance  in  this  crater  is  comparable  to  an  equal 
weight,  of  high  explosive 

Empirically,  it  has  been  found  that  the  radius  of 
a  crater  p  formed  in  the  earth  by  the  detonation  of  a  mass  M 
of  high  explosive  at  the  surface  is  given  by 
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that  is 
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APPENDIX  D 
SHOCK  RELATIONS 
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The  mechanical  shock  relations,  derived  from 
the  laws  of  conservation  of  mass  and  momentum  are 


and 
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where  U  is  the  shock  velocity,  v  the  particle  velocity 
v  the  specific  volume  and  P  the  pressure  respectively, 
and  (  )  refers  to  condition  ahead  of  the  shock,  (  )^  to 
conditions  behind  the  shockc 
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Experiments'^^  have  shown  that  the  relation 

between  shock  velocity  and  particle  velocity  behind  the  shock 
relative  to  the  particle  velocity  ahead  of  the  shock  ai&y  be 
represented  by  a  simple  linear  relation  to  surprising  accuracy 

U  -  v  =  c  +■  5  (V;  -  v0)  (d»3) 

where  c  may  be  identified  as  the  adiabatic  bulk  sonic 
velocity,  and  S  is  related  to  the  change  in  compressibility 
with  pressure o 

Eliminating  v  between  Eq<>  (D,l)  and  (Do 2)  and 
using  Eq0  (D„3).,  a  relation  between  interface  pressure  P 
and  velocity  V  may  be  devised  by  considering  the  shock 
moving  into  the  stationary  undisturbed  target;., 

P'  =  ft,  V*  (ct  +  st  v»)  (D.4) 

where  ih  ~  1/v  represents  density..  Similarly  by  considering 
the  shock  moving  into  the  undisturbed  projectile  material  which 
has  a  velocity  V 

p  *  -  fFo  (v  -  V  j[  *  5P  (V  -  V  *)]  (D.S) 

Equation  (D.4)  and  (D,5)  are  two  simultaneous 
equations  for  interface  pressure  and  velocity. 

Explicit  solutions  for  P  and  V  are  difficult 
to  obtain,  but  the  solution  can  be  obtained  very  simply  by 
plotting  Eq .  (D.4)  and  (D5 )  and  obtaining  the  solution  graphi¬ 
cally,  Since  the  intercept  of  Eq,  (D„5)  on  the  abscissa  is 
the  projectile  velocity  V  ,  it  is  customary  to  draw  the  two 
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curves  on  overlays,  so 
that  they  may  be  positioned 
to  obtain  solutions  for  any 
projectile  velocity. 

For  the  acoustic  approxi¬ 
mation,  i.e.  neglecting  the 
change  of  compressibility  with 
pressure,  3=0  and  the  above 
curves  become  straight  lines . 

An  explicit  solution  may  be 
found 
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The  condition  that  the  shock  in  the  projectile 
moves  below  the  original  target  surface  is  Up  >  0  « 
Applying  Eq.  (D.l),  (D„2)  and  (D..3)  to  the  shocks  in  the 
projectile  and  target,  and  eliminating  the  specific  volume, 
the  condition  reduces  to 


V  > 


(Do  8) 


which  may  be  evaluated  with  the  use  of  Eq.  (D  4)  and  (D.5). 
Again  a  semi-graphical  method  is  useful.  Using  Eq.  (D.5) 


V  >  <zf  <■  s,(v-  V*) 
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and  a  simple  trial  and  error  procedure  employing  the 
curves  previously  used  to  find  P  and  V  may  be  used. 

When  the  projectile  and  target  materials  are 
identical,  V  =  %V  by  symmetry,  and 

V  >  - - —  (I).  10) 

/  - 1  sP 

In  the  acoustic  approximation,  =  0  and 

Eq*  (D.9)  and  (D,10)  reduce  to 

V  >  CP  (Doll) 

as  may  be  expected. 

Interface  velocities  and  pressures  obtained 
graphically  from  P  -  v  plots  are  plotted  in  Fig,.  5,1a 
and  b  as  a  function  of  projectile  velocity..  The  V  vs„  V 
curves  are  very  nearly  straight  lines,  over  the  velocity 
range  shown.  In  fact  they  may  be  approximated  by  straight 
lines  given  by 


above  a  projectile  velocity  of  about  3  km/sec  with  a  maximum 
error  of  about  207,  for  most  materials.  Corresponding  approxi¬ 
mate  values  of  P  may  then  be  obtained  from  Eq„  (D.4) , 

While  the  approximation  is  not  too  accurate,  it  does  afford 
a  simple  means  of  estimating  interface  pressures  and  densities 
rapidly. 
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TABLE  IV 

EMPIRICAL  FITS  TO  THE  PENETRATION  LAW 


,2  1/3 


p/d  s 

p_  2/3 

■U(^) 

( Pt  ) 

(  Ht  * 

Target 

Pt 

(gm/cc) 

Ht  „ 

(kg/nrni^) 

Max.  Vel 
(km/ see) 

Al.  HOOF 

2.79 

•k 

4.37 

Al.  2024  &  24ST 

2.68 

114-120 

3  =  98 

Cu«  (36) 

8.88 

36 

3.33 

Cu.  (65) 

8.88 

65 

6.27 

Pb. 

11.34 

4 

5.00 

Stl,.  1020 

7.69 

111 

3.95 

sti.  1030 

7.69 

123 

3-96 

Stl.  30  RC 

7.74 

302 

3.52 

Cd. 

8  „  66 

23 

5=  01 

Ag. 

10.45 

* 

2  =  54 

Zn„ 

7.13 

45 

2  =  51 

a.  2/3  H  1/3 
k  -  (|£)  (-•%)  (- 

pp  ptv 


0.413 
0.366 
0.419 
0.312 
0.378 
0.371 
0.279 
O': 


0.371 


*• 

+  0.007 
+■  0.009 
±  0.031 
+  0.020 
+  0.010 
+  0.012 
+  0  c  004 
+  0.046 

m» 

in 

+  0. 014 


H  for  these  materials  is  unknown..  The  values  of  k/H^ 
t 


1/3 


f  o> 


these  are :  (i)  HOOF  Al.  target,  (0.287  +  0.008)  x  10  cm/ dyne 

and  (ii)  Ag„  target  (0.277  +  0=024)  x  10^  cm/dyneA^  . 
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OiTO 


TABLE  V 

EMPIRICAL  FITS  TO  TEE  LOGARITHMIC  PENETRATION  LAW 


p/d  »  kL  log  (1  +  } 


Proj . 

Target 

Ht 

Pp/Pt 

Max,  Impact 
Vel. (km/sec) 

\r 

1 

tr 

*2 

Deviation 

(p/d) 

Al,. 

HOOF  Al„ 

* 

1,0 

h  “7 

To  J/ 

0,603 

* 

0.043 

Al. 

2024  Al. 

120 

lo0 

2,00 

0,602 

3,815 

0.091 

Ale 

24ST  Al. 

114 

1,0 

3»98 

0.694 

4.863 

0,049 

Pe. 

2024  Al. 

120 

2,93 

1,65 

0.672 

0.765 

0.099 

Cu. 

2024  Al. 

120 

3,32 

1.86 

1.096 

1.649 

0.146 

Sn  „ 

n 

5°  3 

1„0 

2.27 

Q„421 

9.261 

0.102 

Alo 

1020  Stlo 

ui 

0,342 

2.19 

0.182 

6,089 

0.040 

Fe« 

1.020  St!. 

111 

1.0 

1.71 

0.521 

4,305 

0.106 

Cu  „ 

1020  Sfcl. 

111 

1.13 

1.64 

0,550 

4.874 

0.067 

JPb. 

1020  Sfcl, 

111 

1.46 

1  90 

1.116 

11,639 

0..  049 

Mg» 

1030  Stlo 

123 

0.226 

3»70 

0.270 

17.795 

0.032 

Alo 

1030  Sfcl.. 

123 

0.360 

4.00 

0,287 

4.181 

0.020 

Feo 

1030  Stlo 

123 

lo0 

3,96 

0.501 

3,416 

0.043 

Fe  o 

4140  Stlo 

24b 

1.0 

2,30 

0.359 

1,137 

0.062 

F6o 

30RC  St).. 

3O2 

1.0 

2,26 

0.632 

3.730 

0.020 

CUo 

C-H. 

36 

1.0 

3,33 

0,535 

4.145 

0.029 

Alo 

Ca  0 

65 

Oo  301 

Co  .  9  7 

0.427 

14.695 

0.058 

Cu. 

Cu  0 

65 

1.0 

3 » 06 

0.421 

1.555 

0.080 

Alo 

Pb. 

n 

Y 

0.237 

5.00 

0.454 

63.909 

0.108 

Cu. 

Fb, 

4 

0.786 

2.19 

0,539 

6.847 

0.102 

H  for  this  material  is  unknown,  The  value  of  (k9Hfc)  from 
the  least  squares  fit  is  122.. 5  *  10^  dynes/cm^ 
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The  Penetration  of  Aluminum  Spheres  into  24  ST  Aluminum  Targets 
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1*6  The  Penetration  of  Copper  into  Aluminum  2024  1-3  and  T-4  Targets 


Target:  Aluminum  2014  GRIMMINGER 


The  Penetration  of  Tungsten  Carbide  Spheres  into  Aluminum  2014  Targets 


Fig.  1.9  The  Penetration  of  Alirainum  Projectiles  into  1020  Steel' Targets 


Velocity  (1000  feet  per  sec.) 

1.10  The;  Penetration  of  Aluminum  Projectiles  into  1030  Steel  Targets 


p/d  UoTIJBJQSUOd 


p  /d  uoxT!?a"(auoci 


The  Penetration  of  Lead  Projectiles  into  Lead  Targets 


Velocity  (1000  feet  per  sec.) 

ig.  1.25  The  Penetration  of  Tungsten  Carbide  Projectiles  into  Lead  Targets 


rinell 


0  20  40  60  80  100  120  140  160  180  200 

Yield  Stress  (1000  Ibs/ing^) 

Fig.  4.1  Brine 11  Hardness  versus  Yield  Strength  for 

a  Number  of  Representative  Target  Materials 
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Fig.  5.2  Craters  in  bead  Targets  formed  by  1/8"  Spheres 
at  Approximately  700  £t/see 
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Projectile 


Fig,  5.2  continued  Cratci'S  in  Lead  Targets  formed  by 
1/8"  Spheres  at  Approximately  1500  ft/sec 
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5.2  continued  Craters  in  Lead  Targets 
1/8"  Spheres  at  Approximately  2000 
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Lead 

Proj eetile 


Copper 

Projectile 


Fig*  5.2  continued  Craters  in  Lead  Targets  formed  by 
1/8"  Spheres  at  Approximately  3000  ft/sec 
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ProiocH  1  e 
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Fig.  5.2 


I  MTIIITMl 


ncluded  Craters  in  Lead  Targets  formed  by 
8"  Spheres  at  Approximately  5000  ft/sec 

2T1 


'  ’!! 


ll"  ■’nbl'MllM  fcp 
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Fig.  5.5  Illustration  of  Cratering  in  the  High  Velocity  Region. 

Copper  Target  Struck  by  Copper  Sphere.  The  Projectile 
Material  has  been  Lifted  to  Illustrate  the  Plating 
Effect.  (Ref.  30) 
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Unloosed  I of  PuljUfflllnn  by 
R  A|  |  IS  TICS  MF  *»f  ARCH  I  AllodArORY 
AMf  UOFtN  PnOVINU  (WOUND.  MAHfl  AMO 


High  Speed  Framing  Camera  Pictures  of  Penetration 


of  a  Lucite  Target  by  a  Steel  Projectile  (from  Ref.  7  ) 
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c.  -  165  1-1  SEC  d.  -  ?.80  ll  SEC 

R  n  loot  Ad  lor  Publication  by 
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Fig.  5.8  Flash  X-Radiograph  of  the  Spray  Issuing  from  a 
Crater  in  the  Process  of  Formation  in  a  Lead 
Target  (from  Ref.  7) 
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INITIAL  CONDITIONS 


P  =  e  =  O  EVERYWHERE 
p  =  p0  EVERYWHERE 
i  =  o 

VELOCITY  =  0  IN  TARGET 

VELOCITY  *  V  IN  Y  DIRECTION  IN  PROJECTILE 
PROJECTILE  HAS  JUST  CONTACTED  TARGET 

Fig.  5.9  Numerical  Solution  of  the  Penetration  of  an  Iron 

Projectile  into  an  Iron  Target,  Neglecting  Material 
Strength  (from  Ref.  128  ) 
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PRESSURE  CONTOURS  AND  VELOCITY  FIELD 


a  \  v  m  / 


Fig.  5.9  Continued,  (Mote  Changing  Scale) 
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PRESSURE  CONTOURS  AND  VELOCITY  FIELD 


Fig.  5.9  Continued ,  (Note  Changing  Scale) 
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VELOCITY  AND  PRESSURE  CONTOURS 


Fig.  5.9  Concluded,  (Note  Changing  Scale) 
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Fig.  5.10  Numerical  Solution  for  the  Penetration  of  a 
Meteorite  into  Tuff  (from  Ref.  152) 
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Fig.  5.10  Continued  (Note  Changing  Scale) 
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Fig.  5.10  Continued  (Note  Changing  Scale) 
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Fig.  5.10  Continued  (Note  Changing  Scale) 
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Fig.  S.10  Continued  (Note  Changing  Scale) 
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Fig.  5*10 


Concluded  (Note  Changing  Scale) 
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Velocity  (1000  Lfc.  per  sec.) 


Fig.  5.13  Penetration  in  Copper  Targets 
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Fig.  5.14  Penetration  in  Lead  Targets 
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ig,  5.29  Depth- Diameter  Ratios  for  Craters  in  Aluminum  Targets 
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Fig.  5.33  Energy-Volume  Ratios  for  Aluminum  Targets 


Kneivy- Volume  Ratio  Dependence  on  Density 
Rat  io 
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Fig,  5-33  Initial  Interface  Velocity  Dependence  en  Projectile  Density 
at  Constant  Projectile  Mass  and  Velocity 
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Fig.  6.1  The  Penetration  of  Stainless  Steel  Spheres  (lOOp 

and  150(i  diameter)  into  Various  Targets 
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rig.  6./  Craters  formed  in  Load  Targets  by  Steel 
Projectiles  at  Oblique  Angles  of  Impact 
(from  Ref.  7) 
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Fig.  6.21  Crater  formed  by  a  Slender  Unyawed  Ductile 
Rod  at  Low  Velocity 
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Fig.  6.22  Crater  formed  by  a  Slender  Rod  Yawed  at  10° 
to  Direction  of  Motion 
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Fig.  6.25  Photographs  of  Projectile  Penetrating  a  Thin  Target. 

The  target  is  at  the  Loft,  and  Projectile  Motion  is 
from  Right  to  Left.  (0.22  inch  Aluminum  Sphere 
Impacting  1/8  inch  Aluminum  Target) 
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Fig.  6. 28  Total  J’ent-U  atioii  of  a  Shielded  Target  in  Terms  of 
Penetration  in  an  Unshielded  Ta  rget  an  a  Function 
of  Shield  Thickness 
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